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THE  WEST  SPITSBERGEN  CURRENT: 
TRANSPORT,  FORCING,  AND  VARIABILITY 

by  Dennis  Joe  Hanzlick 


Chairman  of  Supervisory  Committee:  Knut  Aagaard 

Research  Professor 
School  of  Oceanography 

The  West  Spitsbergen  Current  provides  the  major  oceanic  influx  of 
heat  and  mass  into  the  Arctic  Ocean.  Recent  measurements  show  that  the 
flow  has  significant  lateral  structure  and  exhibits  large  variability 
over  a  broad  range  of  time  scales.  Wind  forcing  appears  to  be  particu¬ 
larly  Important. 

During  1976-77,  the  mean  volume  transport  was  5.6  Sv  northward  and 
ranged  from  1.4  Sv  southward  to  11.9  Sv  northward.  The  largest  trans¬ 
ports  occurred  In  late  fall  and  winter,  and  the  minimum  was  In  March. 
The  corresponding  mean  heat  transport  was  46  TW  (one  tera-watt  (TW)  * 
10*2W),  with  a  range  of  8  TW  southward  to  122  TW  northward  and  a 
seasonal  pattern  similar  to  that  of  the  volume  transport.  The  along¬ 
shore  and  cross-isobath  turbulent  heat  fluxes  were  respectively  1  and  2 
orders  of  magnitude  smaller  than  the  alongshore  advectlve  mean  flux. 
Three-year  long  current  records  showed  consistent  seasonal  fluctuations 
but  little  Interannual  variability  between  1976  and  79.  However,  typi¬ 
cal  velocities  were  only  about  one-half  as  large  as  those  measured  In 
1971-72,  Indicating  that  significant  Interannual  changes  can  occur. 

Current  records  Indicate  lateral  structure  over  the  upper  slope 
west  of  Spitsbergen  with  cross-isobath  length  scales  of  10-20  km. 

The  regional  wind  stress  curl  can  account  for  the  northward  volume 
transport  of  the  West  Spitsbergen  Current,  and  the  curl  and  flow  are 
related  at  the  annual  time  scale.  The  oceanic  response  can  be  In  the 
Sverdrup  mode  only  If  forcing  Is  at  time  scales  much  longer  than  a 
month,  and  this  response  will  lag  the  forcing  by  about  six  months. 


•  Numerical  model  results  indicate  that,  while  barotropic  instabil¬ 
ity  is  not  an  important  factor  in  the  variability  of  the  flow,  baro- 
clinic  instability  may  contribute  to  the  variability  at  time  and  length 
scales  of  3-4.5  days  and  30-50  km. 

0  The  dispersion  characteristics  of  the  second  mode  coastal  trapped 

wave  which  could  be  supported  by  the  stratification  and  topography  west 
of  Spitsbergen  include  corresponding  wavelength  and  period  ranges  of 
260-1000  km  and  4-14  days.  The  associated  velocity  structure  is  essen¬ 
tially  barotropic  and  contains  a  node  over  the  slope.  These  features 

®  are  consistent  with  observations.  Both  the  first  and  third  trapped 

wave  modes  are  not  important  in  explaining  variability  in  the  range  of 
a  few  days  to  a  few  weeks. 

Meandering  of  the  banded  flow  may  contribute  to  the  short  period 

•  variability,  and  it  Is  also  likely  that  the  northeastward  passage  of 
cyclones  across  the  region.  In  conjunction  with  bottom  topography  and 
the  deep  extent  of  atmospheric  forcing,  contributes  to  the  variability 
of  the  flow  at  time  scales  corresponding  to  the  frequency  of  their 

A  passage. 
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1.  Introduction 


\ 

The  West  Spitsbergen  Current  (hereafter  referred  to  as  WSC)  is 
commonly  considered  to  be  the  northern  extension  of  the  Norwegian 
Atlantic  Current,  and  it  forms  one  of  the  branches  of  the  circulation 
pattern  of  the  Greenland-Norwegian  sea  system.  It  flows  northward 
along  the  continental  margin  west  of  Spitsbergen,  entering  the  Arctic 
basin  on  the  east  side  of  the  Greenland-Spitsbergen  passage  (Fram 
Strait).  Previous  work  (cfv  Aagaard  and  Greisman,  1975)  has  shown  that 
it  provides  the  major  component  in  the  mass  and  heat  balances  of  the 
Polar  Basin. 

Beginning  in  summer  1976,  moored  current  and  temperature  recorders 
have  been  deployed  in  the  WSC  at  79°N,  and  the  data  records  from  1976- 
79  and  1980-81  provide  the  basis  for  this  work.  The  data  have  been 
used  to  estimate  the  volume  and  heat  transports  by  the  current  and  to 
define  the  structure  and  variability  of  the  flow.  In  addition  to  the 
presentation  of  basic  results,  there  is  a  discussion  of  the  dynamics  of 
the  current  with  emphasis  on  explaining  the  length  and  time  scales  of 
the  observed  variability  and  on  ascertaining  the  nature  of  the  forces 
driving  the  flow. 


1.1  Description  of  the  Green 1 and -Norweg i an  sea  region 


1.1.1  Bathymetry 

Figure  1  shows  the  geographic  setting  and  place  names  for  the 
Greenland-Norwegian  sea  region.  Deep  basins  separated  by  submarine 
ridges  are  the  dominant  features  of  the  bottom  topography  of  the  region 
(Fig.  2).  Mohn  Rise,  a  section  of  the  mid-Atlantic  ridge,  separates 
the  Greenland  Sea  from  the  Norwegian  Sea  to  the  southeast.  The  sill 
depth  of  Nohn  Rise  Is  about  2500  m.  The  Greenland  Sea  itself  is 
divided  by  the  Greenland  Fracture  Zone  into  a  northern  basin  ranging  to 
3100  m  deep  and  a  southern  basin  averaging  about  3600  m  deep.  Glacia¬ 
ted  continental  shelves  form  the  lateral  boundaries  of  the  region.  Off 
Greenland,  the  shelf  break  generally  lies  at  200-400  m,  while  on  the 
eastern  side  of  the  basin  the  shelf  break  is  typically  shallower,  100- 
300  m.  The  slope  region,  subtended  approximately  by  the  2500  m  and  300 
m  Isobaths,  narrows  from  a  width  of  about  350  km  off  northern  Norway 


Figure  1.  Geographic  setting  and  place  names  for  the  Greenland-Nor 
weglan  sea  region,  (adapted  from  Trangeled,  1974) 
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down  to  50-60  km  west  of  Spitsbergen.  This  represents  a  seven-fold 
increase  In  the  average  bottom  slope. 

Figure  3  shows  the  depth  profile  west  of  Spitsbergen  at  79°N,  and 
the  mooring  locations  where  the  current  records  which  form  the  basis  of 
this  investigation  were  obtained.  The  shallow  shelf  extends  from  the 
Spitsbergen  coast  westward  to  about  8°45'E  where  the  depth  increases 
sharply  from  about  200  m  to  nearly  1100  m  at  8°E,  a  very  steep  slope  of 
about  6  x  10"2.  Between  8°E  and  6°45*E  the  depth  increases  only 
slightly  to  about  1300  m,  and  the  bottom  slope  is  a  magnitude  less  than 
that  of  the  segment  nearer  shore.  West  of  6°45'E,  the  bottom  slope 
increases  again,  to  about  2  x  10"2,  as  the  seabed  descends  to  a  depth 
of  roughly  3000  m  at  the  foot  of  the  continental  slope.  For  the  most 
part,  the  isobaths  west  of  Spitsbergen  run  north-south.  In  the  mooring 
area  at  79°N  they  are  oriented  along  352°  true. 

1.1.2  Hydrographic  features 

The  primary  focus  of  this  study  is  on  the  WSC,  and  therefore  I 
only  note  the  salient  features  of  the  entire  Sreenland-Horwegian  sea 
region  as  they  bear  on  the  character  of  the  WSC.  In  the  references  can 
be  found  much  more  detail  of  the  hydrography  and  circulation  of  the 
larger  area. 

While  the  WSC  Is  a  major  component  of  the  Greenland  Sea  circula¬ 
tion,  It  has  long  been  recognized  that  its  ultimate  source  is  much 
farther  south.  Consequently,  the  water  properties  of  the  WSC  differ 
from  those  of  the  rest  of  the  Greenland  sea  region. 

Several  earlier  references  offer  excellent  descriptions  of  the 
hydrography  of  the  Greenland  and  Norwegian  seas  as  well  as  of  the 
Barents  Sea.  Among  them  Is  the  classic  work  of  Helland-Hansen  and 
Nansen  (1909).  Using  observations  dating  from  about  1900,  they  pro¬ 
duced  the  first  comprehensive  treatise  on  the  oceanographic  regime  of 
the  area.  Much  later,  Klllerlch  (1945)  synthesized  data  collected 
between  1891  and  1933  Into  a  cohesive  discussion  of  the  hydrographic 
features  and  circulation.  Lee  (1963)  reviewed  the  hydrography  of  the 
Norwegian,  Greenland  and  Barents  seas,  summarizing  the  work  up  to  that 
time.  His  discussion  centered  on  the  water  masses  and  currents,  and 
included  fluctuations  in  the  hydrographic  regime  and  variations  In 
volume  transport.  Trangeled  (1974)  provided  an  exhaustive  bibliography 
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Bottom  profile  and  mooring  locations  west  of  Spitsbergen 
along  79°N.  The  letters  A,  B,  C,  and  D  Identify  the  re¬ 
spective  mooring  locations. 


and  a  brief  summary  of  the  1870-1970  literature  on  the  oceanography  of 
the  Norwegian,  Greenland,  and  Barents  seas  as  well  as  the  waters  around 
Iceland.  While  their  primary  emphasis  was  on  the  North  Polar  Basin, 
Coachman  and  Aagaard  (1974)  also  reviewed  the  oceanography  of  the 
Greenland  Sea.  Their  discussion  covered  the  major  water  masses,  fea¬ 
tures  of  the  circulation,  and  both  wind  and  thermohaline  forcing  of  the 
motions.  Of  recent  works  dealing  with  the  hydrography  of  the  Greenland 
Sea,  probably  the  most  comprehensive  is  that  of  Carmack  (1972),  upon 
which  the  following  summary  draws  heavily.  Additionally,  Dietrich's 
(1969)  atlas,  containing  hydrographic  sections  obtained  in  the  Green¬ 
land  and  Norwegian  seas  during  IGY  investigations  in  1958,  provides  a 
fine  visual  complement  to  these  other  sources.  In  another  context. 
Swift  (1980)  also  presented  a  recent  overview  of  the  hydrography  of  the 
Greenland  and  Norwegian  seas. 

Carmack  recognized  three  basic  water  masses  in  the  Greenland  Sea: 
Alantic  water.  Polar  water,  and  Deep  water.  He  assumed  all  Greenland 
Sea  water  masses  either  to  be  mixtures  of  these  or  to  be  results  of 
local  modifications  due  to  surface  processes.  He  identified  the  deep 
water  as  typically  having  temperatures  less  than  0C  and  salinities  of 
34.85-34.95.  This  deep  water  mass,  delineated  by  the  OC  isotherm, 
forms  a  huge  domed  feature  centered  on  the  prime  meridian  at  75°N  in 
the  Greenland  Sea.  The  WSC,  both  where  it  flows  northward  and  where  a 
portion  of  It  curves  westward  near  77°-80°N,  depresses  the  OC  surface 
along  the  eastern  and  northern  perimeters  of  the  Greenland  Sea,  forming 
the  eastern  and  northern  boundaries  of  the  deep  water. 

Together,  the  Atlantic  water  (T>0C,  S  >34.90)  and  polar  water  (S 
<34.60)  constitute  the  upper  source  waters  of  the  Greenland  Sea.  The 
upper  waters  may  In  turn  be  categorized  as  belonging  to  one  of  four 
regions  according  to  the  respective  current  systems.  These  are  1)  the 
East  Greenland  Current  area,  2)  the  Greenland  gyre  area,  3)  the  Nor¬ 
wegian  Current  area,  and  4)  the  WSC  area,  which  Is  the  primary  concern 
In  this  paper. 

The  East  Greenland  Current  flows  southward  mainly  over  the  conti¬ 
nental  margin  off  the  east  coast  of  Greenland.  It  embodies  two  water 
masses:  1)  polar  water,  whose  cold  temperature  (<0C)  and  relative 
freshness  (30.0*34.0)  reflect  Its  polar  origin,  and  2)  Atlantic  water 
which  exhibits  high  temperature  and  salinity  with  maxima  at  the  core  of 


>2.0C  and>  35.0,  respectively.  The  layer  of  Polar  water  extends  down¬ 
ward  from  the  surface  to  about  150  m.  The  Atlantic  layer  occurs 
between  150-800  m,  the  result  of  subsurface  contributions  of  water  from 
both  the  WSC  and  the  Polar  Basin. 

The  upper  water  of  the  Greenland  gyre  forms  a  thin  layer  (50-200 
m),  and  according  to  Kiilerlch  (1945)  Is  composed  of  a  mixture  of  Polar 
water  and  Atlantic  water.  Carmack. and  Aagaard  (1973)  characterized  the 
upper  water  of  the  gyre  by  its  salinity,  between  34.60  and  34.90.  Just 
prior  to  winter,  the  overall  density  of  the  upper  water  near  the  center 
of  the  gyre  is  high,  while  the  vertical  stability  is  fairly  low.  Car¬ 
mack  and  Aagaard  viewed  the  increase  in  density  and  the  subsequent 
sinking  of  this  water  due  to  winter  cooling  as  an  important  part  of  the 
deep  water  formation  process.  Carmack  (1972)  pointed  out  that  near  the 
northern  and  southern  edges  of  the  gyre,  T-S  correlations  show  a  con¬ 
sistent  salinity  maximum,  providing  evidence  of  the  zonal  movement  of 
Atlantic  water. 

The  West  Spitsbergen  Current  has  its  Immediate  origins  in  the 
Norwegian  Atlantic  Current  which  flows  northward  into  the  Norwegian  Sea 
basin.  Extending  to  about  800  m,  it  is  characterized  by  a  monotonic 
decrease  of  both  temperature  and  salinity  with  depth.  Polar  water  is 
absent  as  a  component  of  the  upper  water  in  this  region.  Because  water 
temperatures  remain  warm,  4-5C  even  in  winter,  the  eastern  Greenland- 
Norwegian  sea  region  Is  ice  free  year  round.  This  warm,  saline  water 
forms  the  major  component  of  surface  waters  In  the  Norwegian  Sea. 
While  the  Norwegian  Atlantic  Current  Is  commonly  considered  to  be  an 
extension  of  the  North  Atlantic  Current,  Reid  (1979)  has  proposed  that 
outflow  from  the  Mediterranean  may  add  a  major  subsurface  contribution. 

1.1.3  General  circulation 

The  general  pattern  of  water  movement  In  the  Greenland -Norwegian 
sea  basin  has  been  known  for  hundreds  of  years.  Early  in  this  century, 
Helland-Hansen  and  Nansen  (1909)  charted  the  basic  circulation  using 
modem  hydrographic  methods.  While  there  have  been  some  modifications 
and  refinements  In  their  scheme  since  then,  there  has  not  been  any 
drastic  change  In  the  picture.  Figure  4  shows  the  prominent  features 
of  the  circulation  In  the  area.  The  southward  flowing  East  Greenland 
Current  and  the  northward  flowing  WSC  form  the  western  and  eastern 
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Prominent  circulatory  features  In  the  Greenland-Norweglan 
sea  region  (from  Trangeled,  1974). 
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boundaries  of  the  cyclonic  Greenland  gyre.  Zonal  flows  close  the  gyral 
circulation  on  the  northern  and  southern  edges.  Earlier  conceptualiza¬ 
tions  Indicated  a  major  cyclonic  gyre  in  the  Norwegian  Sea,  while  more 
recent  work  (Fuglister,  1954;  Alekseev  and  Istoshin,  1960;  Metcalf, 
1960;  Aagaard,  1970)  has  shown  that  the  northward  flow  in  the  southern 
and  eastern  portions  of  the  Norwegian  Sea  is  much  broader  than  depicted 
earlier.  This  implies  that  the  Norwegian  Atlantic  current  extends 
farther  west  than  thought,  and  that  the  gyral  feature  is  smaller  and 
displaced  farther  to  the  west. 

1.1.4  Forcing  of  the  circulation 

There  have  been  several  attempts  to  determine  the  mechanisms 
responsible  for  the  general  circulation.  Krauss  (1955)  argued  that  the 
variability  of  the  winds  In  the  area  precluded  the  possibility  of  their 
producing  the  observed  circulation.  Instead,  he  proposed  a  thermoha¬ 
line  forcing  mechanism  which  relied  upon  coastal  freshwater  runoff. 
This  Implies,  for  Instance,  that  the  Norwegian  Atlantic  Current  is  a 
secondary  phenomenon  primarily  driven  by  the  coastal  current.  Some 
support  for  the  Idea  may  be  found  in  variations  of  the  coastal  currents 
around  Iceland,  which  Stefansson  (1962)  felt  were  likely  maintained  by 
spring  and  summer  freshwater  runoff.  Coachman  and  Aagaard  (1974) 
pointed  out,  however,  that  rainfall  of  about  100  m  per  year  ("...high, 
even  for  western  Norway")  would  be  needed  to  drive  even  a  modest  circu¬ 
lation  of  the  type  envisioned  by  Krauss.  Worthington  (1970),  consid¬ 
ering  the  Norwegian  Sea  as  a  mediterranean  basin,  suggested  that 
replacing  the  deep,  dense  outflow  to  the  south  from  the  Norwegian  Sea 
necessitated  the  Inflow  of  surface  waters  from  the  North  Atlantic. 
Thermohaline  effects  probably  do  contribute  to  the  circulation,  but  it 
may  be  difficult  to  ascertain  their  exact  role. 

There  Is  also  evidence  that  wind-forcing  does  play  a  major  role  in 
driving  the  circulation,  contrary  to  Krauss*  (1955)  argument. 
Aagaard's  (1970)  calculated  annual  mean  circulation  driven  by  wind 
stress  curl  was  qualitatively  similar  to  the  observed  circulation.  He 
found  that  Sverdrup  dynamics  could  account  for  transports  of  the  magni¬ 
tudes  seen  In  the  Greenland  and  Norwegian  seas;  specifically,  his 
calculations  yielded  a  southward  return  transport  of  35  Sv  In  the  East 
Greenland  Current,  In  agreement  with  the  transport  Inferred  from  direct 
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current  measurements.  Also,  his  calculated  wind-driven  transports  of 
Atlantic  water  into  the  Polar  Basin  through  the  Green land-Spitsbergen 
passage  agreed  well  with  transports  calculated  from  dynamic  sections. 
While  he  termed  these  agreements  fortuitous,  his  interpretation  was 
that  the  calculated  annual  mean  wind-driven  transports  were  of  realis¬ 
tic  magnitudes.  He  also  noted  that  the  wind-driven  circulation  is 
largely  internal  to  the  Greenland  and  Norwegian  seas. 

Semtner  (1976)  studied  the  long-term  circulation  in  the  Arctic 
Ocean  and  Greenland  Sea  using  a  three-dimensional,  atmospherically 
forced  numerical  model.  He  used  the  same  wind  field  over  the  Greenland 
Sea  that  Aagaard  (1970)  used  for  the  mean  annual  case.  His  simulation 
produced  a  northward  flowing  Norwegian  Current,  an  intense  southward 
flowing  East  Greenland  Current,  a  southward  return  flow  of  Atlantic 
water  beneath  the  surface  East  Greenland  flow,  and  a  large  cyclonic 
gyre  over  the  Greenland  and  Norwegian  basins.  All  these  features 
qualitatively  resemble  observed  features  of  the  circulation  in  the 
Greenland  and  Norwegian  seas  (Fig.  4). 

In  connection  with  observed  variability  in  the  WSC,  Greisman  and 
Aagaard  (1979)  proposed  that  wind  stress  curl  acting  over  the  Green- 
land-Norwegian  sea  region  drives  a  northward  Sverdrup  flow  which  is 
ultimately  concentrated  by  topographic  effects,  and  is  Identifiable  as 
the  WSC.  Thus,  they  contended  that  the  WSC  and  its  fluctuations  are  not 
driven  locally,  but  by  wind  forcing  over  the  Green land-Norwegi an  sea. 

1.2  Description  of  the  West  Spitsbergen  Current 

1.2.1  Hydrographic  features 

The  most  obvious  characteristics  of  the  WSC  are  its  relatively 
high  temperature  and  salinity,  reflecting  Its  origin  far  to  the  south. 
Water  from  the  North  Atlantic  enters  the  southeastern  Norwegian  Sea 
through  the  Faroe-Shetland  channel,  flowing  northward  along  the  conti¬ 
nental  margin  as  the  Norwegian  Atlantic  Current.  Mosby  (1970)  assigned 
typical  temperature  and  salinity  values  of  about  9C  and  35.3.  to  the 
Atlantic  water  core  In  the  southeastern  Norwegian  Sea.  Farther  north, 
near  the  southern  tip  of  Spitsbergen,  a  contribution  from  the  East 
Spitsbergen  Current  slightly  freshens  the  upper  layer  of  the  flow  and 
strengthens  Its  shallow  (~50  m)  halocllne.  Continuing  northward,  the 
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flow  loses  both  heat  and  salt.  Even  so,  summer  temperatures  in  the 
upper  layers  reach  6C  and  rarely  fall  below  OC  in  the  winter. 

Several  of  the  salient  features  are  apparent  in  Fig.  5,  which 
shows  contours  of  temperature  and  salinity  in  vertical  sections  along 
79°N  in  late  summer  1968  and  1969.  Notably,  in  1969  a  core  region  of 
high  salinity  (S  >35.2)  and  high  temperature  (T>  4C)  lay  just  seaward 
of  the  shelf  break  at  about  200  m  depth  at  7°30'E.  Within  the  halo- 
cline  in  the  top  100  m,  the  salinity  increased  from  32  to  35,  while  the 
temperature  increased  from  about  2C  to  4C.  The  cooler,  fresher  water 
near  the  surface  was  a  contribution  of  Polar  water  from  the  East  Spits¬ 
bergen  Current.  In  1968,  the  center  of  the  salinity  core  (S>  35.1)  lay 
at  about  the  same  position  as  in  1969.  However,  in  1968  the  halocline 
region  was  only  about  half  as  thick,  and  the  salinity  step  was  only 
about  1.  Clearly,  a  marked  influence  of  the  East  Spitsbergen  Current 
present  in  1969  was  absent  in  1968  at  the  same  location  and  time  of 
year.  Figure  6  (from  Grelsman,  1976;  based  on  hydrographic  data 
obtained  between  1954-72)  shows  a  T-S  envelope  for  the  WSC,  demon¬ 
strating  the  high  temperature  and  salinity. 

Vertical  profiles  of  temperature,  salinity  and  sigma-t  obtained  in 
early  September,  1977  (Fig.  7)  show  the  marked  variation  in  properties 
in  the  upper  50  m.  The  homogeneity  in  the  upper  25  m  gives  way  to  a 
sharp  gradient  region  between  25  m  and  35  m  in  which  the  temperature 
decreases  about  1.5C  and  the  salinity  Increases  about  0.2.  Below  50 
the  water  is  nearly  homohaline,  while  the  temperature  decreases  almost 
linearly.  A  sharp  pycnocllne  lies  at  30  m,  and  the  increase  in  density 
in  the  top  50  m  Is  roughly  1.5  times  greater  than  between  50  and  1000 
m.  Below  50  m,  the  density  Increases  by  about  0.25  sigma-t  units 
almost  linearly  down  to  1000  m,  and  the  static  stability,  E,  is  low, 
with  values  in  the  range  2-3  x  10“6  m"1.  The  extreme  changes  In 
vertical  structure  In  the  WSC  are  therefore  confined  for  the  most  part 
to  the  upper  100  m. 


1.2.2  Previous  estimates  of  volume  and  heat  transports 

Grelsman  (1976),  summarized  the  results  of  earlier  investigations 
of  the  WSC.  Of  necessity  my  own  review  contains  some  of  the  same 
information. 
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In  a  study  of  the  volume  transport  of  the  WSC,  Hill  and  Lee  (1957) 
and  Lee  (1961a)  utilized  hydrographic  observations  beginning  in 
1949.  The  data  were  from  a  section  extending  west  of  Bear  Island  about 
80  km  across  the  eastern  branch  of  the  current.  Referenced  to  a  level 
of  750  db,  their  calculated  transport  above  400  m  depth  varied  between 
3.2  Sv  to  the  north  and  0.2  Sv  to  the  south.  These  values  are  very 
likely  low  because:  1)  contributions  below  400  m  were  neglected,  2) 
their  survey  did  not  extend  westward  across  the  entire  stream,  and  3) 
their  calculations  did  not  include  the  barotropic  contribution.  They 
determined  a  semi-annual  cycle  with  maxima  in  January-February  and  July 
and  minima  in  April  and  October.  They  also  found  a  correlation 
between  the  transport  above  400  m  and  the  strength  of  the  southerly 
wind  component  during  the  preceding  10  days. 

Kislyakov  (1960)  published  transport  values  across  74°30'N  varying 
between  6.7  Sv  and  0.6  Sv  to  the  north,  based  on  a  reference  level  of 
1000  db  from  data  obtained  between  1954-59.  He  also  established  a 
correlation  between  the  volume  transport  above  200  m  and  the  longshore 
(north-south)  atmospheric  pressure  gradient.  Both  Kislyakov  and  Hill 
and  Lee  found  a  residual  northward  transport  when  their  respective 
forcing  mechanisms  were  zero. 

Using  13  years  of  hydrographic  measurements,  Timofeyev  (19$51 
reported  a  maximum  transport  across  78°N  of  about  3  Sv  fn  December- 
January  and  a  minimum  of  about  2  Sv  in  May,  relative  to  1000  decibars. 
Mandel  (1976)  estimated  a  mean  volume  transport  of  Atlantic  water  in 
the  depth  interval  0-500  m  of  4.3  Sv.  with  a  maximum  of  5.8  Sv  in 
November  and  a  minimum  of  2.5  Sv  in  April.  He  estimated  this  to 
represent  about  75%  of  the  transport  to  a  reference  depth  of  1000  m. 

•His  calculations  were  based  on  a  linear  relation  between  the  water 
level  at  Barentsburg  (~78°N,  12.5°E)  and  the  volume  flow  of  the  WSC. 

Greisman  (1976)  was  the  first  to  estimate  transports  using  long¬ 
term  direct  current  records  obtained  west  of  Spitsbergen.  From  current 
meter  data  obtained  in  1971-72  on  two  moorings  west  of  Spitsbergen  at 
79°N,  he  estimated  a  mean  northward  transport  of  about  7  Sv.  His  method 
Involved  drawing  Isotachs  of  northward  velocity  In  the  cross-section 
during  11  weeks  when  all  current  meters  were  operating,  computing  the 
regression  relation  between  the  Integrated  transport  and  the  velocity 
at  one  of  the  current  meters  over  the  slope  west  of  Spitsbergen  for  the 
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same  period,  and  then  using  that  regression  relation  along  with  the 
entire  single  current  record  to  estimate  a  monthly  time  series  of 
transport  during  1971-72.  He  identified  a  maximum  of  10.7  Sv  north¬ 
ward  in  September  and  a  minimum  of  about  3  Sv  in  February.  Using 
hydrographic  data  from  two  cruises  in  1971  and  1972,  he  estimated  the 
baroclinic  component  of  the  transport  at  about  3.5  Sv,  i.e.,  one-half 
the  total.  One  reason  for  the  disparity  between  Greisman's  transport 
figures  and  the  others  above  is  that  the  direct  current  measurements 
allowed  an  estimate  of  the  total  transport,  while  the  reliance  upon 
dynamic  calculations  limited  earlier  investigators  to  the  baroclinic 
portion  only. 

Greisman  (1976)  also  used  temperature  records  concurrent  with  the 
velocity  records  to  estimate  the  amount  of  heat  advected  into  the 
Arctic  basin  by  the  WSC.  He  calculated  a  yearly  average  of  about  67  TW, 
a  summer  value  of  100  TW,  and  a  winter  value  of  50  TW  (one  tera-watt 
(TW)  *  10*2W).  Mandel  (1978)  calculated  a  mean  heat  transport  of  28  TW 
in  the  0-500  m  depth  interval  along  78°N  between  2°E  and  the  coast  of 
Spitsbergen  over  the  period  1949-74.  He  identified  a  seasonal  pattern 
with  a  maximum  in  late  July  and  a  minimum  in  April. 

Aagaard  and  Greisman  (1975)  reviewed  previous  budgets  and  formu¬ 
lated  new  mass  and  heat  budgets  for  the  Arctic  Ocean.  They  showed  that 
the  WSC  plays  a  major  part  in  the  heat  budget,  transporting  about  two- 
thirds  of  the  oceanic  heat  flux  into  the  Arctic.  It  is  likely  then, 
that  large  variations  in  this  flow  at  very  long  time  scales  (months  to 
years)  as  suggested  by  earlier  investigators  may  be  reflected  in 
changes  in  the  oceanic  regime  in  the  Arctic  Basin. 

1.2.3.  Temporal  variability 

The  previous  discussions  of  the  WSC  have  been  largely  in  terms  of 
mean  values  of  mass  and  heat  transports,  but  there  is  ample  evidence  of 
significant  low-frequency  variability  about  the  mean  values.  Figures 
8,  9,  10  and  Table  1  all  demonstrate  the  probable  existence  of  such 
variability. 

Figure  8,  taken  from  Lee  (1961b),  shows  2-3  month  time  scale  fluc¬ 
tuations  of  volume  transport  as  derived  from  hydrographic  data  across 
the  West  Bear  Island  section  over  the  interval  1949-1959.  Figure  9, 
from  Mandel  (1979)  shows  the  year-to-year  fluctuations  of  volume  and 
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Figure  8.  Volume  transport  estimates  (Sv)  for  the  West  Spitsbergen 
Current  across  the  Bear  Island  section  (74°30'N)  in  the  0- 
400  m  layer.  The  curve  represents  means  over  all  years^ 
Positive  sign  Indicates  northward  transport.  One  Sv  *  10° 
nrs"1.  (from  Lee,  1961b). 
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Volume  (Sv,  solid  line  1)  and  heat  (TW(  dashed  line  2) 
transports  of  the  West  Spitsbergen  Current  northward  across 
the  Barentsberg  section  at  78°N.  One  Sv  *  10b  nrs“A;  one 
TW  -  10I?W.  (from  Mandel,  1979). 


Northward  volume  transport  (Sv)  of  the  West  Spitsbergen 
Current  in  1971-72  computed  from  direct  current  measure¬ 
ments  In  conjunction  with  velocity  contours  determined 
from  hydrographic  data.  One  Sv-  10°  nrs  .  (from  Grels- 
raan,  1976) 
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Table  1.  Annual  mean  volume  and  heat  transports,  and  transport 
temperature  of  Atlantic  water  entering  the  Arctic  basin 
•  across  the  Barentsburg-ice  edge  section  west  of  Spits- 

bergen  along  78°N,  (from  Timofeyev,  1962). 

One  Sv  *  10B  nrs-1;  one  TW  *  10rzW. 
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heat  transports  northward  across  the  Barentsburg  section  at  78°N  from 
1949-77.  These  results  were  obtained  from  a  linear  relationship  formu¬ 
lated  between  the  monthly-averaged  temperature  and  water  level  at 
Barentsburg  and  the  temperature  and  geostrophic  transport  of  the  USC. 
Table  1  lists  the  volume  and  heat  transports  for  various  years  between 
1933  and  1960  determined  by  Timofeyev  (1962)  from  hydrographic  data 
obtained  along  78°H  from  Spitsbergen  westward  to  the  ice  edge.  The 
largest  year-to-year  variations  are  on  the  order  of  50%.  The  monthly 
transport  estimates  in  Fig.  10,  from  Greisman  (1976),  were  derived  from 
direct  current  measurements,  and  exhibit  monthly  to  seasonal  fluctua¬ 
tions.  In  contrast  to  the  long-term  variability  noted  -above.  Hill  and 
Lee  (1957)  provided  evidence  of  variability  with  a  time  scale  on  the 
order  of  a  few  days.  They  reported  on  four  workings  of  a  section  west 
of  Bear  Island  within  48  hours  which  gave  transport  values  from  about 
0.1  Sv  southward  to  1  Sv  northward.  Neglecting  detailed  comparisons, 
it  is  clear  from  the  above  that  the  volume  flow  may  be  variable  over 
time  scales  from  days  to  years. 

Using  hydrographic  data  taken  in  late  summer  of  successive  years 
during  1965-1969,  Dickson  and  Doddington  (1968,  1970)  detected  a  trend 
of  Increasing  temperature  and  salinity  in  the  flow  west  of  Spitsbergen 
at  79°N.  However,  they  also  pointed  out  the  presence  of  short-term 
variability.  For  example.  In  1969  the  high-salinity  core  (>35.1)  in 
the  section  at  79°40'N  was  larger  than  that  at  79°N,  and  the  tempera¬ 
tures  associated  with  the  core  were  higher  than  at  79°N.  They  there¬ 
fore  questioned  the  validity  of  regarding  the  hydrographic  changes 
between  successive  annual  surveys  as  true  interannual  changes,  because 
of  aliasing  by  the  short-term  variability  of  the  flow. 

The  temperature  of  the  WSC  also  appears  to  fluctuate  in  time, 
which  In  turn  Influences  heat  transport  estimates.  Solonltsina  (1976, 
1978)  has  described  low-frequency  temperature  fluctuations  In  the  WSC, 
both  seasonal  and  Interannual.  Table  1  Includes  temperatures  used  by 
Timofeyev  (1962)  to  calculate  heat  transports  of  the  WSC  at  78°N  for 
various  years  between  1933  and  1960.  These  temperature  values  repre¬ 
sent  essentially  the  depth-averaged  temperatures  over  the  upper  1000  m. 
The  greatest  apparent  interannual  change,  0.80C,  occurred  between  1933 
and  1934,  but  other  years  also  showed  large  variability. 
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1.2.4  Spatial  variability 

Many  investigators  have  presented  evidence  of  spatial  variability 
in  the  flow  on  the  eastern  side  of  the  Greenland  and  Norwegian  seas. 
Helland-Hansen  and  Nansen  (1909)  recognized  a  pattern  of  waviness  in 
the  isollnes  of  their  vertical  sections  in  the  Norwegian  Atlantic 
Current.  They  termed  the  undulations  "puzzling  waves",  and  postulated 
that  they  were  due  either  to  large  lateral  variations  of  the  currents 
or  to  Internal  waves,  favoring  the  former.  Later  observations 
(Helland-Hansen,  1934}  seemed  to  support  the  idea  that  the  waviness 
corresponded  with  lateral  variations  in  the  flow.  Lee  (1959)  concluded 
that  the  wav i ness  was  probably  due  to  vortices  about  10  km  in  diameter 
and  that  the  system  of  vortices  changed  only  slowly.  Saelen  (1963) 
believed  that  the  waviness  was  due  to  eddies,  both  stationary  and 
moving,  and  Dickson  (1972)  suggested  that  the  Norwegian  Atlantic  Cur¬ 
rent  was  made  up  of  a  number  of  streams  fairly  constant  in  position, 
producing  a  pattern  of  wavy  isolines  with  a  scale  of  40-80  km  which 
approximately  followed  bottom  contours. 

Earlier,  Alekseev  and  Istoshin  (1956)  had  published  a  chart  of 
surface  currents  using  primarily  data  obtained  from  post-World  War  II 
Soviet  hydrographic  work.  Although  they  gave  no  estimates  of  the  lat¬ 
eral  scale,  the  lower  limit  would  correspond  to  their  station  spacing, 
which  was  typically  20-40  km.  Their  chart  portrayed  the  Norwegian 
Atlantic  Current  as  splitting  into  eastern,  western  and  central 
branches  which  reunite  farther  north  to  form  the  WSC.  Furthermore, 
they  asserted  the  existence  of  elongated  eddies  between  the  various 
streams  and  branches,  and  they  concluded  that  a  close  link  existed 
between  the  patterns  of  flow  and  the  bottom  topography.  Kislyakov 
(1960),  looking  at  hydrographic  observations  at  intervals  of  24  km 
between  7°-16°E  across  74°30'N,  identified  four  different  northward 
flowing  streams.  According  to  his  findings,  the  number  of  streams 
remained  constant,  but  their  positions  varied  laterally  16  to  24  km  to 
either  side  of  a  mean  position.  Interestingly,  he  also  found  that 
velocity  fluctuations  among  the  streams  tended  to  be  out  of  phase. 
Further  evidence  of  the  banded  structure  comes  from  a  chart  of  geopo- 
tential  topography  by  Eggvln  (1961)  which  clearly  shows  three  streams 
of  the  WSC  with  a  typical  scale  of  50  km  west  of  Bear  Island.  Rudels 
and  Andersson  (1982)  presented  a  velocity  section  across  Fram  Strait  at 
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79°N  which  was  derived  from  hydrographic  data  and  referenced  by  direct 
current  measurements.  A  region  of  reversed  velocity  in  the  WSC  sug¬ 
gested  the  scale  of  the  lateral  structure  in  the  flow  to  be  about  15 
km.  In  addition,  a  separate  data  set  containing  data  obtained  at 
intervals  of  about  10  km  shows  core-like  features  in  temperature  con¬ 
tours  in  a  section  along  80°N.  These  features  indicate  lateral  struc¬ 
ture  in  the  flow,  as  does  the  20/400  db  dynamic  height  anomaly  which 
exhibits  three  obvious  maxima  between  2°-ll°E  (Fig.  11). 

With  respect  to  the  WSC  and  the  Norwegian  Atlantic  Current,  the 
implication  of  these  investigations  is  that  the  entire  system  is  com¬ 
posed  of  a  filamentous  structure  of  northward  flowing  streams  separated 
by  elongated  eddies.  The  length  scales  associated  with  the  lateral 
structure  appear  to  range  from  kilometers  to  tens  of  kilometers. 
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Figure  11.  Temperature  contours  and  dynamic  height  anomaly  (dyn.  m, 
20/400  db)  west  of  Spitsbergen  along  80°N  obtained  Novem¬ 
ber  26-27,  1977.  Polarsirkel  cruise,  (prepared  by  S. 
Hillman). 


2.  Experimental  design,  data  sources 


2.1  Mooring  details 

During  successive  years  beginning  in  July  1976,  direct  current 
measurements  were  made  at  various  locations  in  the  WSC.  Figure  3  shows 
the  positions  of  the  four  moorings  deployed  the  first  year.  The  fol¬ 
lowing  year,  after  recovery  of  the  previous  year's  moorings,  new  moor¬ 
ings  were  deployed  at  locations  B  and  C  and,  after  that,  only  at 
location  C.  The  taut-line  moorings  (using  Kevlar  line)  carried  two 
current  meters  each,  suspending  the  upper  instrument  at  a  depth  of 
about  100  m  and  the  lower  at  about  400  m.  Sampling  intervals  were  40 
minutes  for  the  first  year  and  60  minutes  for  the  following  years. 
Table  2  gives  the  basic  deployment  details  for  all  the  current  meters 
used  in  this  study. 

In  order  to  minimize  the  effect  of  the  recirculating  flow  on 
estimates  of  volume  transport,  i.e.,  to  restrict  the  estimates  to  that 
which  actually  passes  into  the  Arctic  basin,  the  mooring  section  was 
located  along  79°N,  which  is  north  of  the  region  where  the  recircula¬ 
ting  flow  curves  westward.  The  varying  lateral  spacing  of  the  moorings 
was  designed  to  allow  the  estimation  of  several  different  length  scales 
across  the  flow,  while  the  300  m  vertical  separation  of  current  meters 
permitted  some  investigation  of  the  vertical  structure  in  the  upper 
part  of  the  water  column. 

Each  mooring  used  an  AMF  model  242  release,  Aanderaa  RCM-4  current 
meters,  and  subsurface  steel  sphere  flotation.  The  internally  record¬ 
ing  current  meters  used  a  modified  Savonius  rotor  magnetically  coupled 
through  the  pressure  case  to  count  revolutions  over  the  sampling  inter¬ 
val,  providing  an  average  speed  for  the  interval.  A  thermistor  and 
magnetic  compass  gave  Instantaneous  temperature  and  direction  values  at 
the  end  of  each  sampling  Interval.  A  large  vane  oriented  the  current 
meter  with  the  flow. 

The  current  meter  rotors  were  calibrated  by  towing  at  known  speeds 
in  a  tow  tank  operated  by  the  U.S.  Army  Corps  of  Engineers  at  Bonne¬ 
ville,  Washington.  Thermistors  were  calibrated  at  the  NOIC  facility  in 
Bellevue,  Washington,  and  magnetic  compasses  were  calibrated  at  Sand 
Point  Naval  Air  Station  In  Seattle,  Washington.  The  calibration  proce¬ 
dures  were  the  same  as  detailed  by  Grelsman  (1976). 


The  Aanderaa  data  tapes  were  read,  edited  and  converted  to  7-track 
computer-compatible  tape  at  the  Department  of  Oceanography  Computer 
Services  facility.  An  initial  machine  scan  produced  a  listing  of  the 
data  and  flagged  the  speed  values  considered  unreasonable.  The  data 
listing  was  used  to  establish  start  and  end  times  for  the  record  man¬ 
ually  and  to  correct  obviously  erroneous  speed  values  by  linearly 
interpolating  between  the  nearest  acceptable  data  points.  Subsequent 
analysis  of  the  edited  data  was  done  on  the  computers  at  the  University 
of  Washington  Academic  Computer  Center. 

2.2  Measurement  errors 

Aanderaa  RCM-4  current  meters  have  been  used  widely,  and  their 
performance  is  well  documented.  After  calibration,  errors  in  speed 
indications  may  range  up  to  2-3  cm  s"1,  and  compass  errors  are  probably 
in  the  vicinity  of  5°.  Thermistors  are  accurate  to  within  +  0.1C  and 
appear  stable  over  time.  Among  the  current  meters,  the  maximum  speeds 
recorded  ranged  to  75  cm  s“*,  while  the  means  were  about  5-15  cm  s_1, 
so  that  errors  were  approximately  an  order  of  magnitude  greater  than 
the  expected  instrument  error.  With  the  upper  current  meters  at  100  m 
depth  and  with  subsurface  flotation,  contamination  due  to  surface 
effects  or  to  pumping  by  vertical  mooring  motions  is  within  the  error 
limits  of  the  data. 

2.3  Hydrographic  data  sources 

Table  3  lists  ship  cruises  which  provided  hydrographic  data  used 
in  this  investigation.  Hydrographic  coverage  in  the  mooring  area  is 
sparse,  especially  in  winter.  Some  stations  were  taken  in  conjunction 
with  mooring  deployment  and  recovery,  providing  hydrographic  data  con¬ 
current  with  direct  measurements.  In  general,  suitable  geographic 
position  and  confidence  in  the  accuracy  of  the  data  were  the  two  pri¬ 
mary  considerations  in  selecting  the  hydrographic  stations,  although 
the  most  serious  problem  was  probably  wide  lateral  spacing. 

2.4  Designation  convention  for  current  records 

A  simple  convention  is  adopted  for  Identifying  the  year,  mooring 
location,  and  depth  of  the  various  current  records  mentioned  in  the 
subsequent  discussions.  The  format  used  is  My-ddd;  where  M  is  the 


Table  2.  Details  of  current  meter  moorings  deployed  near  79°N  in  the 
West  Spitsbergen  Current. 


Current 

Record 

Water 

Depth 

(m) 

Latitude 

(°N) 

Longitude 

(°E) 

Depth 

(m) 

Start 

Stop 

Al-98 

1248 

78°58.8‘ 

006°54.1' 

98 

07/16/76 

09/02/77 

Al-398 

398 

07/16/76 

09/04/77 

Bl-398 

1198 

78°57.0' 

007°38.0' 

398 

07/15/76 

09/04/77 

Cl-92 

1142 

78°57.8‘ 

007°53.6‘ 

92 

07/15/76 

09/01/77 

Cl-392 

392 

07/15/76 

08/31/77 

Dl-93 

343 

78°59.2' 

008° 32.1* 

93 

07/15/76 

08/17/77 

B2-110 

1210 

78°57.2' 

007°35.6' 

110 

09/05/77 

08/19/78 

B2-410 

410 

09/05/77 

06/20/78 

C2-100 

1150 

78°58.3‘ 

007° 50. 6' 

100 

09/05/77 

08/19/78 

C2-400 

400 

09/05/77 

08/19/78 

C3-100 

1150 

78°58.8' 

007°50.0' 

100 

08/20/78 

06/23/79 

C3-400 

400 

08/20/78 

06/07/79 

C5-84 

1124 

78°47.7' 

007°33.8' 

84 

09/09/80 

07/29/81 

C5-384 

384 

09/09/80 

04/20/81 

C5-984 

984 

09/09/80 

06/08/81 

mooring  location,  either  A,  B,  C  or  D;  y  is  the  year  number  with  y  =  1 
for  1976-77,  y  =  2  for  1977-78,  etc.;  and  ddd  is  the  depth  in  meters  at 
which  the  current  record  was  obtained.  Thus  Al-98  refers  to  the  cur¬ 
rent  record  obtained  at  mooring  location  A  during  1976-77  at  a  depth  of 
98  m,  and  C3-400  refers  to  mooring  location  C,  1978-79,  and  depth  of 
400  m.  The  designation  CB  refers  to  bridged  current  records  from 
location  C. 


Table  3.  Cruises  which  provided  hydrographic  data  used  in  this  study 


Ship 

Year 

Month 

NOOC  ID 

Armauer  Hansen 

1926 

8 

580259 

Atka 

1954 

1 

310445 

G.O.  Sars 

1956 

10 

580999 

Ernest  Holt 

1957 

5,  7 

740709 

G.O.  Sars 

1958 

4,  5,  6 

580013 

Johan  Hjort 

1958 

9 

580957 

Atka 

1962 

oo 

310965 

Edisto 

1964 

9 

310296 

G.O.  Sars 

1965 

5 

580064 

Arvakur 

1971 

9 

Meteor 

1972 

7 

060157 

G.O.  Sars 

1976 

7 

Johan  Hjort 

1977 

9 

Polarsirkel 

1977 

11 

.  Volume  transport  of  the  West  Spitsbergen  Current 


3.1  Data  preparation  and  transport  calculation  method 

The  current  records  from  moorings  A,  B,  C  and  D  were  used  to 
calculate  the  volume  transport  of  the  WSC  during  the  1976-77  deployment 
interval.  Tides  and  higher  frequency  signals  were  first  removed  by 
applying  a  filter  with  a  cutoff  of  40  hours.  Filtering  was  accom¬ 
plished  by  Fourier  transforming,  zeroing  the  coefficients  which  corres¬ 
ponded  to  those  frequencies  beyond  the  desired  cutoff,  cosine-tapering 
the  truncated  end,  and  then  transforming  back  to  the  time  domain  to 
construct  the  filtered  data  series.  The  resulting  time  series  were 
then  sub-sampled  at  24-hour  intervals.  To  arrive  at  a  suitable  averag¬ 
ing  interval  for  the  data,  the  autocorrelation  function  was  calculated 
for  each  of  the  north-south  velocity  records.  After  Davis  (1976)  and 
Bryden  and  Pillsbury  (1977),  the  time  scale  of  independence  (TSI)  for 
the  respective  records  is: 

x=+mAt 

TSI  «  E  C2(t)  At  , 

T=-mAt 

where  mU  is  the  greatest  lag  t  for  which  the  autocorrelation  C(t  ) 
differs  from  zero  at  the  95X  confidence  level.  The  TSI  is  the  time 
interval  beyond  which  successive  velocity  measurements  may  be  con¬ 
sidered  to  be  statistically  independent  of  each  other.  The  TSI  for  the 
downstream  velocity  components  varied  from  about  8  days  to  19  days  and 
averaged  13.7  days.  The  final  data  series  used  to  estimate  volume 
transports  were  then  obtained  by  averaging  the  24-hour  data  over  14-day 
intervals. 

Figure  12  shows  the  cross-section  used  for  calculating  the  trans¬ 
port  of  the  flow.  The  spatial  limits  for  the  flow  were  chosen  according 
to  criteria  outlined  below,  and  within  the  zero  isotach  the  current 
meters  were  assigned  sub-areas  to  which  their  respective  measurements 
applied.  The  deeper  current  records  at  moorings  B  and  C  were  highly 
correlated  (see  Table  4),  and  so  those  velocities  were  averaged  before 
making  the  transport  calculations. 

The  western  limit  of  the  current  was  taken  to  be  3°30'E  after 
considering  hydrographic  data  and  the  direct  current  measurements  from 
1971-72.  The  hydrographic  data  showed  a  sharp  decrease  of  the  tempera- 


Cross  section  used  in  this  study  for  calculating  volume 
and  heat  transports  of  the  West  Spitsbergen  Current.  Seg¬ 
ment  designations  correspond  with  mooring  location  desig¬ 
nations.  Flow  was  considered  barotropic  below  600  m. 
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Table  4.  Cross-correlation  coefficients  for  records  of  alongshore 
velocities  at  four  mooring  sites  in  the  West  Spitsbergen 
Current  during  1976-77.  Ru  and  R^  respectively  represent 
the  cross -correlation  at  zero  lag  between  upper  and  lower 
pairs  of  current  meters  on  the  various  moorings. 


R1  *  -.41 


Ru  *  -.42 
Ri  *  -.52 


Ru  =  +.16 


Ri  *  +.71 


*RU  *  -.03** 


♦Maximum  value  was  +.30  with  D  leading  C  by  37  days. 

♦♦This  value  not  significantly  different  from  zero 
at  the  95%  confidence  level. 
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ture  at  30  m  depth  between  3°-4°E.  This  temperature  transition  is  one 
indication  of  the  frontal  boundary  between  the  WSC  and  the  East  Green¬ 
land  Current,  and  the  location  is  approximately  the  same  as  that  deter¬ 
mined  by  Greisman  (1976)  by  linear  interpolation  between  sets  of 
current  meter  data  taken  both  near  the  southward  flowing  East  Greenland 
Current  and  near  the  highest  velocity  core  of  the  WSC. 

The  eastern  limit  was  set  at  9°30‘E,  well  shoreward  of  the  shelf 
break.  Transports  east  of  9°30’E  calculated  from  geostrophic  veloci¬ 
ties  referenced  to  250  db  were  about  0.3  Sv,  so  that  excluding  this 
region  probably  introduces  an  error  of  less  than  about  10-15%,  even  if 
a  barotropic  component  of  equal  magnitude  is  included. 

Figure  13  shows  contours  of  dynamic  height  referenced  to  600  db. 
The  dynamic  height  values  were  obtained  by  averaging  the  available  data 
from  various  cruises  near  79°N  by  one-degree  bins  between  2°-4°E  and  by 
half-degree  bins  between  4°-9°E.  The  number  of  stations  averaged 
ranged  from  2  to  11,  with  the  densest  coverage  between  6°-9°E.  Mooring 
locations  A  and  C  appear  to  lie  in  separate  velocity  cores  (Fig.  13), 
and  the  dynamic  field  indicates  the  presence  of  a  third  velocity  core 
between  5°-6°E.  With  the  contour  plot  of  dynamic  heights  as  a  guide, 
the  transport  section  was  partitioned  into  five  segments.  Three  of 
these  corresponded  to  the  three  mooring  locations  A,  C  and  D,  the 
fourth  encompassed  the  third  core  between  5°-6°E,  and  the  fifth 
included  the  westernmost  fringe  across  which  the  velocity  was  assumed 
to  linearly  approach  zero  at  3°30'E.  The  vertical  boundary  between 
mooring  A  and  C  segments  was  assumed  to  be  the  midpoint  between  them, 
and  that  between  mooring  C  and  0  segments  was  assumed  near  the  shelf 
break.  Preliminary  transport  calculations  for  only  mooring  A  and  C 
segments  yielded  values  of  about  1  to  2  Sv  in  each.  The  transport  of 
the  indicated  core  between  5°-6°E  was  assumed  to  be  comparable.  Spe¬ 
cifically,  the  segment  for  that  core  was  assumed  to  be  the  same  size  as 
mooring  segment  A,  and  the  mooring  A  current  records  were  taken  to  be 
representative  of  the  velocities  in  the  assumed  segment.  Essentially, 
mooring  segment  A  was  doubled  in  area  to  include  the  westernmost  core 
at  5°-6°E  indicated  by  the  hydrographic  data. 

Baroclinlc  contributions  to  the  transport  were  considered  negligi¬ 
ble  below  600  m.  This  depth  was  determined  by  averaging  and  contouring 
the  available  hydrographic  data  by  one-half  degree  intervals  between 


Contours  of  dynamic  height  anomaly  (dyn.  m,  0/600  db)  west 
of  Spitsbergen  along  79aN.  Values  were  obtained  by  avera¬ 
ging  available  hydrographic  data  by  one-degree  bins 
between  2°-4°E  and  by  half -degree  bins  between  4^-9°E. 
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5°-9°E,  calculating  average  dynamic  depths,  and  then  selecting  the 
geometric  depth  at  which  the  vertical  gradient  of  dynamic  depth  became 
nearly  uniform  across  the  flow.  The  600  db  reference  level  corresponds 
with  that  determined  by  Palfrey  (1967)  and  Greisman  (1976).  Calcula¬ 
tions  using  more  recent  hydrographic  data  from  various  stations  across 
the  current  indicated  that  the  contribution  of  geostrophic  transport 
from  the  600-1000  m  depth  interval  could  amount  to  10-25*  of  that  above 
600  m;  this  was  judged  sufficiently  small  that  baroclinicity  below  600 
m  was  neglected.  The  velocity  between  600  m  and  the  lower  boundary  of 
the  section  was  assumed  to  be  entirely  barotropic  and  was  assigned  a 
constant  fraction  of  the  velocity  measured  at  400  m  depth.  This  frac¬ 
tion,  0.86,  was  determined  by  assuming  a  linear  velocity  profile, 
performing  a  linear  regression  between  the  upper  and  lower  northward 
current  speeds  at  moorings  A  and  C,  and  using  the  resulting  slope  to 
extrapolate  the  velocity  to  600  m. 

From  9°30'E,  the  zero  isotach  coincides  with  the  sea  bottom  west¬ 
ward  to  about  7°E  where  the  bottom  slope  increases  by  a  factor  of 
three.  From  about  7°E  the  zero  isotach  extends  horizontally  at  a  depth 
of  about  1250  m  westward  to  5°20'E,  which  marks  the  eastern  edge  of  the 
assumed  transition  zone.  From  this  point  the  zero  isotach  angles 
upward  to  the  west  and  Intersects  the  surface  at  3°30'E.  This  shape 
resembles  that  used  by  Greisman  (1976,  p.  38),  although  it  was  arrived 
at  differently.  It  also  incorporates  an  assumption  that  speeds  in  the 
WSC  diminish  across  the  transition  region  toward  the  East  Greenland 
Current,  essentially  providing  a  conservative  transport  estimate. 
Placing  the  zero  Isotach  at  1250  m  depth  admittedly  contradicts  the 
assumption  that  below  600  m  the  flow  Is  entirely  barotropic.  Including 
the  area  in  the  section  between  1250  m  depth  and  the  bottom  between 
3°30*  and  7°E  and  basing  the  velocity  on  the  mooring  A  current  measure¬ 
ments  would  at  most  add  slightly  more  than  1  Sv  to  the  total  transport. 
Direct  current  measurements  do  show  a  mean  northward  velocity  component 
of  3.9  cm  s"1  during  1980-81  at  984  m  depth  at  mooring  location  C,  but 
there  are  no  direct  current  data  available  for  the  deep  region  west  of 
about  7°E.  Omitting  this  western  deep  region  by  placing  the  zero 
Isotach  at  1250  m  thus  reflects  the  choice  to  seek  a  conservative 
estimate  of  the  total  transport. 
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3.2  Lateral  structure  of  the  flow 

Rather  than  simply  being  spatially  uniform,  or  even  uniformly 
varying,  the  north-south  (v-component)  velocity  field  in  the  mooring 
section  west  of  Spitsbergen  shows  significant  variation  across  the 
stream.  Note  the  waviness  in  the  geostrophic  field  (Fig.  13),  with  its 
implied  geostrophic  velocity  structure  between  5°-9°E.  Since  one 
degree  of  longitude  at  79°N  corresponds  to  about  21  km,  the  dynamic 
height  contours  indicate  a  fairly  small  spatial  scale  of  about  20  km 
for  the  variability.  Figure  11  shows  similar  evidence  of  lateral  struc¬ 
ture  at  80°N. 

The  weekly  mean  velocities  for  1976-77  plotted  as  stick  diagrams 
(Fig.  14)  illustrate  the  spatial  and  temporal  variability.  One  imme¬ 
diate  observation  is  that  the  speeds  varied  markedly  among  the  moor¬ 
ings,  more  so  laterally  than  vertically.  Sudden  reversals  occurred  at 
moorings  B  and  C  around  January  1,  1977  which  are  not  common  to  the 
records  from  A  and  0.  Over  the  last  three  months  of  the  deployment, 
speeds  at  mooring  C  increased  by  a  factor  of  three.  No  such  strong 
trend  appeared  at  any  of  the  other  moorings,  even  B  which  was  only  5.5 
km  away.  In  addition,  the  speeds  at  Cl-392  actually  exceeded  those  at 
Cl-92  at  C  over  those  final  3  months.  As  well,  the  directional  varia¬ 
bility  differed  markedly  among  the  moorings;  the  current  was  most 
nearly  unidirectional  at  mooring  D  and  least  so  at  mooring  A.  This  is 
likely  due  to  several  factors:  1)  the  strong  influence  of  topography 
at  mooring  D  constrained  the  flow  there  to  follow  bottom  contours;  2) 
the  greater  distance  of  mooring  A  than  mooring  C  from  the  very  steep 
portion  of  the  slope  region,  which  allowed  for  less  of  a  channeling 
influence  at  mooring  A;  and  3)  the  location  of  mooring  A  near  a  sharp 
discontinuity  in  bottom  slope  made  for  increased  sensitivity  to  small 
perturbations  which  increased  the  incidence  of  fluctuations  in  direc¬ 
tion  as  the  flow  adjusted  to  conserve  potential  vorticity. 

The  results  of  cross-correlation  calculations  (Table  4)  show  the 
cross-stream  structure  to  be  rather  complicated.  Records  of  C1-92V  and 
D1-93V  are  essentially  uncorrelated  at  zero  lag,  so  that  the  14  km 
separation  between  the  two  is  greater  than  the  local  correlation  length 
scale  (refer  to  Fig.  3  for  mooring  locations).  In  contrast,  the  analo¬ 
gous  correlation  between  A1-98V— D1-93V,  separated  by  35  km,  is  0.16, 
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Figure  14.  Stick  diagrams  of  weekly  mean  velocity  for  all  current 
records  obtained  In  1976-77.  Current  meter  locations  are 
shown  at  the  right. 
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significantly  different  from  zero  at  the  95%  confidence  level.  Inter¬ 
estingly,  the  cross-correlations  between  A1-398V— B1-398V,  and  between 
A1-98V— C1-92V  and  A1-398V— C1-392V  are  consistently  and  significantly 
negative.  This  means  that  periods  when  the  velocity  at  mooring  C 
exceeded  the  local  northward  mean  at  C  corresponded  with  periods  when 
the  velocity  at  A  was  less  than  the  northward  mean  at  A,  and  vice 
versa.  One  effect  of  this  complicated  cross-stream  structure  is  to 
increase  the  error  in  estimates  of  the  volume  transport. 

3.3  Estimation  of  rms  error 

Fandry  and  Pillsbury  (1979)  presented  an  objective  method  of 
estimating  the  barotropic  transport  across  an  oceanographic  section, 
and  they  applied  it  to  the  transport  of  the  Antarctic  Circumpolar 
Current  through  Drake  Passage.  With  this  method,  the  values  for  the 
rms  error  in  the  objective  estimate  of  the  transport  depended  upon  the 
number  of  current  meters  at  a  given  level,  on  the  correlation  length 
scale  of  the  flow,  and  on  the  noise.  Specifically,  using  a  spatial 
correlation  function  everywhere  positive,  they  found  that  for  current 
meter  spacing  less  than  or  equal  to  the  correlation  length  and  for 
measurement  noise  variance  less  than  10%,  the  rms  error  in  the  trans¬ 
port  estimate  is  less  that  20%.  For  current  meter  separations  larger 
than  the  correlation  length  or  for  a  spatial  correlation  function  which 
crosses  zero,  the  error  can  be  2  to  3  times  greater. 

Applying  this  method  rigorously  to  the  WSC  is  not  possible  for 
several  reasons:  1)  the  spatial  correlation  function  is  complicated 
and  does  cross  zero,  2)  the  current  meter  moorings  are  not  equidis- 
tantly  spaced,  so  that  the  velocities  cannot  simply  be  linearly  aver¬ 
aged,  and  3)  the  current  meter  mooring  separation  is  not  always  less 
than  the  correlation  length  scale.  Nevertheless,  the  Fandry  and  Pills¬ 
bury  method  provides  at  least  a  guide  line  for  estimating  the  probable 
error  in  the  transport.  Their  Fig.  2a  (p.  451)  shows  that  for  a 
spatial  correlation  function  which  crosses  zero  and  for  data  spacing 
somewhat  greater  than  the  correlation  length  scale,  the  rms  error  in 
calculated  transports  is  in  the  neighborhood  of  0.8,  i.e.,  about  80%  of 
the  transport  value.  Thus,  the  probable  errors  in  the  transport  esti¬ 
mates  of  the  WSC  are  roughly  as  large  as  the  estimates  themselves. 
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3.4  Volume  transport  based  on  alj.  current  records.  1976-77 

The  solid  curve  of  Fig.  15  is  a  plot  of  the  volume  transport  of 
the  WSC  computed  from  direct  current  records  obtained  during  1976-77. 
Calculations  were  based  on  the  cross  section  portrayed  in  Fig.  12. 


3.4.1  Mean  and  range 

During  1976-77,  the  transports  averaged  by  two-week  periods  over 
the  deployment  interval  ranged  from  1.4  Sv  southward  to  11.9  Sv  north¬ 
ward.  The  average  was  5.6  Sv  northward  and  of  this,  approximately  1.9 
Sv  was  below  600  m  depth.  Assuming  that  the  barotropic  flow  per  unit 
area  of  transport  cross  section  above  and  below  600  m  was  about  the 
same,  then  the  ratio  of  cross-sectional  areas  above  and  below  600  m 
gives  a  barotropic  component  above  600  m  of  2.7  Sv,  and  together  with 
the  contribution  from  deeper  than  600  m,  the  total  barotropic  transport 
for  the  entire  section  is  4.6  Sv.  The  difference  between  the  total 
transport  and  the  barotropic  transport  is  then  1.0  Sv  and  corresponds 
to  the  baroclinic  contribution.  This  is  only  one-half  the  baroclinic 
transport  value  of  2  Sv  between  3°30'-8°30,E  referenced  to  600  db 
calculated  from  hydrographic  data  obtained  during  a  single  cruise  in 
November,  1977.  However,  it  does  agree  more  closely  with  a  value  of 
0.7  Sv  calculated  from  the  bin-averaged  hydrographic  data. 

Timofeyev  (1963)  calculated  an  absolute  value  of  about  3.0  Sv  for 
the  volume  transport  of  the  WSC.  This  figure  is  less  than  the  total 
transport  of  5.6  Sv  estimated  from  the  present  current  meter  data,  but 
roughly  three  times  the  estimated  value  of  1.0  Sv  for  the  baroclin’c 
contribution  during  1976-77.  Timofeyev's  hydrographic  data  were 
sparsely  distributed  among  months  and  years,  so  that  in  9  of  the  13 
years,  transport  estimates  were  drawn  from  only  a  single  section.  Man- 
del  (1976)  estimated  a  mean  volume  transport  of  the  WSC  in  the  0-1000 
m  layer  to  be  5.7  Sv.  His  computational  method  involved  formulating  a 
relationship  between  the  water  level  at  Barentsburg  and  the  geostrophic 
velocity  across  78°N  as  determined  from  hydrographic  data  and  refer¬ 
enced  to  1000  meters.  He  then  used  monthly  mean  water  level  data  to 
estimate  the  inflow  of  water  into  the  Arctic  basin  over  the  interval 
1949-73.  Greisman  (1976)  calculated  a  total  mean  annual  value  of  7  Sv. 
This  is  the  only  other  estimate  based  on  direct,  long-term  current 

measurements. 


Figure  15.  Estimates  of  volume  (Sv)  and  heat  (TW)  transports  of  the 
West  Spitsbergen  Current  in  1976-77.  Estimates  are  based 
on  two-week  averages  of  direct  measurements  of  alongshore 
velocity  and  temperature  obtained  west  of  Spitsbergen 
along  79°N*  Pjositive  value  means  northward  transport. 
One  Sv  *  106  nrs  ,  one  TW  *  lO^W. 
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Various  other  investigators:  Lee  (1961b),  Kislyakov  (1960),  Mosby 
(1938),  and  cited  in  Mosby  (1962);  Laktionov  (1959),  Zaitsev  (1960), 
Zaitsev,  et  al.  (1961),  Treshnikov  (1960),  have  estimated  the  mean 
volume  transport  of  the  WSC  using  hydrographic  data.  Their  estimates 
ranged  from  about  1  to  more  than  4  Sv,  but  they  were  based  on  data  from 
various  latitudes  and  seasons,  on  sections  of  different  lateral  extent 
across  the  flow,  and  on  differing  reference  levels,  so  that  a  compari¬ 
son  does  not  seem  possible. 

After  investigating  the  reference  levels  used  in  the  WSC  by  pre¬ 
vious  investigators.  Palfrey  (1967)  reached  the  conclusion  that  either 
600  db  or  1000  db  could  be  used  (Timofeyev  chose  the  1000  db  reference 
level).  However,  as  shown  in  Table  5,  the  mean  velocity  shears  in  the 
upper  part  of  the  water  column  calculated  from  direct  current  measure¬ 
ments  show  that  the  degree  of  baroclinicity  does  vary  markedly.  Over 
the  same  depth  range,  the  vertical  shear  ranged  from  0.003  to  0.0163  cm 
s"1  per  m  between  1976  and  1981.  Specifically,  the  vertical  shear 
during  1976-77  was  about  one-half  that  observed  in  any  other  year,  and 
the  values  of  0.84  and  0.86  for  the  ratios  of  mean  velocities  (Al- 
98V/A1-398V  and  C1-92V/C1-392V,  respectively)  were  noticeably  higher 
than  for  other  years.  This  is  consistent  with  only  a  small  baroclinic 
contribution  to  the  transport  during  1976-77. 

However,  note  that  in  1980-81,  the  only  period  which  included  a 
deep  moored  current  meter  near  1000  m,  the  mean  vertical  velocity  shear 
between  84  m  and  384  m  was  greatest,  but  there  was  also  sianificant 
shear  between  384  m  and  984  m  (0.007  cm  s_1  per  m)  as  well  as  a  mean 
velocity  of  3.9  cm  s-1  at  984  m.  In  that  year,  neither  the  assumption 
of  no  baroclinicity  below  600  m  nor  the  assumption  of  a  zero-velocity 
reference  level  at  1000  db  would  have  been  correct.  The  direct  current 
evidence  shows  that  major  differences  between  various  estimates  of  the 
baroclinic  transport  because  of  temporal  variability  in  the  baroclinic 
contribution  are  to  be  expected. 

A  comparison  of  other  measured  velocities  also  supports  the  possi¬ 
bility  that  the  natural  variability  of  the  WSC  admits  large  dispari¬ 
ties  in  mean  transport  estimates.  At  nearly  the  same  mooring  location 
as  C  (79°N,  ~8°E),  the  average  velocities  at  the  upper  and  lower  cur¬ 
rent  meters  in  autumn  1971-72  shown  by  Greisman  (1976)  were  roughly 
four  times  those  during  autumn  1976-77.  The  year-long  mean  velocity  at 
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Table  5.  Mean  vertical  velocity  shear,  ratios  of  lower  (V^  to  upper 
(V,,)  velocities,  and  ratios  of  velocity  differences  to  upper 
velocities  for  moorings  with  vertically  separated  current 
meters.  Positive  shear  indicates  a  higher  northward  V..  than 
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the  upper  current  meter  was  about  twice  as  great  in  1971-72  than  in 
1976-77,  so  that  Greisman's  larger  transport  in  1971-72,  about  7  Sv,  is 
reasonable  for  that  time  (Forcing  by  wind  stress  curl  is  discussed  in 
Section  5,  and  one  result  is  that  the  magnitude  of  the  wind  stress  curl 
in  1971-72  was  nearly  twice  that  in  1976-79,  so  that  there  is  a  dynamic 
basis  for  the  difference  in  mean  velocity  during  these  respective 
periods).  On  the  other  hand,  only  one  of  the  1971-72  moorings,  with 
two  current  meters,  was  near  the  core  of  the  WSC,  and  thus  the  data 
coverage  was  much  sparser  than  in  1976-77.  Furthermore,  the  deeper  of 
the  current  meters  ceased  operating  about  75  days  after  it  was  deployed 
in  early  September,  so  that  Greisman  essentially  relied  upon  a  single 
year-long  current  record  to  estimate  the  annual  transport.  This  is 
important  because  of  another  point  of  comparison:  for  a  75-day  interval 
beginning  September  9  during  each  of  1971  and  1976,  the  ratios  of  mean 
upper  to  lower  northward  velocities  were  0.60  and  0.66,  respectively 
(Table  6).  The  analogous  ratio  of  year-long  means  during  1976-77  was 
0.86.  At  least  during  1976-77,  therefore,  the  vertical  shear  in  fall 
was  different  from  that  for  the  year  as  a  whole.  (Note,  however,  that 
the  two  sets  of  ratios  are  more  nearly  equal  in  1977  and  1978.)  This 
raises  the  question  of  whether  or  not  the  transport  based  on  '.he  fall 
currentpattern  in  1971-72  was  truly  representative  of  the  entire  year. 
Indeed,  geostrophic  transports  across  79°N  referenced  to  500  db  between 
about  5°30'-8°30'E  calculated  from  data  obtained  on  cruises  in  Septem¬ 
ber  and  November,  1977  were  0.02  and  1.2  Sv,  respectively,  indicating 
significant  fluctuations  on  a  time  scale  of  only  a  few  months. 
Finally,  application  of  the  objective  error  limits  (see  earlier)  shows 
that  the  expected  error  for  Greisman's  transport  estimate  is  greater 
than  the  estimate  itself.  Nevertheless,  based  on  current  record 
comparison,  it  seems  likely  that  the  transport  in  1971-72  actually  was 
greater  than  in  1976-77. 

3.5  Volume  transport  using  bridged  current  records,  1976-79 

Only  at  mooring  location  C  were  current  meters  maintained  for  all 
three  years  1976-79,  and  therefore  the  three-year  bridged  current 


records  from  mooring  location  C  were  used  to  estimate  volume  transport 
for  the  entire  period.  The  gaps  of  about  4  days  between  the  first  and 
second  years'  data  and  of  9  hours  between  the  second  and  third  years' 


Table  6  Mean  northward  velocity  (cm  s“^)  and  vertical  velocity  dif¬ 
ference  at  mooring  location  C  for  75  days  beginning  9  Sep¬ 
tember  in  respective  years.  Vu  and  Vi  correspond  to  upper 
and  lower  current  meters,  typically  at  100  and  400  m  depth. 


1971* 

1976 

1977 

1978 

vu 

28.6 

7.0 

18.9 

16.6 

Vl 

17.3 

4.6 

12.6 

10.6 

Vvi 

11.3 

2.4 

6.3 

6.0 

vl/vu(75  da) 

0.60 

0.66 

0.67 

0.64 

Vvu(l  yr) 

0.86 

(1976-77) 

0.72 

(1977-78) 

0.71 

(1978-79) 

♦Greisman  (1976) 
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data  were  bridged  by  inserting  the  appropriate  linear  trend  in  the  gap, 
graphically  copying  the  ends  of  each  data  series  halfway  across  the 
gap,  and  then  picking  off  data  values  at  the  proper  sampling  interval. 
Thus,  the  data  which  filled  the  gaps  most  closely  resembled  the  mea¬ 
sured  data  nearest  in  time  on  either  side  of  the  gaps. 

A  regression  relation  between  the  total  transport  and  segment  C 
transport  during  1976-77  was  used  to  estimate  transport  for  all  three 
years  1976-79.  The  regression  relation  was  calculated  after  deleting 
transport  values  during  two  intervals,  one  in  August  1976  and  the  other 
in  late  January/early  February  1977  (Fig.  14),  when  transports  of 
segments  A  and  C  were  in  significant  opposition.  It  was  decided  that 
the  segment  C  transport  was  not  a  good  index  of  the  total  transport 
during  these  two  anomalous  periods.  No  other  such  extreme  reversals 
occurred  at  C  in  subsequent  years.  Omitting  the  anomalous  values,  the 
correlation  coefficient  between  total  transport  and  segment  C  transport 
during  1976-77  was  0.69,  and  the  regression  relation  was: 

total  transport  (Sv)  *  1.09  x  segment  C  transport  (Sv)  +  2.56. 

Bi-weekly  values  of  total  transport  for  1976-79  were  calculated  from 
this  relation,  and  the  results  are  plotted  in  Fig.  16.  For  comparison, 
the  graph  includes  the  first  year's  transport  estimates  which  were 
calculated  using  all  six  current  records. 

3.5.1  Comparison  with  1976-77  transport  estimates 

The  mean  total  transport  of  5.3  Sv  during  1976-77  based  only  on 
the  transport  of  segment  C  is  nearly  the  same  as  the  corresponding 
value  of  5.6  Sv  calculated  by  using  all  segments.  Over  their  common 
span,  the  two  transport  curves  (Fig.  16)  show  moderate  similarity  both 
in  magnitude  and  timing  of  fluctuations.  The  largest  discrepancy 
occurs  during  late  July  and  August  1977  when  the  total  transport  calcu¬ 
lated  from  segment  C  exceeded  that  derived  from  all  six  current  records 
by  about  5  Sv.  The  stick  diagrams  in  Fig.  14  show  that  the  velocity 
increased  markedly  during  this  period  at  mooring  location  C,  clearly 
biasing  upward  any  transport  estimate  based  only  on  segment  C. 
Although  it  appears  that  the  regression  relation  does  yield  approximate 
estimates  of  the  total  mean  transport  and  fluctuations,  the  preceding 
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Figure  16.  Estimates  of  volume  transport  (Sv)  of  the  West  Spitsbergen 
Current  in  1976-79.  Estimates  are  based  on  two-week  aver¬ 
ages  of  bridged  records  of  alongshore  current  measured  at 
mooring  location  C  (78°58'N  7°54'E)  west  of  Spitsbergen. 
Estimates  based  on  6  current  records  at  four  mooring  loca- 
•  tions  in  1976-77  are  included  for  comparison.  Dotted  seg¬ 

ments  indicate  periods  when  the  current  records  at  mooring 
C  were  not  representative  of  the  entire  flow.  .Positive 
value  means  northward  transport.  One  Sv  *  10°  irrs  . 
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comparison  demonstrates  the  importance  of  proper  spatial  coverage  when 
measuring  flow  having  variable  lateral  structure. 

3.5.2  Seasonal  and  interannual  variability 

In  general,  the  plot  of  the  calculated  volume  transport  during 
1976-77  (Fig.  15)  shows  a  maximum  in  late  fall  and  winter  and  then  a 
decrease  through  spring.  This  pattern  agrees  qualitatively  with  find¬ 
ings  by  Timofeyev  (1962)  based  on  hydrographic  data  from  13  different 
years  obtained  from  a  section  along  78°H  from  Barentsburg  westward  to 
the  ice  edge.  Note  however,  that  there  is  also  a  large  peak  in  August 
1976  and  a  minimum  in  March  1977.  Mandel  (1976)  presented  a  seasonal 
pattern,  based  on  several  years  of  data,  of  a  maximum  in  November  (7.7 
Sv)  and  a  minimum  in  April  (3.3  Sv).  Sreisman's  (1976)  calculations 
showed  a  maximum  in  September  with  nearly  the  same  value  in  October  and 
November  1971,  then  a  sharp  decrease  to  a  low  value  in  February,  and 
high  values  in  May  and  April  (Fig.  10).  The  minimum  in  August  1972 
contrasts  with  the  large  value  in  August  1976.  The  basic  seasonal 
pattern  interpreted  from  all  these  data  appears  to  be  of  larger  trans¬ 
ports  occurring  in  fall  and  winter  and  smaller  transports  in  spring. 
This  is  a  pattern  which  corresponds  to  that  of  the  strength  of  atmos¬ 
pheric  circulation  in  the  region;  the  relationships  between  atmospheric 
forcing  and  the  flow  are  treated  in  Section  5. 

Timofeyev's  (1962)  analyses  provided  evidence  for  year-to-year 
variations  In  transport.  Using  hydrographic  sections  along  78°N,  the 
mean  transport  in  1956  was  3.9  Sv,  while  in  1960  it  was  only  2.4  Sv. 
Likewise,  Kislyakov  (1960)  reported  mean  transports  in  several  years 
ranging  from  3.2  Sv  In  1954  to  4.9  Sv  in  1957  (excluding  a  suspiciously 
low  value  of  0.8  Sv  In  1955).  Mandel  (1979)  showed  variations  on  about 
a  three  year  cycle  which  appeared  to  be  superimposed  on  a  much  longer 
cycle  of  at  least  30  years  (Fig.  9). 

In  an  attempt  to  provide  a  consistent  basis  for  the  following 
discussion  of  the  long-term  transports  and  patterns  of  fluctuation, 
only  the  transport  values  computed  by  using  the  regression  relation  are 
used.  For  the  entire  three-year  Interval,  the  transport  averaged  5.3 
Sv,  only  slightly  less  than  the  first  year's  average  of  5.6  Sv. 

The  mean  transports  for  the  first  52  weeks,  the  second  52  weeks, 
and  the  remaining  46  weeks  were  respectively,  5.3,  5.7,  and  5.0  Sv. 
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Thus,  during  1976-79,  there  appeared  to  be  little  Interannual  change  in 
the  mean  transport.  There  appeared  to  be  a  weak  seasonal  variation 
with  maximum  values  in  November-December  and  minimum  values  in  May, 
which  corresponds  generally  with  Timofeyev's  (1962)  results.  Compar¬ 
ing  the  three  year-long  segments  shows  that  the  cycle  was  not  smooth. 
During  the  second  year  the  transport  completed  a  cycle  upon  which 
short-term  variations  with  time  scales  of  6  to  8  weeks  and  magnitudes 
to  +3  Sv  were  superimposed.  In  contrast,  in  late  fall  and  winter  of 
the  first  and  third  year,  the  transport  fluctuated  sharply  over  some¬ 
what  longer  time  scales  and  with  greater  magnitudes  than  during  the 
second. 

Some  caution  should  be  used  regarding  the  transport  values  calcu¬ 
lated  only  from  the  flow  at  mooring  location  C.  Recall  that  during  the 
first  year's  deployment,  the  current  records  from  mooring  C  were  nega¬ 
tively  correlated  with  those  of  mooring  A.  This  means  that  the  flows 
at  locations  A  and  C  tended  to  compensate  each  other.  If  this  is 
generally  true,  inclusion  of  mooring  location  A  current  data  in  the 
transport  calculating  scheme  for  all  three  years  would  have  diminished 
the  magnitudes  of  the  transport  peaks  and  smoothed  out  the  transport 
curve.  Obviously  no  such  tempering  can  occur  using  only  mooring  C  data, 
and  the  results  may  reflect  only  the  temporal  charcteri sties  of  this 
particular  current  core  rather  than  those  of  the  total  flow.  Thus, 
inferences  about  the  total  transport  determined  only  from  the  single 
location  rest  heavily  upon  the  assumption  that  to  some  degree,  the 
magnitude  and  variability  of  the  core  flow  at  mooring  C  reflect  those 
of  the  total  transport. 

In  light  of  the  possible  biasing  of  seasonal  and  interannual  varia¬ 
bility  caused  by  basing  inferences  on  data  from  only  one  mooring  loca¬ 
tion,  it  bears  repeating  that  the  transport  measurements  based  on  all 
four  locations  during  1976-77  did  show  a  distinct  seasonal  pattern 
which  agreed  with  the  results  of  previous  investigators. 


4.  Heat  Transport  of  the  West  Spitsbergen  Current 

4.1  Data  preparation  and  transport  calculation  method 

The  temperature  records  were  used  in  conjunction  with  the  current 
records  to  estimate  the  heat  flux  accomplished  by  the  WSC  during  the 
1976-77  deployment.  As  with  the  northward  volume  transport,  the  mean 
heat  flux  was  calculated  for  each  segment  of  the  entire  section,  and 
these  values  were  then  summed  to  yield  the  total  heat  transport. 

A  representative  temperature  profile  was  arrived  at  by  referring 
to  hydrographic  data  obtained  west  of  Spitsbergen  near  79°N.  The  water 
column  was  considered  in  two  parts,  above  and  below  600  m,  as  was  done 
with  the  velocity.  The  temperature  in  the  upper  part  was  sampled  by 
instruments  typically  located  at  100  m  and  400  m.  The  former  was 
assumed  to  represent  the  temperature  between  the  surface  and  200  m,  and 
the  latter  between  200  m  and  600  m.  This  is  sufficient  in  the  verti¬ 
cally  averaged  sense,  as  is  apparent  from  the  discussion  of  error 
estimation  in  4.2. 

The  temperature  below  600  m  was  inferred  from  the  temperature 
measured  at  400  m.  Figure  17  shows  two  temperature  profiles  obtained 
by  averaging  data  from  different  seasons  over  several  years.  One 
profile  derives  from  temperature  data  between  7°-8°E,  and  the  other 
from  data  between  about  5°-8°E.  While  the  profiles  differ  consistently 
by  about  0.3C,  they  are  remarkably  linear  and  have  nearly  identical 
slopes.  In  each  case  the  temperature  at  600  m  is  about  0.9C  colder 
than  that  at  400  m,  i.e., 

T600  =  T400  "  °*9C* 

From  600  m,  the  temperature  was  taken  to  decrease  linearly  to  -1C  at 
the  depth  of  the  zero  isotach  used  in  the  volume  transport  calcula¬ 
tions.  Inspection  of  several  temperature  profiles  shows  this  is  not 
true  in  every  case.  Some  temperature  profiles  reached  -1C  shallower 
than  the  zero  Isotach,  and  some  were  not  linear,  being  either  concave 
upward  or  downward.  However,  the  assemblage  of  profiles  below  600  m  is 
very  nearly  as  described  by  the  simple  model. 

To  calculate  the  mean  heat  transport  for  1976-77,  the  successive 
14-day  averaged  temperatures  were  used  to  derive  the  vertical  tempera- 
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Figure  17.  Average  temperature  profiles  in  the  West  Spitsbergen  Cur¬ 
rent  near  79°N.  Profile  W  was  obtained  by  averaging  all 
available  temperature  data  between  5°-8°E;  N  by  averaging 
between  7°-8°E. 
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ture  structure.  These  temperature  profiles  were  then  multiplied  by  the 
corresponding  volume  transports  segment-by-segment  to  yield  a  time 
series  of  29  bi-weekly  heat  transport  values.  A  reference  temperature 
of  -0.1C  was  used,  as  in  Aagaard  and  Greisman  (1975).  The  product  of 
density  and  specific  heat  was  assumed  to  be  unity  in  CGS  units. 

4.2  Error  estimation 

Errors  in  the  heat  transport  estimates  arise  both  from  errors  in 
the  volume  transport  estimates  and  from  errors  due  to  the  simplified 
vertical  temperature  distribution.  The  first  of  these  was  discussed 
previously,  and  it  was  found  that  errors  in  the  volume  transport  esti¬ 
mates  may  be  as  large  as  the  estimates  themselves.  The  errors  intro¬ 
duced  by  the  simplified  temperature  representation  depend  on  the  com¬ 
plexity  of  the  actual  temperature  profile.  Comparisons  were  made 
between  the  transport  temperatures  derived  from  two  original  tempera¬ 
ture  profiles  and  those  derived  from  the  simplified  representations  of 
the  same  two  profiles.  Data  were  from  hydrographic  stations  taken  at 
very  nearly  the  same  geographic  location,  but  with  one  in  early  Septem¬ 
ber  1977  and  the  other  in  November  1977.  These  stations  were  selected 
as  worst-case  examples  since  they  showed  extremely  high  gradients. 
Figure  18  shows  these  profiles  and  their  representation  for  heat  trans¬ 
port  calculations.  For  profile  A,  the  difference  in  transport  tempera¬ 
tures  between  the  two  is  only  about  5%,  and  for  profile  B  it  is  about 
9%.  These  discrepancies  are  comparable  to  those  which  would  be 
expected  from  uncertainties  of  0.1C  in  the  readings  of  the  thermistors 
measuring  a  temperature  of  about  2C.  The  conclusion  is  therefore  that 
errors  in  the  volume  transport  estimates  are  principally  responsible 
for  errors  in  the  heat  transport  estimates,  and  that  errors  in  the 
temperature  representation  contribute  less  than  10X.  Errors  in  the 
heat  transport  estimates  may  thus  approach  the  estimates  themselves. 

4.3  Heat  transport  b^  mean  flow,  1976-77 

Figure  15  (p.  40)  shows  a  plot  of  heat  transport  by  mean  flow  of 
the  WSC  averaged  over  two-week  intervals  during  1976-77. 

4.3.1  Mean,  range  of  variation 

The  mean  heat  transport  during  1976-77  was  46  TW  northward,  with  a 
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maximum  of  122  TW  northward  and  a  minimum  of  8  TW  southward  (one  tera- 
watt  (TW)  a  10^W).  Table  7  lists  earlier  estimates  of  mean  heat 
transport  for  the  WSC.  These  estimates  were  derived  by  various  methods 
such  as  direct  measurement,  budget  calculations,  and  geostrophic  calcu¬ 
lations  in  conjunction  with  temperature  data,  and  they  were  generally 
based  on  a  reference  temperature  of  OC.  The  range  is  a  wide  one,  from 
about  18  TW  to  130  TW  northward.  The  maximum,  which  was  the  residual 
needed  to  balance  a  heat  budget,  is  nearly  twice  the  next  largest  value 
and  three  to  four  times  as  large  as  the  rest.  If  it  is  excluded  as 
dubiously  large,  the  average  of  the  remaining  six  estimates  is  34  TW, 
with  which  the  1976-77  estimate  of  46  TW  compares  favorably. 

4.3.2  Seasonal  fluctuations 

The  variable  nature  of  the  heat  transport  during  1976-77  is  evi¬ 
dent  in  Fig.  15.  It  ranges  from  122  TW  northward  in  late  December  1976 
to  8  TW  southward  in  early  March  1977.  The  figure  suggests  that  the 
variability  during  the  first  half  of  the  deployment,  July  1976  -  Feb¬ 
ruary  1977,  is  greater  than  that  during  the  second  half,  and  in  fact 
the  variance  during  the  first  half  is  2.4  times  that  during  the  second 
half.  Furthermore,  the  heat  flux  was  generally  higher  in  late  summer, 
fall  and  winter  than  in  spring  and  early  summer.  Figure  19  shows  the 
bi-weekly  heat  transports  by  spatial  segments  over  the  1976-77  deploy¬ 
ment,  and  it  is  apparent  that  the  same  trends  of  greater  variability 
and  generally  greater  heat  transport  over  the  first  half  hold  for  the 
individual  segments  as  well.  In  addition,  on  two  occasions,  large 
fluctuations  in  the  segment  A  and  C  mean  heat  transports  opposed  each 
other.  Generally  the  heat  transport  of  segment  C  was  greater  than  that 
of  either  segment  A  or  D.  Since  the  areas  assigned  to  A  and  C  were 
approximately  the  same,  this  was  due  to  higher  velocity  and  temperature 
at  mooring  C.  The  same  is  not  true  for  the  heat  transport  of  segment 
D.  Here  the  velocity  was  actually  higher,  but  the  total  cross  section 
assigned  to  segment  D  was  smaller  than  that  assigned  to  either  A  or  C. 
This  biasing  is  aptly  demonstrated  by  comparing  the  mean  heat  trans¬ 
ports  per  unit  area  for  segments  A,  C,  and  D  which  were  respectively 
26,  39,  and  100  in  arbitrary  units.  Thus  if  the  total  flow  had  been 
represented  by  the  velocity  and  temperature  measured  at  mooring  D,  the 
heat  transport  estimates  would  have  been  larger  by  about  a  factor  of 
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Table  7.  Previous  estimates  of  the  mean  heat  transport  of  the  West 
Spitsbergen  Current.  One  TW  =  lO^W. 


Source 


Zaitsev,  et  al.  (1961) 

Mosby  (1962) 

Timofeyev  (1963) 

Vowinckel  and  Orvig  (1970) 
Aagaard  and  Greisman  (1975) 
Mandel  (1978) 

Rudels  and  Andersson  (1982) 


Heat  Transport  (TW) 
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Figure  19.  Bi-weekly  mean  heat  transport  (TW)  of  individual  transport 
cross-sectional  segments  corresponding  to  mooring  loca¬ 
tions  A,  C,  and  D  in  the  West  Spitsbergen  Current  during 
1976?Z7.  Positive  value  means  northward  transport.  One  TW 
=  101ZW. 


three,  again  illustrating  the  importance  of  adequate  determination  of 
the  lateral  current  structure. 

Virtually  all  investigations  of  the  volume  and  heat  transports  of 
the  WSC  reveal  seasonal  fluctuations,  generally  with  similar  timing  and 
characteristics.  Hill  and  Lee  (1957)  interpreted  hydrographic  data 
obtained  during  1949-56  west  of  Bear  Island  along  74°25'N  to  show  that 
the  volume  (and  presumably  also  heat)  transport  of  the  WSC  executed  a 
seasonal  cycle  with  maxima  in  January  and  July  and  minima  in  April  and 
October.  However,  the  scatter  of  the  data  points  used  by  Hill  and  Lee 
(Fig.  3,  p.  109)  leaves  so  much  latitude  in  fitting  a  curve  that 
their  interpretation  is  questionable.  Timofeyev  (1963)  showed  a  pat¬ 
tern  for  the  mean  heat  transport  which  had  a  minimum  of  about  12  TW  in 
April,  a  maximum  of  about  48  TW  in  October  (both  northward),  and  an 
annual  mean  of  31  TW  over  the  time  span  1933-60.  However,  as  Lee 
(1963)  pointed  out,  all  but  three  of  his  eighteen  hydrographic  sections 
fell  in  the  period  May-September  so  that  seasonal  fluctuations  were  not 
really  detectable.  For  the  0-500  m  layer  and  relative  to  0C,  Mandel 
(1978)  determined  a  seasonal  pattern  with  a  minimum  heat  transport  of 
13  TW  in  late  April  sharply  increasing  to  a  maximum  of  46  TW  in  late 
July.  The  mean  was  28  TW  northward.  His  method  involved  linearly 
correlating  the  water  level  and  temperature  at  Barentsburg  with  the 
geostrophic  velocity  and  temperature  of  the  WSC  at  78°N,  and  then  using 
historical  water  level  and  temperature  observations  at  Barentsburg  to 
calculate  mean  monthly  and  mean  yearly  volume  and  heat  transport  esti¬ 
mates  dating  back  to  1949.  Statistically,  his  correlation  was  signif¬ 
icant  at  a  confidence  level  of  95%,  but  it  is  difficult  to  assess  the 
actual  degree  to  which  the  flow  velocity  and  temperature  are  repre¬ 
sented  by  conditions  at  Barentsburg,  which  lies  in  an  inlet  on  the 
west  side  of  Spitsbergen.  Furthermore,  Pattullo,  et  al.  (1955)  showed 
that  at  high  latitudes,  seasonal  steric  effects  may  account  for  nearly 
50%  of  the  observed  change  in  recorded  sea  level  over  a  yearly  cycle. 
Whether  or  not  Mandel  took  this  into  account  is  not  mentioned.  How¬ 
ever,  it  seems  that  even  if  the  accuracy  of  the  absolute  magnitudes 
determined  by  Mandel  may  be  uncertain,  the  seasonal  pattern  may  still 
be  representative.  Greisman  (1976)  calculated  a  heat  transport  of  99 
TW  for  July-November  and  about  50  TW  for  December-June,  a  pattern 
consistent  with  that  reported  by  Timofeyev  (1963)  and  Mandel  (1978). 


In  general,  the  seasonal  pattern  of  heat  transport  of  the  WSC 
appears  to  be  greater  transport  in  fall  and  winter  than  in  spring  and 
simmer,  similar  to  the  volume  transport. 

Figure  15  shows  shorter-term  fluctuations  imposed  upon  the  sea¬ 
sonal  pattern  of  heat  transport  for  1976-77.  Especially  between  July 
1976  -  March  1977,  the  heat  transport  executed  several  sharp  fluctua¬ 
tions  with  peak  values  separated  by  about  10  weeks  and  with  changes  of 
100%  and  greater  in  times  as  short  as  two  weeks.  Neither  of  the 
seasonal  heat  transport  patterns  determined  by  Timofeyev  and  Mandel 
showed  such  short-term  fluctuations,  but  this  is  likely  due  to  averag¬ 
ing  over  several  years. 

Visually  comparing  the  1976-77  volume  and  heat  transport  plots  in 
Fig.  15  shows  that  they  are  unmistakably  correlated,  so  it  appears  that 
at  least  on  time  scales  of  a  few  months  or  so,  fluctuations  in  the  heat 
transport  were  largely  due  to  fluctuations  in  the  volume  transport. 
The  absence  of  major  changes  in  the  temperature  records  simultaneous 
with  the  major  changes  in  the  heat  transport  support  this  conclusion. 

The  temperature  records  themselves  do  show  a  strong  annual  signal 
however,  especially  at  about  100  m  (Fig.  20).  Maxima  of  5-6C  generally 
occurred  in  August-September,  and  minima  of  0  to  -1C  in  February-March. 
At  about  400  m,  the  magnitude  of  the  seasonal  fluctuation  is  much 
smaller  (Fig.  20),  the  temperature  ranging  from  slightly  greater  than 
3C  to  slightly  less  than  1C.  As  well,  the  maxima  and  minima  are  not 
nearly  as  sharp  as  those  at  100  m,  and  so  the  seasonal  signal  is  much 
less  distinct  at  400  m. 


4.3.3  Transport  (integrated)  temperature 

One  way  of  considering  the  seasonal  temperature  fluctuations  of 
the  flow  is  to  calculate  the  transport  temperature  (integrated  tempera¬ 
ture)  by  dividing  the  total  heat  transport  by  the  corresponding  volume 
transport.  The  results  for  1976-77  are  plotted  in  Fig.  21.  The  average 
transport  temperature  for  the  entire  interval  was  about  2.0 C.  Because 
of  the  broad  summer  maximum  spanning  the  first  26  weeks,  followed  by  a 
fairly  sharp  transition  to  a  broad  winter  minimum  spanning  about  26 
weeks,  the  temperature  appeared  to  show  a  decreasing  trend.  Also 
plotted  in  Fig.  21  is  the  transport  temperature  for  segment  C  for  the 
same  interval,  which  had  a  mean  of  about  1.9C  for  the  entire  interval 


Figure  20.  Temperature  records  from  depths  of  92  and  392  m  at  mooring 
location  C  (78°58'N  07°54'E)  during  1976-77.  Fluctuations 
with  periods  less  than  168  hours  have  been  removed  by 
filtering. 
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Figure  21.  Transport  (integrated)  temperature  of  the  West  Spitsbergen 
•  .  Current  (all  cross-sectional  segments),  and  of  cross  sec¬ 

tional  segment  C  only,  during  1976-77.  Transport  tem¬ 
perature  was  calculated  by  dividing  bi-weekly  estimates  of 
heat  transport  by  the  corresponding  estimates  of  volume 
transport. 
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and  showed  the  same  basic  pattern  as  the  transport  temperature  for  the 
entire  flow  section.  If  it  is  assumed  that  the  transport  temperatures 
for  the  flow  section  and  for  segment  C  are  similar  in  other  years  also, 
then  segment  C  can  be  used  to  estimate  the  temperature  history  of  the 
entire  flow  over  1976-79.  The  results  are  plotted  in  Fig.  22,  and  it 
appears  that  the  decreasing  trend  in  1976-77  actually  was  a  part  of  a 
trend  of  slightly  decreasing  temperature  over  the  three  year  interval. 
Three  broad  maxima  occurred  in  autumn  of  successive  years,  and  the 
minima  occurred  in  March -April.  The  correspondence  of  the  transport 
temperature  for  1976-77  and  the  average  transport  temperature  of  the 
flow  through  segment  C  over  the  entire  interval  1976-79  suggests  little 
long-term  temperature  change  in  the  WSC  during  this  three-year  period. 
Because  both  the  heat  transport  and  the  transport  temperature  depend  on 
the  volume  transport,  it  is  not  possible  to  isolate  the  role  which 
water  temperature  changes  play  in  determining  long-term  variations  in 
heat  transport  of  the  WSC.  Mandel  (1979)  claimed  that  there  was  a  2-3 
year  cycle  in  volume  and  heat  transports  between  1949-76  (cf.  Fig.  9). 
However,  the  data  series  used  in  this  present  study  are  not  suffi¬ 
ciently  long  to  examine  this  question. 

4.4  Heat  transport  by  turbulent  processes 

4.4.1  Alongshore  transport,  seasonal  fluctuations 

The  heat  flux  effected  by  mean  advective  processes  may  be  comple¬ 
mented  by  heat  flux  due  to  turbulent  processes,  that  is,  by  the  corre¬ 
lated  fluctuations  of  velocity  and  temperature  about  their  respective 
means. 

For  1976-77,  the  turbulent  heat  transport  was  estimated  by  calcu¬ 
lating  the  products  of  daily  averaged  velocity  and  temperature  fluctua¬ 
tions,  then  averaging  these  daily  products  in  groups  of  14  days,  and 
finally  multiplying  the  14-day  averaged  values  by  the  appropriate  sub- 
area  of  the  transport  cross  section.  The  results  are  shown  in  Fig.  23. 

The  average  northward  turbulent  heat  transport  over  the  entire 
1976-77  interval  was  3  TW,  less  than  a  tenth  of  that  due  to  mean 
advection.  The  maximum  turbulent  heat  transport  of  about  16  TW  north¬ 
ward  occurred  both  in  late  December  1976  and  in  early  March  1977,  and 
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Figure  22.  Transport  (integrated)  temperature  of  cross-sectional 
segment  C  during  1976-79.  Transport  temperature  was  cal¬ 
culated  by  dividing  bi-weekly  estimates  of  heat  transport 
by  the  corresponding  estimates  of  volume  transport. 
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Figure  23.  Bi-weekly  total  and  segment-by-segment  meridional  turbu¬ 
lent  heat  transport  (TW)  of  the  West  Spitsbergen  Current 
during  1976-77.  Values  represent  the  averaged  product  of 
fluctuations  of  temperature  and  alongshore  velocity  mul¬ 
tiplied  by  the  appropriate  area  of  the  transport  cross 
section,.  Positive  value  means  northward  transport.  One 
TW  -  1012W. 
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the  minimum  of  -12  TW  (southward)  occurred  a  month  later  in  early 
April.  Particularly  during  the  first  half  of  the  deployment,  the  mean 
and  turbulent  heat  transports  showed  very  similar  patterns  of  varia¬ 
tions.  Agreement  between  the  two  was  much  poorer  after  the  early  part 
of  January,  about  midway  through  the  deployment  interval,  which  is 
roughly  when  the  temperatures  at  the  moorings  first  approached  or  fell 
below  their  respective  record-length  means.  This  also  corresponded  to 
an  interval  when  the  velocity  records  showed  higher  variability.  For 
example,  in  early  March  1977  the  maximum  in  the  turbulent  heat  trans¬ 
port  coincided  with  the  minimium  in  the  mean  advective  heat  transport. 
During  the  interval  April  27-May  11,  1977,  in  each  of  the  transport 
section  segments  but  especially  at  mooring  C,  negative  velocity  and 
temperature  fluctuations  combined  to  produce  a  large  northward  turbu¬ 
lent  heat  flux.  At  the  same  time  the  mean  flow  showed  very  little 
change.  The  minimum  turbulent  transport  during  the  interval  March  25- 
April  7,  1977  resulted  from  the  combination  of  positive  velocity  fluc¬ 
tuations  and  negative  temperature  fluctuations. 

The  sign  of  the  velocity  fluctuations  was  particularly  important 
in  determining  the  direction  of  the  turbulent  heat  transport  during  the 
1976-77  deployment.  Between  late  February  and  early  April  1977,  the 
velocities  at  the  various  moorings  changed  rapidly  from  negative  to 
positive  values  in  phase  with  each  other  to  produce  the  maximum  north¬ 
ward  and  maximum  southward  turbulent  heat  flux  values. 

Only  during  late  July  and  early  August,  1976  were  the  turbulent 
heat  transport  contributions  from  the  different  moorings  in  significant 
opposition.  The  transport  for  segment  C  was  about  22  TW  southward 
while  that  for  segment  A  was  17  TW  northward.  For  both  segments,  the 
temperatures  were  2.5-3.5C  above  their  respective  means.  However,  the 
velocity  in  segment  A  underwent  a  sharp  positive  change  to  30  cm  s”* 
greater  than  its  mean,  and  at  the  same  time  the  velocity  in  segment  C 
showed  a  negative  change  to  about  25  cm  s'*  less  than  its  mean.  Very 
similar  sharp  velocity  fluctuations  in  January  1977  produced  no  such 
striking  opposition  in  turbulent  heat  flux  because  the  temperatures  in 
segments  A  and  C  stood  very  nearly  at  their  mean  values  and  varied 
little  at  that  time. 

The  year-long  plot  of  alongshore  turbulent  heat  transport  (Fig. 
23)  displays  five  peaks  which  occurred  at  an  interval  of  about  70  days. 
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Whether  or  not  this  signal  is  real  is  difficult  to  determine.  The  only 
other  evidence  of  such  a  long-period  signal  during  1976-77  is  a  varia¬ 
tion  with  a  period  of  about  60  days  which  shows  up  in  the  alongshore 
velocity  at  mooring  D.  Except  for  the  first  two  months  of  the  deploy¬ 
ment,  it  appears  the  variation  of  the  total  turbulent  heat  flux  in  the 
alongshore  direction  was  dominated  by  segment  C.  This  follows  from  the 
fact  that  mooring  C  had  higher  variance  than  moorings  A  and  B,  and 
variance  comparable  to  that  at  mooring  D. 

4.4.2  Cross-isobath  flux 

The  cross-isobath  turbulent  heat  flux  was  -0.17  W  cm-2  (offshore) 
averaged  over  all  six  current  records  during  1976-77.  In  comparison, 
the  alongshore  turbulent  heat  flux  was  1.20  W  cm-2  northward.  The 
standard  deviations  for  the  flux  values  for  each  current  record  were 
roughly  ten  times  the  mean  values  for  both  the  cross-isobath  and  along¬ 
shore  fluxes.  The  largest  offshore  heat  flux,  -1.98  W  cm-2,  occurred 
at  mooring  C.  Other  large  values  were  +1.01  and  +0.75  W  cm'2  respec¬ 
tively  at  moorings  A  and  0  (positive  values  indicate  onshore,  eastward 
turbulent  heat  flux).  One  possible  explanation  is  that  the  average 
local  horizontal  temperature  gradient  had  opposite  signs  at  A  and  C, 
and  this  is  possible  if  A  and  C  were  on  opposite  sides  of  a  velocity 
extremum.  In  this  case,  assuming  in  general  that  geostrophy  obtains, 
the  isotherms  would  slope  oppositely  on  either  side  of  the  velocity 
extreme.  The  current  records  and  hydrographic  data  imply  that  the  WSC 
does  possess  lateral  structure,  so  that  such  a  situation  may  obtain. 
The  location  over  the  shelf  break  and  closer  access  to  the  coastal  flow 
over  the  shallow  shelf  preclude  a  similar  speculation  for  mooring  D. 

The  mean  advective  alongshore  heat  flux  is  roughly  an  order  of 
magnitude  greater  than  the  alongshore  turbulent  heat  flux  and  about  two 
orders  of  magnitude  greater  than  the  cross-isobath  turbulent  heat  flux. 
Thus,  the  turbulent  fluxes  appear  to  have  little  significance  for  the 
WSC  as  far  as  heat  transport  is  concerned. 


5.  Forcing  of  the  West  Spitsbergen  Current 

The  earlier  sections  have  been  devoted  to  describing  the  WSC  and 
quantifying  the  volume  and  heat  transports.  In  this  section  attention 
is  focused  on  dynamic  aspects  related  to  forcing  and  long-term  fluctua¬ 
tions  of  the  flow. 

5.1  Wind  forcing  of  the  mean  flow 

The  Icelandic  low  is  the  major  atmospheric  circulation  feature  of  the 
Green land -Norwegian  sea  region  (Vowinckel  and  Orvig,  1970).  In  particu¬ 
lar,  contours  of  annual  mean  air  pressure  (Fig.  24)  show  a  huge  elon¬ 
gated  low  pressure  trough  which  extends  northeastward  from  Iceland, 
passing  north  of  Norway  and  along  the  northern  periphery  of  the  Soviet 
Union.  The  feature  is  even  more  pronounced  in  January  (Fig.  25),  but 
fades  by  July  (Fig.  26)  so  that  rather  weak  pressure  gradients  prevail. 
Thus,  the  winter  circulation  is  much  more  intense  than  the  summer 
circulation.  The  winter  pattern  usually  persists  until  March,  when 
both  the  Siberian  high  and  the  Icelandic  low  begin  to  weaken  dramati¬ 
cally,  leading  into  the  transition  to  the  summer  pattern.  By  October, 
the  Siberian  high  and  the  Icelandic  low  have  re-formed  and  approach 
their  wintertime  strengths. 

The  April  and  October  air  pressure  distributions  (not  shown)  more 
closely  resemble  the  January  than  the  July  distribution  so  that  the 
Icelandic  low  is  a  significant  factor  for  roughly  8-9  months  of  the 
year.  Consequently,  the  annual  mean  winds  reflect  the  influence  of 
the  low,  and  are  directed  toward  the  southwest  over  the  central  Green- 
land-Norweglan  sea  area. 

As  one  means  of  addressing  atmospheric  forcing  of  the  WSC,  daily 
surface  pressure  data  were  taken  from  National  Weather  Service  (NWS) 
pressure  maps  for  the  northern  hemisphere.  These  were  used  to  investi¬ 
gate  the  relationship  between  velocity  measurements  in  the  WSC  between 
1976-79  and  the  atmospheric  pressure  distribution  over  the  Greenland 
and  Norwegian  seas.  Surface  locations  for  a  5-point  pressure  grid  were 
chosen  at  75°N  15°W,  75°N  0°E,  75°N  15°E,  80°N  0°E,  and  70°N  0°E.  The 
pressure  records  extended  from  July  1,  1976  to  June  30,  1979  and  con¬ 
stituted  a  time  series  1095  days  long.  They  were  used  variously  to 
calculate  time  series  of  simple  pressure  differences,  wind  stress,  and 


Annual  mean  air  pressure  (mbar)  over  the  north  polar 
region.  A  value  of  10  corrresponds  to  1010  mbar,  etc. 
(from  Vowinckel  and  Orvig,  1970). 


curl  of  the  wind  stress.  The  grid  provided  for  the  calculation  of 
north-south  wind  stress  at  75°N  7.5°W  and  75°N  7.5°E  and  of  east-west 
wind  stress  at  77.5°N  0°E  and  72.5°N  0°E.  The  wind  stress,  t  ,  was 
calculated  with  the  square-law: 


t  =  p Cuiui  , 

where  p  is  the  density  of  air,  1.27  x  10“8  gm  cm"8;  C  is  the  drag 
coefficient,  1.3  x  10“^;  and  u  is  the  surface  wind  speed,  assumed  equal 
to  the  geostrophic  wind.  This  last  assumption  is  not  strictly  true 
since  surface  winds  are  proportional  to  geostrophic  winds  by  a  factor 
usually  less  than  unity  and  are  deflected  to  the  left,  typically  15-20° 
(Aagaard,  1969).  However,  the  main  goal  in  comparing  atmospheric  and 
water  flow  data  is  to  investigate  statistical  relationships  and  asso¬ 
ciated  time  scales,  and  these  relationships  are  insensitive  to  the  use 
of  geostrophic  rather  than  surface  winds. 

Table  8  lists  the  three-year  average  of  daily  air  pressure  at  the 
five  grid  points  from  July  1976  to  June  1979.  This  pattern  corresponds 
well  with  the  mean  pressure  distribution:  it  reflects  the  influence  of 
the  Icelandic  low  in  the  southern  and  eastern  areas,  and  the  higher 
pressure  to  the  north  and  west  of  the  central  grid  point  represents  the 
influence  of  the  Greenland  high.  The  mean  large-scale  geostrophic 
winds  associated  with  this  pattern,  and  thus  the  surface  stress,  are 
directed  toward  the  southwest.  This  alone  makes  the  prospect  of  forc¬ 
ing  the  mean  flow  of  the  WSC  by  direct  wind  stress  unlikely,  simply 
because  the  wind  stress  is  in  the  wrong  direction. 

The  major  contributions  to  the  curl  of  the  wind  stress  arise 
mainly  from  the  high  gradient  regions  in  the  north  and  west  of  the 
grid;  the  mean  curl  at  75°N  0°E  is  positive  (cyclonic),  and  this  is 
characteristic  over  most  of  the  Greenland  and  Norwegian  seas. 

Between  July  1976  and  June  1979,  the  wind  stress  curl  averaged  1.4 
x  10"8  dynes  cnT^cm"*.  Assuming  a  zonal  section  860  km  wide  (30°  of 
longitude)  at  75°N,  the  meridional  volume  transport  according  to  the 
Sverdrup  relation: 

[  M,.  dx  =  [v  x  r)/e  , 
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Table  8.  Mean  surface  atmospheric  pressure  (mbar)  derived  from  daily 
•  National  Weather  Service  synoptic  pressure  maps  for  five 

locations  in  the  Greenland-Norwegian  sea  region  between 
July  1976  and  June  1979. 


• 

Location 

Mean  Pressure 
(mbar) 

80°N 

00°E 

1014 

9 

75°N 

00°E 

1010 

70°N 

00°E 

1009 

75°N 

15°W 

1013 

9 

75°N 

15°E 

1010 

was  about  20  Sv  northward  (b  =  6  x  10"^  cm“^s“^  at  79°N).  A  second 
transport  estimate  calculated  similarly  for  the  same  zonal  width  but 
based  on  a  mean  value  of  0.56  x  10"®  dynes  cm~®cm"^  derived  from  wind 
stress  curl  data  provided  by  the  National  Marine  Fisheries  Service 
(NMFS)  yields  about  8  Sv.  Either  figure  can  thus  account  for  the  total 
average  northward  transport  of  about  5  Sv  for  the  WSC  from  1976-79. 
Further  support  for  the  role  of  wind  stress  curl  in  forcing  the  mean 
flow  is  that  during  1971-72,  the  magnitude  of  the  curl  was  2.4  x  10“® 
dynes  cm~®cm"^,  nearly  twice  that  during  1976-79.  While  the  mean 
transport  estimate  of  7  Sv  in  1971-72  is  slightly  larger  than  5  Sv  for 
1976-79,  the  mean  velocity  (unbiased  by  computational  method,  in  con¬ 
trast  to  transport  estimates)  at  approximately  the  same  location  was 
also  nearly  twice  as  large  in  1971-72. 

The  difference  between  the  two  values  for  the  wind  stress  curl 
during  1976-79  stems  in  part  from  the  grid  extent  and  in  part  from  the 
conversion  of  the  NMFS  wind  data  from  geostrophic  to  surface  values. 
The  western  and  northern  gridpoints  used  in  calculating  the  larger 
value  lay  in  the  high  gradient  region  near  Greenland.  A  deviation  of  1 
mbar  in  the  mean  pressure  at  either  point  would  change  the  estimate  of 
the  curl  and  thus  of  the  volume  transport  by  about  17%.  The  smaller 
value  for  the  curl  computed  from  NMFS  data  was  based  on  a  smaller,  more 
local  grid  around  75°N  0°E.  Also,  the  magnitude  of  the  geostrophic 
wind  vector  was  multiplied  by  a  factor  of  0.7  and  rotated  15°  left  to 
yield  surface  winds  (cf.  Bakun,  1973).  Since  wind  stress  is  propor¬ 
tional  to  the  square  of  the  speed,  this  conversion  reduced  the  magni¬ 
tude  of  both  the  wind  stress  and  its  curl  to  0.49  of  their  values  had 
geostrophic  winds  been  used.  Aagaard  (1969)  found  that  for  this  high 
latitude  area,  the  conversion  factor  between  geostrophic  and  surface 
winds  was  about  1  rather  than  0.7.  Transports  calculated  from  NMFS 
wind  stress  curl  data  based  on  a  conversion  factor  of  1  would  thus 
yield  a  value  comparable  to  that  based  on  the  pressure  map  data. 

The  inexactness  of  the  knowledge  of  the  atmospheric  pressure  field 
in  the  region  admits  uncertainties  in  calculations  of  wind  stress  and 
its  curl,  but  even  so  it  is  probably  safe  to  say  that  the  magnitude  of 
the  mean  northward  volume  transports  calculated  in  a  previous  section 
can  be  accounted  for  by  the  wind  stress  curl.  Earlier  work  by  Aagaard 


(1970)  suggested  that  transports  as  high  as  35  Sv  could  be  maintained 
by  the  wind  stress  curl  over  the  Greenland  and  Norwegian  seas. 


5.2  Eastern  intensification  and  vorticit.y  considerations 

Greisman  and  Aagaard  (1979)  hypothesized  that  the  northward  Sver¬ 
drup  transport  might  be  steered  eastward  by  the  steep  topography  in  the 
Greenland  Sea  and  concentrated  into  a  relatively  intensive  feature,  the 
WSC,  as  it  followed  converging  isobaths  which  curved  northward.  This 
is  essentially  a  mechanism  for  forcing  water  of  the  Norwegian  Atlantic 
Current  farther  northward,  which  Metcalf  (1960)  showed  to  be  quite 
broad  at  70°N,  extending  to  nearly  10°W. 

Northward  transport  due  to  cyclonic  wind  stress  curl  has  thus  been 
postulated  as  a  driving  mechanism  for  the  WSC,  but  it  is  important  at 
this  point  to  note  that  this  flow  is  not  a  return  current  in  the  sense 
of  the  Gulf  Stream.  The  Gulf  Stream  maintains  continuity  in  the  North 
Atlantic  by  transporting  back  to  the  north  water  which  was  driven 
southward  in  the  interior  by  anti-cyclonic  wind  stress  curl,  and  vor- 
ticity  considerations  dictate  that  the  return  flow  be  an  intensified 
western  boundary  current.  The  WSC  is  essentially  an  intensified  pole- 
ward  eastern  boundary  current.  Were  it  a  return  flow,  then  those  same 
vorticity  arguments  would  dictate  that  it  appear  as  a  southward  current 
along  the  east  coast  of  Greenland.  Such  a  current  exists,  the  East 
Greenland  Current,  but  it  is  probably  not  simply  a  return  flow  which 
closes  the  flow  in  the  Greenland  and  Norwegian  seas.  The  East  Greenland 
Current  is  driven  in  part  by  continuity  constraints  for  the  Arctic 
Ocean,  but  it  appears  that  buoyancy  forcing  is  also  important  (Grif¬ 
fiths  and  Linden,  1981).  While  the  vorticity  conservation  constraint 
does  not  prohibit  concentrated  eastern  boundary  currents,  return  trans¬ 
ports  must  be  on  tho  westward  side  of  the  ocean  (Stewart,  1964),  and 
any  concentrated  eastern  transports  must  be  directly  driven.  The 
implication  is  that  the  mean  flow  of  the  WSC  is  directly  driven  by  the 
wind  (and  possibly  by  thermohaline  forcing),  in  this  case  by  the  curl 
of  the  wind  stress  rather  than  by  the  stress  itself. 

It  is  useful  to  recall  Stommel's  (1965)  classic  heuristic  explana¬ 
tion  for  asymmetric,  western  intensified  flow  in  the  North  Atlantic  and 
to  comment  on  the  Green  land -Norwegian  sea  circulation  and  the  WSC  in 
that  light.  Stommel  asserted  that  the  negative  wind  stress  vorticity 


tendency  could  not  be  balanced  in  both  the  western  and  eastern  portions 
of  the  ocean  if  the  circulation  were  symmetrical.  Assuming  no  net 
meridional  flow,  the  negative  and  positive  planetary  vorticity  tenden¬ 
cies  due  to  northward  and  southward  moving  water  columns  in  the  west 
and  east  respectively  would  balance,  and  unless  the  velocities  were 
many  times  as  fast  as  observed  in  the  real  ocean,  lateral  friction 
could  not  produce  a  positive  vorticity  tendency  on  the  west  side  of  the 
basin  strong  enough  to  balance  the  negative  wind  stress  vorticity 
tendency.  However,  in  the  case  of  asymmetric,  western  intensified 
flow,  the  frictional  and  planetary  vorticity  tendencies  are  greatly 
enhanced  on  the  western  side  of  the  basin  due  to  increased  shear  and 
velocity  in  the  narrow  northward  flow.  Thus  at  the  western  boundary, 
the  primary  balance  is  between  the  large  negative  planetary  and  posi¬ 
tive  frictional  vorticity  tendencies,  while  the  balance  on  the  eastern 
side  is  between  the  small  negative  wind  stress  and  positive  planetary 
vorticity  tendencies. 

Now  consider  the  Greenland  and  Norwegian  seas  where  the  mean  wind 
stress  vorticity  tendency  is  positive.  On  the  eastern  side,  both  the 
planetary  and  the  frictional  vorticity  tendencies  are  negative  and  will 
be  enhanced  by  the  intense  northward  flow  of  the  WSC,  although  the 
frictional  vorticity  tendency  may  dominate  over  the  planetary  vorticity 
tendency  because  of  the  small  value  for  beta  at  high  latitudes.  Thus 
they  cannot  balance  each  other,  and  since  they  are  together  much 
greater  than  the  wind  stress  vorticity  tendency,  they  cannot  be  bal¬ 
anced  by  it.  The  point  is  that  a  vorticity  argument  analogous  to  that 
for  the  North  Atlantic/Gulf  Stream  system  does  not  yield  a  balance  on 
the  eastern  side  of  the  Greenland-Norwegian  sea  region.  This  in  turn 
implies  that  eastern  intensification  cannot  result  from  a  constraint 
that  the  vorticity  tendencies  arising  from  input  by  wind  stress  curl, 
from  frictional  effects,  and  from  meridional  flow  all  balance.  In 
fact,  eastern  intensification  serves  to  further  destroy  the  possibility 
of  balancing  the  negative  and  positive  vorticity  tendencies,  with  the 
net  result  that  negative  vorticity  would  be  transported  into  the  Arctic 
Ocean  by  the  WSC.  Whether  or  not  this  is  actually  true  will  not  be 
pursued  here,  although  the  shallow  shelves  and  extensive  ice  cover  of 
the  Arctic  Ocean  certainly  have  the  potential  to  dissipate  any  vortic¬ 
ity  that  may  be  transported  in. 
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More  recently,  western  intensification  has  been  understood  in 
terms  of  the  reflection  of  westward-propagating  planetary  waves  at  a 
western  boundary.  Dissipation  of  the  resulting  eastward  propagating 
waves  near  the  boundary  forms  an  intense  boundary  layer  (cf.  Anderson 
and  Gill,  1979).  Thus,  vorticity  waves  generated  in  the  interior  give 
rise  to  a  western  boundary  current.  Over  flat  topography,  dynamical 
effects  constrain  these  free  waves  to  possess  only  westward  phase 
velocity,  and  therefore  this  mechanism  could  not  produce  an  intense 
eastern  boundary  current.  The  inclusion  of  latitudinally-varying  topo¬ 
graphy  can  allow  eastward-propagating  free  planetary  waves  if  the 
topographic  vorticity  effect  opposes  and  dominates  the  planetary  vorti¬ 
city  effect  (LeBlond  and  Mysak,  1978).  This  condition  is  met  in  the 
Greenland  and  Norwegian  seas,  but  the  orientation  of  the  topography  is 
such  that  topographic  waves  would  not  propagate  all  the  way  to  the 
eastern  boundary.  Thus,  reflection  of  planetary  waves  and  subsequent 
dissipation  of  energy  at  the  eastern  boundary  of  Greenland  and  Nor¬ 
wegian  seas  is  highly  unlikely  as  a  mechanism  to  produce  the  observed 
intensification. 

The  most  probable  explanation  for  the  eastern  intensification  is 
the  steering  effect  exerted  on  the  flow  by  bottom  topography.  If  the 
flow  is  constrained  to  conserve  potential  vorticity,  it  will  tend  to 
follow  contours  of  f/H,  where  f  is  the  Coriolis  parameter  and  H  is  the 
depth.  This  is  true  regardless  of  the  type  of  forcing  which  initiates 
the  motion.  In  high  latitude  regions,  f  is  very  nearly  constant, 
changing  only  about  4%  between  70°N  and  80°N,  so  that  the  flow  will 
very  nearly  coincide  with  bottom  contours.  The  major  topographic 
feature  in  the  region  is  the  Mohn  Rise  which  trends  slightly  north  of 
east  from  Jan  Mayen,  separating  the  Greenland  and  Norwegian  sea  basins 
at  about  75°N.  The  typical  depth  of  the  top  of  Mohn  Rise  is  2  km,  and 
its  southern  flank  ascends  northward  with  a  slope  of  roughly  1/100  from 
a  depth  of  3  km  in  the  Norwegian  sea.  Scale  analysis  of  the  barotropic 
vorticity  equation  shows  that  the  effects  exerted  on  the  circulation  by 
topography  and  by  changes  in  latitude  are  comparable  when: 

Hy/H  -  6/f  , 

where  the  subscript  denotes  differentiation,  and  y  is  the  meridional 


coordinate.  The  expressions  for  6  and  f  are: 


6  =  (2^/R)  cos  0  ,  and  f  =  2^sin  0  , 

where  fl  is  the  angular  rotation  of  the  earth,  0  is  the  latitude,  and  R 
is  the  earth's  radius,  " 6.4  x  108  cm.  We  have  at  75°N  for  the  Mohn 
Rise: 

Hy/H  =  ( 1/100 ) / ( ( 3  x  105  cm)  =  3  x  10"8  cm-1  , 

e/f  =  1/(6. 4  X  108  cm  tan  75°)  =  4  x  10‘10  cm"1  . 

The  vorticity  effect  of  the  Mohn  Rise  is  two  orders  of  magnitude 
stronger  than  the  planetary  effect  due  to  latitudinal  changes,  so  that 
a  general  northward  drift  induced  over  the  relatively  flat  Norwegian 
basin  would  tend  to  turn  and  follow  depth  contours  eastward  along  the 
southern  flank  of  Mohn  Rise.  Alternatively,  if  we  define  a  topographic 
scale  length  1  such  that  Hy/H  =  0{1_1)  then  it  follows  that  topographic 
and  planetary  effects  are  equal  for  1  ~  f/o  =  2.4  x  104  km  at  75°N.  In 
practice,  the  scale  length  of  the  topography  is  much  shorter  than  this, 
so  that  topographic  effects  will  in  general  dominate. 

Three  numerical  studies  have  shown  the  potentially  strong  influence 
of  bottom  topography  at  high  latitudes.  Galt  (1973)  ran  a  barotropic 
model  of  the  Arctic  Ocean,  basically  wind  driven,  but  using  only  5%  and 
10%  topography  in  order  to  simulate  the  possible  reduction  of  topo¬ 
graphic  control  by  stratification.  Even  so,  the  model  topography  still 
dominated  the  streamfunction  pattern.  Semtner  (1976)  included  the 
effects  of  stratification  directly  and  obtained  a  streamfunction  pat¬ 
tern  for  the  Arctic  Ocean  similar  to  Galt's.  Creegan  (1976)  modeled 
the  Greenland  and  Norwegian  seas  using  a  two-layer  scheme  which  was 
forced  by  the  wind  and  which  included  the  effects  of  topography  and 
exchange  with  adjacent  seas.  With  flat  topography,  the  circulation 
essentially  reflected  the  wind  field.  With  closed  boundaries  and 
topography  which  she  modeled  as  30%  of  actual  slopes,  the  flow  pattern 
clearly  showed  the  influence  of  the  bottom.  However,  she  obtained  the 
most  realistic  results  when  she  further  included  inflow  through  the 
Faroe-Shetland  Channel  in  the  south  and  allowed  exchange  through  Fram 
Strait  at  the  northern  edge  of  her  model.  None  of  these  three  models 
included  any  thermoha line  forcing. 
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In  summary,  a  plausible  mechanism  has  been  presented  by  which 
atmospheric  forcing  and  topographic  effects  can  give  rise  to  the  WSC. 
Briefly,  it  is  a  northward  Sverdrup  transport  in  the  flat  interior 
being  steered  eastward  along  f/H  contours  by  the  Mohn  Rise,  and  concen¬ 
trated  by  converging  isobaths  which  curve  northward.  It  is  essentially 
an  elaboration  of  the  hypothesis  advanced  by  Greisman  and  Aagaard 
(1979).  Available  atmospheric  and  oceanographic  data  are  insufficient 
to  carry  the  argument  further,  but  there  are  also  other  shortcomings  of 
this  conceptual  model.  One  omission  is  the  role  of  the  Norwegian 
Atlantic  Current,  which  flows  over  the  steep  slope  west  of  Norway. 
This  flow  constitutes  one  portion  of  the  WSC,  and  may  be  affected 
dynamically  by  processes  farther  to  the  south.  Another  consideration  is 
the  degree  to  which  the  response  of  the  ocean  to  the  wind  stress  curl 
is  barotropic.  For  much  of  the  year,  low  stratification  and  weak 
stability,  appropriate  to  a  barotropic  response,  characterize  the 
Greenland-Norwegian  sea  region.  However,  some  barocl inicity  exists, 
and  the  resulting  tendency  toward  confining  the  response  to  the  upper 
ocean  both  diminishes  the  transport  and  weakens  the  effects  of  bottom 
topography.  Nonetheless,  in  light  of  the  earlier  discussion  on  the 
relative  importance  of  topographic  and  planetary  effects,  and  of  the 
numerical  model  results,  it  appears  that  the  barocl inic  effects  on  the 
oceanic  response  are  minor.  A  third  concern  is  the  neglect  of  thermo¬ 
haline  forcing,  and  this  merits  further  discussion. 

5.3  Thermohaline  forcing  of  the  mean  flow 

The  Mediterranean  and  Red  seas  are  two  examples  of  semi-enclosed 
basins  in  which  density  differences  due  to  thermohaline  forcing  drive 
the  circulation.  In  both,  water  flows  in  at  the  surface  and  denser 
water  out  at  depth.  Increase  of  salinity  of  waters  within  the  basin  is 
the  primary  contributor  to  the  density  difference  between  the  inflowing 
and  outflowing  waters. 

Worthington  (1970)  discussed  the  Norwegian  Sea  as  a  similar  medi¬ 
terranean  basin  within  which  inflow  of  surface  water  was  necessary  to 
balance  the  production  and  subsequent  export  of  dense  deep  water.  A 
thermohaline  analog  for  driving  the  mean  flow  of  the  WSC  would  require 
the  production  in  the  Arctic  Basin  of  water  which  would  exit  as  a 
density-driven  outflow  and  which  would  be  replaced  by  a  corresponding 
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inflow  in  order  to  maintain  continuity.  Aagaard,  Coachman  and  Carmack 
(1981)  and  Melling  and  Lewis  (1982)  have  shown  that  cooling  and  salini¬ 
zation  associated  with  the  freezing  of  ice  can  produce  density-driven 
flows  over  the  shelf  and  slope  regions  on  the  perimeter  of  the  Arctic 
Basin  which  then  inject  themselves  into  the  water  column  at  the  appro¬ 
priate  density  level.  Aagaard  et  a_U  (1981)  suggested  this  to  be  a 
mechanism  for  maintaining  the  ubiquitous  halocline  in  the  Arctic  Ocean, 
and  they  estimated  the  production  rate  of  this  water  to  be  about  2.5 
Sv.  The  export  of  such  dense  water  would  represent  a  thermohaline 
forcing  for  the  WSC.  With  respect  to  the  deep  circulation,  Aagaard 
(1981)  differentiated  two  primary  components  of  the  deep  water  in  the 
Arctic  Ocean;  upper  deep  water  less  saline  than  34.92-34.93,  with 
probable  origin  in  the  Greenland  or  Norwegian  Sea,  and  slightly  more 
saline  lower  deep  water  which  has  had  its  salinity  augmented  within  the 
Arctic  Ocean.  Whether  this  lower  deep  water  results  from  brine  rejec¬ 
tion  or  from  cooling  of  saline  Alantic  water  upwelled  near  the  margins, 
it  represents  the  dense  product  of  a  (thermohaline)  water  mass  forma¬ 
tion  process  in  the  Arctic  Ocean,  capable  of  driving  an  exchange  with 
the  Greenland  Sea.  The  presence  of  ice  cover  complicates  the  investi¬ 
gation  of  deep  outflows  on  the  Greenland  side  of  Fram  Strait,  so  it  is 
not  yet  possible  to  verify  a  southward  flow  of  lower  deep  water  out  of 
the  Arctic  Ocean. 

The  discussion  above  of  thermohaline  forcing  has  been  in  terms  of 
the  transformation  of  waters  within  the  Arctic  Basin  and  their  subse¬ 
quent  export.  Stigebrandt  (1981)  developed  an  estuarine  model  for  the 
Arctic  Ocean  based  on  buoyancy  sources  and  sinks  related  to  freshwater 
input  and  removal  through  runoff,  melting  and  freezing  processes, 
vertical  mixing,  and  heat  exchange  with  the  atmosphere.  The  signifi¬ 
cant  result  was  that  by  using  published  data  and  established  parameter- 
izations,  he  was  able  in  the  model  to  balance  realistic  heat,  salt,  and 
volume  budgets  for  the  Arctic  Ocean.  This  demonstrated  the  possibility 
that  estuarine  dynamics  can  account  for  exchanges  comparable  to  those 
actually  observed.  An  additional  feature  of  a  possible  estuarine-like 
circulation  in  Fram  Strait  is  that  the  inflows  and  outflows  are  sepa¬ 
rated  horizontally  rather  than  vertically,  reflecting  the  combined 
influences  of  rotation  and  topography. 
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The  question  of  atmospheric  versus  thermohaline  forcing  of  the 
mean  flow  of  the  WSC  cannot  be  settled  presently.  There  are  simply  too 
few  concurrent  velocity,  hydrographic,  and  climatic  measurements  to 
isolate  and  compare  the  effects  of  both  mechanisms.  Cases  for  both  can 
be  argued,  but  from  a  practical  standpoint,  some  combination  of  both 
probably  operates,  and  moreover,  the  fraction  contributed  by  each  is 
probably  far  from  uniform  over  time. 

5.4  Fluctuations  of  wind  forcing  and  flow 

I  turn  next  to  the  question  of  the  low-frequency  variability  of 
the  flow  and  its  relation  to  variable  wind  forcing.  Stick  diagrams  for 
the  1976-77  current  records  (Fig.  14),  and  three  successive  year-long 
plots  of  the  alongshore  velocity  at  about  100  m  depth  at  mooring  C 
during  1976-79  (Fig.  27)  clearly  demonstrate  that  the  flow  fluctuated 
markedly  even  at  time  scales  as  short  as  a  week. 

5.4.1  Cross-correlations  between  flow  and  atmospheric  paremeters 

To  investigate  possible  connections  between  observed  fluctuations 
of  the  flow  and  changes  in  the  wind  field  as  determined  from  the  large 
scale  daily  pressure  maps,  cross-correlations  were  calculated.  Wi lie- 
brand  (1978)  pointed  out  that  one  deficiency  in  wind  data  derived  from 
hemispheric  pressure  maps  is  that  variability  of  length  scales  1000 
km  cannot  be  defined.  This  is  important  because,  as  shown  by  Wi  lie- 
brand,  et  al.  (1980),  the  oceanic  response  changes  qualitatively  at 
length  scales  of  0(100  km),  below  which  the  response  may  be  strongly 
trapped  near  the  ocean  surface.  Thus  any  correlations  reflect  an 
oceanic  response  to  forcing  on  a  basin-wide  scale,  which  is  about  1000 
km  for  the  Green land-Norweg Ian  sea  region.  Little  if  any  comment  can 
be  made  on  any  smaller  scale  response  based  on  the  pressure  data  used 
here. 

The  bridged  current  record  from  about  100  m  depth  at  mooring  C 
was  chosen  to  represent  the  flow  characteristics  of  the  WSC,  and  cross¬ 
correlations  were  computed  between  the  velocity  components  and  the 
atmospheric  pressure  differences  across  the  basin,  the  wind  stress,  and 
the  curl  of  wind  stress.  This  particular  current  record  was  selected 
because  it  came  from  the  core  region  of  the  flow  and  because  it  was  the 
longest  record  available,  spanning  nearly  three  years.  The  current 
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Figure  27.  Alongshore  velocity  (cm  s'*)  at  100  m  at  mooring  location 
C  during  1976-79.  Fluctuations  with  periods  less  than  168 
hours  were  removed  by  filtering. 
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record  from  about  400  m  depth  at  mooring  C  was  not  included  in  the 
calculations  because,  even  though  speeds  were  about  302  smaller,  it 
closely  resembled  the  shallower  record  and  therefore  would  have  yielded 
similar  cross-correlation  results.  Daily  values  of  current  and  atmo¬ 
spheric  data  were  used,  and  the  time  series  were  slightly  less  than 
1100  days  long.  All  time  series  were  low-pass  filtered  at  several 
cutoff,  frequencies  to  reduce  the  noise  level  in  each,  thereby  smoothing 
the  resulting  cross  correlation  functions  and  rendering  them  easier  to 
interpret.  Specific  details  are  included  in  the  discussion  of  the 
results  for  each  case. 

Correlation  analysis  revealed  no  significant  relation  between 
fluctuations  in  the  cross-isobath  velocity  component  and  any  of  the 
following  atmospheric  parameters;  1)  curl  of  wind  stress  derived  from 
the  5-point  grid,  2)  meridional  component  of  wind  stress  at  75°N  7.5°E, 
and  3)  pressure  difference  between  the  northern  and  southern  grid 
points  on  the  prime  meridian.  Two  sets  of  calculations  were  made,  one 
in  which  the  time  series  were  filtered  with  a  low-pass  cutoff  of  84 
hours  and  a  second  using  a  cutoff  of  360  hours.  Oort  and  Taylor  (1969) 
presented  spectra  of  horizontal  wind  speed  from  six  northeastern  United 
States  locations.  The  spectra  are  typical  of  latitudes  poleward  of 
40°N  and  show  energy  primarily  in  the  2-7  day  range  with  peaks  near  a 
period  of  about  5  days.  The  cutoff  period  of  84  hours  was  chosen  to 
reduce  the  noise  contributed  by  the  high  frequency  end  of  the  spectrum 
and  still  pass  periods  which  contain  most  of  the  synoptic  atmospheric 
variability.  The  cutoff  of  360  hours  was  chosen  in  order  to  exclude 
most  of  the  atmospheric  variability  and  allow  examination  of  longer- 
term  relations  between  wind  and  flow.  In  all  cases,  the  filtered  time 
series  were  subsampled  daily,  and  the  calculations  were  carried  out  to 
400  lags.  Calculations  using  the  higher  frequency  time  series  are  of 
little  value  because  the  potential  forcing  mechanisms  listed  above  are 
essentially  remote,  and  it  is  quite  unlikely  that  atmospheric  fluctua¬ 
tions  could  be  coherent  at  such  short  time  scales  over  hundreds  of 
kilometers. 

The  absence  of  any  strong  correlations  with  the  cross-isobath 
velocity  component  at  longer  time  scales  can  be  attributed  to  the 
rectifying  effect  of  topography  at  mooring  C.  The  energy  spectrum  for 
the  cross-isobath  flow  peaks  broadly  at  a  period  of  about  6  days,  and 
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the  energy  level  at  a  period  of  15  days  is  nearly  an  order  of  magnitude 
smaller.  There  is  little  energy  in  very-low  frequency  fluctuations, 
and  thus  no  reason  to  expect  any  strong  correlation  with  the  wind  at 
corresponding  time  scales.  An  explanation  for  the  cross-isobath  var¬ 
iability  at  time  scales  on  the  order  of  a  week  is  sought  subsequently 
in  Section  6. 

Correlation  analysis  revealed  no  significant  connection  between 
the  north-south  velocity  fluctuations  and  any  of  the  earlier-mentioned 
atmospheric  parameters  when  the  respective  time  series  were  filtered 
passing  only  signals  with  periods  longer  than  84  hours.  However, 
similar  calculations  using  low-passed  data  with  a  cutoff  period  of  15 
days  did  show  some  correlations  (Fig.  28).  For  example,  in  the  case  of 
the  wind  stress  curl  and  the  meridional  flow  at  mooring  C,  the  maximum 
magnitude  correlation  coefficient  is  -0.35  with  95%  confidence  limits 
of  +  0.05,  and  with  the  wind  stress  curl  leading  the  flow  velocity  by 
47  days.  The  negative  (positive)  correlation  means  that  positive 
(negative)  fluctuations  in  the  wind  stress  curl,  and  therefore  positive 
(negative)  fluctuations  in  the  northward  Sverdrup  transport,  are  fol¬ 
lowed  by  negative  (positive)  fluctuations  in  the  velocity  of  the  WSC 
about  seven  weeks  later.  If  the  WSC  represents  Sverdrup  transport,  the 
correlation  therefore  has  the  wrong  sign.  In  fact,  the  wind  stress 
curl  and  the  flow  do  exhibit  a  positive  correlation,  and  it  reaches  a 
maximum  of  +0.26  with  flow  leading  by  about  85  days.  This  implies  that 
the  response  precedes  the  forcing,  which  is  also  inconsistent  with  the 
Sverdrup  transport  hypothesis. 

Correlation  calculations  using  very-low-frequency  passed  data 
(filter  cutoff  corresponding  to  160  day  period)  yielded  a  similarly 
shaped,  but  smoother  correlation  function  (Fig.  28)  with  extremes  of 
+0.69  and  -0.59.  The  cutoff  of  160  days  was  chosen  in  order  to  pre¬ 
serve  fluctuations  with  approximately  semi-annual  and  longer  period 
while  eliminating  the  "noise"  introduced  by  shorter  period  fluctua¬ 
tions.  Note  that  with  the  wind  stress  curl  leading  the  current  (posi¬ 
tive  lags),  the  first  extremum  of  -0.39  compares  favorably  with  the 
value  of  -0.35  obtained  using  15-day  filtered  dnta.  However,  the  lags 
which  are  respectively  77  days  (160-day  filter)  and  47  days  (15-day 
filter)  differ  considerably.  The  explanation  is  that  the  correlation 
at  47  days  is  near  the  beginning  of  a  lag  interval  of  about  80  days 
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Figure  28.  Cross  correlation  between  the  wind  stress  curl  and  the 
alongshore  velocity  at  100  m  at  mooring  location  C  during 
1976-79.  Wind  stress  curl  was  calculated  from  NWS  daily 
synoptic  pressure  maps.  Positive  (negative)  lag  means 
that  curl  leads  (lags)  velocity.  The  dotted  curve  is  from 
very-low  passed  data  with  a  filter  cutoff  of  160  days. 
The  solid  curve  is  from  data  filtered  with  a  cutoff  of  360 
hours.  Dashed  curves  are  95%  confidence  intervals. 


during  which  the  correlation  is  nearly  always  negative  and  of  even 
smaller  magnitude.  The  occurrence  of  the  first  extremum  at  a  77-day 
lag  using  the  160-day  filtered  data  represents  the  effects  of  smoothing 
the  negative  correlation  over  the  entire  80-day  interval  and  of  remov¬ 
ing  the  frequencies  which  introduce  "bumps"  in  the  correlation  func¬ 
tion.  Put  another  way,  for  a  negative  or  positive  correlation  trend 
bounded  by  zero  crossings,  the  maximum  magnitude  is  more  likely  to 
appear  nearer  the  center  of  the  trend  as  higher  frequencies  are 
removed.  In  view  of  the  effects  of  using  filters  with  drastically 
different  cutoff  periods,  the  lack  of  close  agreement  between  the  lag 
periods  of  47  and  77  days  is  not  of  great  significance. 

The  main  point  is  that  in  spite  of  the  great  difference  in  filter 
cutoff  periods,  the  respective  zero  crossings  which  separate  intervals 
of  negative  and  positive  correlation  trends  agree  fairly  well  between 
the  two  correlation  functions.  Thus  both  correlation  functions  reflect 
the  same  basic  pattern,  but  with  the  higher  frequency  data,  the  short 
term  correlations  on  the  order  of  2-4  weeks  tend  to  obscure  the  pat¬ 
tern.  Similar  results  are  obtained  when  similarly  filtered  wind  stress 
data  rather  than  curl  of  wind  stress  data  are  used.  The  conclusion  is 
that  the  correlation  pattern  which  appears  between  the  wind  stress  curl 
and  the  meridional  flow  at  mooring  C  using  160-day  filtered  data  is  not 
just  an  artifact  of  the  filter. 

The  fundamental  pattern  evident  from  Fig.  28  is  that  over  an  800- 
day  interval  corresponding  to  lags  from  -400  days  (velocity  leads  curl) 
to  +400  days  (curl  leads  velocity),  the  correlation  function  calculated 
from  the  160-day  filtered  data  displays  two  cycles  so  that  there  are 
two  negative  extremes,  two  positive  extremes,  and  five  zero  crossings. 
The  three  possible  sets  of  three  successive  zero  crossings,  A-B-C,  B-C- 
D,  and  C-D-E,  define  three  complete  oscillations,  and  the  average 
period  of  these  three  oscillations  is  389  days.  The  intervals  between 
the  successive  negative  extremes,  F-H,  and  positive  extremes,  G-I  are 
respectively  338  and  430  days.  These  intervals  of  389,  338,  and  430 
days  point  to  a  correlation  at  the  annual  time  scale  between  atmo¬ 
spheric  variability  and  the  flow  of  the  WSC.  Similarly,  the  flow  at 
mooring  C  and  the  meridional  wind  stress  at  75°N  7.5°E  (not  shown)  are 
correlated  at  an  annual  time  scale.  Moreover,  the  first  maximum  corre¬ 
lation  of  +0.60  between  these  two  series  occurs  with  the  wind  stress 
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leading  by  86  days,  corresponding  to  a  situation  where  an  increase 
(decrease)  in  the  northward  wind  stress  preceded  an  increase  (decrease) 
in  the  northward  flow  of  the  WSC.  Neither  the  wind  stress  nor  the  wind 
stress  curl  were  correlated  with  the  flow  at  zero  lag.  An  independent 
verification  of  these  results  is  that  calculations  using  monthly  NMFS 
wind  stress  and  wind  stress  curl  data  yielded  correlation  functions 
similar  in  pattern  and  phase  to  those  obtained  above,  albeit  with 
poorer  time  resolution  and  less  significant  extrema.  The  main  point  is 
that  fluctuations  in  the  northward  flow  of  the  WSC  are  correlated  with 
fluctuations  in  both  the  wind  stress  curl  and  the  meridional  wind 
stress  at  the  annual  time  scale. 

However,  the  data  are  rot  sufficient  to  isolate  the  wind  stress 
curl  as  the  sole  forcing  mechanism  for  fluctuations  of  the  WSC  at  the 
annual  time  scale.  Both  the  wind  stress  and  the  simple  meridional 
pressure  gradient  were  correlated  with  the  meridional  flow,  similarly 
to  the  wind  stress  curl.  Clearly,  these  three  atmospheric  parameters 
are  not  independent  and  share  the  same  annual  cycle.  Even  so,  investi¬ 
gations  of  large-scale  ocean  response  to  atmospheric  forcing  are  gener¬ 
ally  made  in  the  context  of  vorticity  input  by  the  wind  stress  curl. 
This  theory  is  examined  later  in  this  section  and  applied  to  the  Green¬ 
land  and  Norwegian  sea  region. 

To  examine  the  stationarity  of  the  relation  between  forcing  and 
response,  cross-correlations  between  seven  successive  overlapping  seg¬ 
ments  of  wind  stress  curl  and  meridional  velocity  at  mooring  C  were 
also  calculated.  The  chosen  segment  length  of  268  days  was  long  enough 
to  grant  confidence  in  results  at  lags  comparable  to  those  of  50-80 
days  at  which  correlations  were  found  using  the  entire  data  series 
lengths  of  about  1100  points.  Calculations  were  carried  out  to  lags  of 
+  100  days.  The  segment  length  of  268  days  was  intentionally  chosen 
shorter  than  annual  to  lessen  the  chance  of  biasing  by  the  connection 
between  wind  stress  curl  and  water  velocity  at  the  annual  time  scale. 
Also,  the  segments  overlapped  by  50%  such  that  the  start  and  end  dates 
progressed  by  134  days  segment-to-segment.  Thus,  no  particular  season 
interval  could  bias  the  results. 

The  cross-correlation  functions  for  the  268-day  segments  of  the 
current  and  wind  stress  curl  data  differ  greatly  among  the  respective 
segments,  indicating  a  high  degree  of  non-stationarity  in  the  relation 
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between  the  two.  Together,  they  show  correlations  ranging  from  -0.71 
to  +0.56  variously  at  lags  in  the  interval  of  0  to  +  100  days.  The 
only  major  correlation,  after  averaging  over  all  segments,  is  -0.30 
with  the  fluctuations  of  the  wind  stress  curl  leading  fluctuations  of 
the  current  by  about  50  days,  corresponding  with  the  results  when  the 
entire  data  series  lengths  of  about  1100  points  were  used.  The  main 
point  is  that  with  time,  the  relation  between  the  wind  stress  curl  and 
the  meridional  flow  at  mooring  C  varies  greatly  at  time  scales  of  100 
days  or  less.  The  important  implication  is  that  fluctuations  in  the 
flow  of  the  WSC  are  in  general  not  consistently  related  to  fluctuations 
in  atmospheric  forcing  in  the  0-100  day  time  range,  and  probably  even 
longer.  The  fluctuations  of  both  in  this  time  range  are  not  jointly 
deterministic,  but  rather  the  short  term  relation  between  them  appears 
more  stochastic,  so  that  there  may  be  only  episodes  where  the  two  are 
related.  It  is  apparent  that  short  term  fluctuations  in  the  wind 
stress  curl  are  not  solely  responsible  for  short  term  fluctuations  in 
the  flow  of  the  WSC. 

5.4.2  Theory 

The  theory  of  oceanic  response  to  large-scale  atmospheric  forcing 
can  provide  a  basis  for  understanding  the  results  of  the  various  corre¬ 
lation  analyses  between  atmospheric  and  oceanic  fluctuations  in  the 
Greenland  and  Norwegian  seas.  The  response  at  sub-inertial  frequencies 
1$  generally  considered  in  terms  of  the  generation  and  propagation  of 
planetary  (Rossby)  waves  (cf.  Anderson  and  Gill,  1975).  The  verti¬ 
cally-averaged  equations  of  motion  describe  the  barotropic  response  of 
the  ocean,  and  from  these  can  be  derived  the  barotropic  vortlcity 
equation  on  a  beta-plane: 

(v2-f2/c2)*t  +  4.x72*y  +  y2i{(x  +  H[l|,x(f/H)y-*y(f/H)x]  *  F  .  (5.1) 

The  derivation  appears  In  many  texts  and  Is  not  repeated  here.  The 
assumptions  are  that  the  flow  is  frictionless  and  quasi -geostrophic. 
The  symbol  Is  the  transport  streamfunction,  H  Is  the  depth,  g  Is  the 
acceleration  of  gravity,  c  Is  the  shallow-water  wave  speed  equal  to 
(gH)%,  and  F  represents  a  general  forcing  terra.  The  zonal  and  meridi¬ 
onal  velocity  components  are  respectively: 


u  «  -(1/H)*y  ,  and  v  ■  (1/H)fx  , 

and  assuming  spatially  uniform  atmospheric  pressure,  the  surface  eleva¬ 
tion  is: 

n  =  (f/c2)'F  . 

Equation  (5.1)  may  be  simplified  further  by  ignoring  the  influence  of 
topography  and  by  ignoring  the  non-linear  terms  under  the  assumption  of 
small  Rossby  number.  For  the  present,  ignoring  topographic  effects  is 
done  as  a  convenience  to  simplify  the  initial  discussion.  It  was  shown 
earlier  in  this  section  that  the  topographic  influence  in  the  vicinity 
of  the  Nohn  Rise  exceeds  the  beta-effect;  possible  effects  of  topog¬ 
raphy  are  discussed  subsequently.  The  assumption  of  small  Rossby 
number  is  valid:  assuming  a  length  scale  comparable  to  the  basin  width 
of  100C  km  and  a  velocity  of  10  cm  s"1,  it  is  of  0(10"3),  and  even 
assuming  a  local  length  scale  of  10  km  for  the  WSC,  it  is  of  0(10_1). 
Equation  (5.1)  then  reduces  to: 

(v2-  '7c2)*t  +  e*x  *  F  .  (5.2) 

To  provide  a  basis  for  general  discussion,  the  left-hand  side  of  (5.2) 
may  be  scaled  to  yield  a  response  time  scale: 

T  -  L/b(L”2  +  x’2)  ,  (5.3) 

where  x2  «  c2/f2  Is  the  square  of  the  radius  of  deformation,  and  L  is 
the  horizontal  scale  of  the  motion.  Taking  L  «  1000  km  and  H  »  3  km 
(c2  *  gH),  (5.3)  yields  a  value  T  -  4  days  which  corresponds  to  a 
period  P  «  2k T  on  the  order  of  a  month.  According  to  Philander  (1979), 
three  frequency  ranges  characterize  the  response  of  the  system 
described  by  (5.2)  for  forcing  at  time  scale  t: 

1)  If  t  «  P,  the  term  B»x  is  negligible,  and  so  in  this  frequency 
range  the  result  is  forced  waves.  The  fluctuations  are  too  rapid  for 
the  beta-effect  to  be  Important. 

2)  If  t  -  P,  both  terms  on  the  left-hand  side  of  (5.2)  are  compa¬ 
rable,  and  winds  which  vary  on  this  time  scale  can  excite  Rossby  waves. 

3)  If  t»  P,  the  first  term  on  the  left-hand  side  of  (5.2)  is 


negligible,  and  the  response  Is  always  In  Sverdrup  balance  with  the 
slowly-varying  wind  stress. 

Thus,  forcing  at  annual  time  scales  will  tend  to  produce  a  slowly 
varying  Sverdrup  response  In  the  ocean  which  will  manifest  Itself  in 
annual  variations  In  the  USC.  Shorter  term  fluctuations  can  be  attrib¬ 
uted  to  local,  higher  frequency  disturbances  and  appear  superimposed  on 
the  slow  Sverdrup  response.  The  response  at  periods  shorter  than  about 
a  month  will  reflect  the  temporal  nature  of  the  forcing,  which  Is 
largely  aperiodic  In  this  frequency  range.  This  can  account  for  the 
absence  of  any  strong  correlation  in  the  calculations  using  the  higher 
frequency  time  series  of  the  wind  stress  curl  and  observed  velocity  of 
the  USC.  Ulllebrand,  et  al.  (1980)  pointed  out  a  second  factor:  the 
scale  of  the  topography  Introduces  other  response  scales  which  can 
obscure  or  destroy  any  clear  correspondence  between  the  wind  stress 
curl  and  the  flow  velocity.  Use  of  a  sufficiently  long  filter  cutoff 
period  for  oceanic  data  should  remove  the  noise  Introduced  by  the 
effects  of  topography  and  nonstationary  Inhomogeneous  forcing,  leading 
to  stronger  correlations  between  oceanic  and  atmospheric  data  series. 
The  results  of  the  correlation  analyses  reflect  this  conclusion,  which 
In  turn  Implies  that  these  two  factors  (and  possibly  others)  Introduce 
stochastic  variability  In  the  USC  time  scales  much  shorter  than  annual, 
perhaps  on  the  order  of  tens  of  days. 

Solution  of  (5.2)  combined  with  physical  parameters  for  the  Green¬ 
land  and  Norwegian  sea  region  allows  a  rough  estimation  of  the  response 
to  large-scale  atmospheric  forcing  at  the  annual  time  scale.  Assuming 
a  wave  like  form  for  the  streamfunctlon: 

t  ■  *0  cos(kx  ♦  ly  -  «t)  , 

substitution  Into  (5.2)  with  F  ■  0  yields  the  familiar  planetary  wave 
dispersion  relation: 


«  »  -ek(k2  ♦  l2  ♦  f2/c2)  ,  (5.4) 

where  « Is  the  angular  frequency  and  k  and  1  are  the  zonal  and  meridio¬ 
nal  wave  numbers,  respectively.  Since  k  and  1  are  inversely  pro¬ 
portional  to  wavelength,  the  phase  velocity,  defined  as«/k,  Is 


directly  proportional  to  the  square  of  the  wavelength.  The  maximum 
wavelength  possible  in  the  Greenland  and  Norwegian  sea  region  is  the 
basin  length  scale  of  about  1000  km  and  thus  the  corresponding  maximum 
phase  speed  of  a  free  barotroplc  planetary  wave  at  75°N  would  be  about 
4  cm  s’1  to  the  west.  The  phase  speed  of  the  first  baroclinic  mode 
would  be  much  slower  still.  Lighthill  (1969)  showed  that  the  length  of 
time  required  to  establish  a  Sverdrup  balance  after  a  sudden  change  in 
wind  forcing  is  equal  to  the  time  It  takes  a  long  Rossby  wave  to 
propagate  from  the  eastern  boundary  to  the  location  of  interest.  For 
the  Ideal  case  mentioned  above,  a  planetary  wave  with  a  phase  speed  of 
4  cm  s"1  would  take  about  five  months  to  propagate  to  the  center  of  the 
Green land -Norwegian  sea  region.  Accordingly,  the  most  rapid  free 
oceanic  response  would  be  confined  to  time  scales  of  this  order. 

However,  as  shown  previously,  the  topographic  influence  of  the 
Mohn  Rise  cannot  be  neglected.  On  both  Its  northern  and  southern 
flanks,  the  restoring  force  due  to  the  bottom  slope  Is  two  orders  of 
magnitude  greater  than  that  due  to  the  beta-effect,  and  therefore  the 
topographic  wave  phase  speeds  would  be  close  to  4  m  s'1,  leading  to 
much  faster  adjustment  on  the  order  of  1.5  days.  In  this  case,  though, 
the  topography  would  support  eastward  phase  velocities  on  the  northern 
flank  and  westward  phase  velocities  on  the  southern  flank,  and  the 
Integrated  effect  across  the  basin  would  be  small.  Furthermore,  the 
lateral  expanse  of  each  of  the  slope  regions  of  the  Mohn  Rise  is  only 
about  a  tenth  of  the  basin  scale,  and  so  any  response  supported  by 
topography  would  be  spatially  limited. 

The  arguments  above  Imply  that  the  Greenland  and  Norwegian  seas  can 
establish  a  Sverdrup  balance  type  of  response  only  if  the  forcing  by 
wind  stress  curl  Is  at  semi-annual  or  greater  time  scales.  This  con¬ 
clusion,  coupled  with  the  Inconsistency  of  an  initial  negative  correla¬ 
tion  with  wind  stress  curl  leading  current  by  about  50-80  days,  sug¬ 
gests  that  fluctuations  In  the  wind  stress  curl  do  not  account  for 
fluctuations  In  the  current  at  time  scales  of  less  than  several  months. 
However,  It  Is  Important  to  reiterate  that  annual  fluctuations  of  the 
wind  stress  curl  can  account  for  the  variability  of  the  WSC  at  the 
annual  time  scale. 


6.  Variability  with  time  scales  of  a  few  weeks  or  less 

The  previous  section  dealt  with  the  seasonal  and  annual  varia¬ 
bility  of  the  flow  of  the  WSC.  The  velocity  records  also  show  a  high 
degree  of  variability  with  time  scales  of  a  few  weeks  or  less.  This 
variability  Is  consistent  In  that  It  repeats  over  record  lengths  of 
several  multiples  of  the  period  of  the  variability.  In  addition,  the 
velocity  records  exhibit  variability  which  is  not  periodic,  that  is,  it 
appears  more  as  isolated  events.  This  section  addresses  mainly  the 
former  variability  and  discusses  the  dynamics  which  may  explain  it. 
Grelsman  (1976)  treated  the  topic  of  tidal  variability  in  the  WSC,  and 
since  the  focus  in  the  present  work  Is  on  sub-tidal  variability,  tides 
are  not  discussed. 

6.1  Spectral  features  of  velocity  and  temperature  records 

Spectra  were  computed  for  the  cross-isobath  (U)  and  alongshore  (V) 
velocity  components  and  for  the  temperature  for  each  data  record  (see 
Table  2  for  mooring  and  data  record  details)  for  each  of  the  three 
years  1976-79  and  for  1980-81.  In  general,  spectra  computed  for  data 
obtained  In  successive  years  from  the  same  location  and  depth  had 
approximately  similar  shapes,  although  there  were  differences  In  their 
finer  scale  features.  However,  spectra  computed  for  each  of  eight 
successive  segments  (about  50  days  long)  of  selected  data  records 
differed  markedly  within  each  given  record.  Under  the  assumption  that 
spectra  computed  from  longer  records  are  representative  of  dynamical 
features  of  the  flow  which  tend  to  be  more  permanent.  It  was  decided  to 
use  the  three-year  bridged  current  and  temperature  records  from  CB-100 
and  CB-400  to  establish  the  general  spectral  characteristics.  This 
decision  Is  also  a  practical  one  which  reflects,  among  other  things, 
the  limitations  imposed  by  the  lack  of  sufficient  spatial  and  temporal 
resolution  of  wind  and  hydrographic  data  to  determine  the  dynamical 
characteristics  of  the  flow. 

I  begin  with  a  brief  description  of  the  spectral  features  of  the 
velocity  and  temperature  measured  at  mooring  location  C. 


6.1.1  U-commonent  spectra 

Figure  29  shows  the  U-component  spectra  for  the  three-year  bridged 
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Energy-conserving  spectra  for  the  cross-isobath  (U)  compo¬ 
nent  of  velocity  at  100  m  and  400  m  at  mooring  location  C. 
Spectra  were  computed  from  three-year  long  direct  measure¬ 
ments.  Smoothing  Is  done  by  averaging  over  equal  log 
frequency  Intervals. 


current  records  CB-100  and  CB-400.  The  spectra  are  very  similar  In 
shape  and  reflect  the  close  correspondence  of  the  temporal  variability 
of  cross-isobath  flow  at  100  m  and  400  m.  For  periods  shorter  than 
about  1000  hours,  the  energy  at  400  m  is  roughly  half  that  at  100  m. 
The  general  shape  of  the  spectra  shows  a  broad  maximum  at  a  period  of 
about  6  days.  The  decrease  of  energy  with  decreasing  frequency  is  a 
result  of  the  constraint  Imposed  by  the  sloping  bottom  topography, 
which  In  steering  the  flow,  damps  cross-isobath  motions.  Spectra  for 
the  U-component  of  velocity  at  locations  A,  B,  C,  and  D  during  the 
1976-77  deployment  show  an  increase  in  low-frequency  energy  with  dis¬ 
tance  offshore.  However,  the  general  peaked  shape  of  the  spectra  is 
maintained,  and  in  no  case  does  the  energy  at  the  low-frequency  end  of 
the  spectra  approach  the  maximum  energy  level.  Although  two  small  peaks 
appear  at  periods  of  6  days  and  4  days  in  both  CB-100  and  CB-400 
spectra,  they  are  not  significant  within  95%  confidence  limits.  The 
significance  of  any  of  the  lower  frequency  peaks  is  questionable. 

6.1.2  V -component  spectra 

Figure  30  shows  the  V-component  spectra  for  the  three-year  bridged 
current  records  CB-100  and  CB-400.  As  with  the  corresponding  U-compo- 
nent  spectra,  the  shapes  of  the  V-component  spectra  are  very  similar  to 
each  other.  In  general,  the  spectrum  for  the  deeper  current  contains 
less  energy,  especially  at  periods  shorter  than  about  6  days,  for  which 
the  V-component  spectra  have  slopes  and  energy  levels  similar  to  their 
respective  U-component  counterparts.  However,  In  contrast  to  the  U- 
component  spectra  which  show  a  decrease  In  energy  with  decreasing 
frequency,  the  V-component  spectra  tend  to  remain  fairly  level  with 
decreasing  frequency.  Since  the  spectra  are  energy  conserving,  this 
means  that  the  energy  is  distributed  fairly  uniformly  over  frequencies 
less  than  that  which  corresponds  to  a  period  of  about  6  days.  While 
close  Inspection  of  the  V-component  spectra  does  show  a  small  tendency 
to  peak,  centered  broadly  at  about  6  days,  it  is  not  statistically 
significant. 

6.1.3  Temperature  spectra 

The  corresponding  temperature  spectra  both  show  marked  peaks  at 
the  annual  period.  While  they  are  not  statistically  significant,  the 
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Figure  30.  Energy-conserving  spectra  for  the  alongshore  (V)  component 
of  velocity  at  100  m  and  400  m  at  mooring  location  C. 

•  Spectra  were  computed  from  three-year  long  direct  measure¬ 

ments.  Smoothing  Is  done  by  averaging  over  equal  log 
frequency  Intervals. 


distinct  seasonal  temperature  cycle  clearly  apparent  In  Fig.  20  (and 
similarly  In  the  two  succeeding  years,  not  shown)  strongly  Implies  that 
the  peaks  are  not  fictitious.  The  temperature  spectrum  from  100  m 
tends  to  stay  fairly  level  with  frequency  until  the  annual  period 
frequency  Is  approached.  There  appear  to  be  no  other  distinct  features 
in  the  temperature  spectrum.  At  400  m,  the  spectrum  contains  about  an 
order  of  magnitude  less  variability  than  that  at  100  m.  However,  the 
400  m  spectrum  does  have  a  weak  maximum  centered  at  about  10  days  and 
extending  over  a  frequency  range  corresponding  to  periods  from  1  to  70 
days. 

6.1.4  Spectra  of  vertical  differences  of  velocity  and  temperature 

Spectra  were  also  computed  for  three-year  bridged  time  series  of 

the  velocity  and  temperature  differences  between  100  m  and  400  m  depth 
at  mooring  location  C.  The  only  notable  feature  of  the  temperature 
difference  spectrum  is  a  strong  peak  at  the  annual  period.  Both  the  U- 
component  and  V-component  shear  spectra  are  nearly  level  and  coincident 
for  periods  shorter  than  about  10  days.  At  lower  frequencies,  to  about 
the  annual  period,  the  U-component  shear  spectrum  decreases  linearly  by 
almost  two  orders  of  magnitude.  The  V-component  shear  spectrum  reaches 
a  minimum  at  a  period  of  about  80  days  then  increases  with  decreasing 
frequency,  with  a  peak  at  the  annual  period. 

6.1.5  Nonstationaritv 

Spectra  computed  for  successive  67-day  segments  of  the  bridged 
velocity  records  display  marked  differences.  For  each  segment,  the 
mean  and  standard  deviation  of  the  smoothed  spectral  estimates  were 
also  computed  for  each  respective  frequency  band.  Generally,  the  stand¬ 
ard  deviation  Is  about  the  same  as  the  mean.  A  plot  of  these  means 
reproduces  the  shape  of  the  spectra  which  were  computed  for  the  record- 
length  data  series.  The  basic  conclusion  is  that  the  statistics  of  the 
velocity  time  series  are  not  stationary,  which  in  turn  implies  that  the 
dynamics  of  the  flow  change  with  time.  This  is  another  reason  for 
investigating  the  dynamics  by  using  the  statistics  of  the  maximum 
record  lengths  available. 

The  U-component  spectra  offer  the  primary  evidence  of  variability 
at  time  scales  of  a  few  weeks  and  less.  The  absence  of  significant 


spikes  (except  tidal)  indicates  that  such  variability  is  not  confined 
to  a  narrow  range  of  frequencies  but  rather  is  distributed  in  a  broad 
peak  reaching  a  maximum  at  a  period  of  about  6  days.  Small  sharp  peaks 
jut  above  this  general  spectral  shape  at  periods  of  6  days  and  4  days. 
These  two  sharp  peaks  coincide  in  period  with  small  peaks  in  the  tem¬ 
perature  spectrum,  which  implies  some  correlation  between  cross-isobath 
flow  and  temperature,  if  the  peaks  are  real.  Neither  the  temperature 
spectra  nor  the  V-component  spectra  reflect  the  broad  U-component  peak 
at  6  days.  Although  the  V-component  spectra  contain  energy  levels 
similar  to  those  of  the  U-component  spectra  at  6  days,  the  V-component 
spectra  remain  nearly  level  with  decreasing  frequency. 

6.2  Instability  mechanisms 

I  next  consider  dynamical  explanations  for  variability  of  the  flow 
at  a  time  scale  of  about  6  days.  The  first  of  these  are  instabilities. 

6.2.1  Barotropic  instability 

A  characteristic  of  the  current  records  is  that  both  velocity 
components  are  strongly  correlated  in  the  vertical.  It  is  sensible, 
then,  to  examine  first  whether  or  not  a  barotropic  instability  mecha¬ 
nism  can  account  for  the  variability  at  time  scales  of  about  a  week. 
Tareyev  (1971)  investigated  the  stability  of  a  laterally  sheared  baro¬ 
tropic  flow  over  sloping  topography.  He  found  that  for  Instability, 
the  total  gradient  of  the  vorticity  (6); 

6  «  (f  -  vx)(Hx/H)  -  vxx  ,  (6.2.1) 

must  change  sign  within  the  lateral  bounds  of  the  flow,  i.e.  within  the 
flow,  G  »0  at  some  point.  Here,  f  is  the  Coriolis  parameter,  v  the 
velocity,  H  the  depth,  and  the  subscripts  denote  differentiation  (x  is 
cross-isobath).  Scaling  by  a  characteristic  velocity  V  and  length  L 
yields: 

(f  -  V/L ) (H/LH)  -V/L2  »  0  .  (6.2.2) 

Slmpllfyl  •  (6  t)  yields  the  condition  that  the  Rossby  number  V/lf 

should  be  neater  than  0.5.  For  the  WSC,  V  -  0.1  m  s-1  and  L~  100  km, 
so  with  f  -10’*  s*1  the  Rossby  number  is  0.01,  and  barotropic  insta- 


billty  is  unlikely.  Even  assuming  a  length  scale  of  20  km,  comparable 
to  the  width  of  one  of  the  velocity  cores,  and  a  velocity  of  0.6  m  s"* 
corresponding  to  the  largest  observed  lateral  velocity  difference,  the 
Rossby  number  is  0.3,  so  that  instability  is  still  only  a  marginal 
possibility. 

The  scaling  argument  is  only  a  guideline,  so  the  question  of 
whether  the  basic  conditions  for  instability  obtain  is  explored  further 
by  using  a  simple  model.  Niiler  and  Mysak  (1971)  used  a  model  to 
investigate  the  barotropic  stability  of  a  laterally  sheared  flow,  and 
they  applied  it  to  the  Gulf  Stream.  Mysak  and  Schott  (1977)  used  the 
same  model  (model  b  of  Niiler  and  Mysak,  1971)  in  connection  with  an 
investigation  of  variability  of  the  Norwegian  Current  with  time  scales 
of  2-3  days.  Basically,  the  physical  model  is  of  a  barotropic  triangu¬ 
lar  current  profile  (maximum  velocity  at  the  middle)  centered  over  the 
discontinuity  of  a  single  step  bottom  topography  configuration  with 
lateral  boundaries  at  infinity.  The  depth  of  the  lower  step  is  H2  and 
the  upper  Hj_.  The  stability  of  the  triangular  profile  with  maximum 
velocity  V  is  considered  with  respect  to  small  perturbations  of  the 
form: 

f  a  y(x)  exp(ily  +  iwt)  , 

where  vis  the  transport  stream  function,  1  is  the  wavenumber,  andu 
is  the  wave  frequency.  Niiler  and  Mysak  (1971)  derived  the  dispersion 
relation  a*  a(k)  according  to  linear  hydrodynamic  theory,  and  the 
details  are  not  be  repeated  here.  The  dispersion  relation  (their 
equation  8)  is: 

(  o+  R  exp(-k)sinh  k)  (ao^  +  bo  +  c)  *  0  ,  (6.2.3) 

where  a  »  1  ♦  {H2/Hj  , 

b  *  aR(k  -  1  +  exp(-k)cosh  k)  +  2  -  a  , 
c  *  R  exp(-k)  (aR(k  cosh  k  -  sinh  k)  +  (2  -  a)sinh  k)  , 

the  non-dimensional  frequency  a  is  the  dimensional  frequency  u>  divided 
by  f,  the  non-dimensional  wavenumber  k  is  the  dimensional  wavenumber  1 
multiplied  by  the  current  half-width  L,  and  R  is  the  Rossby  number 
V/Lf.  As  Niiler  and  Mysak  (1971)  pointed  out,  the  first  factor  of 


(6.2.3)  describes  a  stable  wave  which  propagates  over  the  depth  discon¬ 
tinuity.  For  a  given  k,  the  second  factor  yields  two  roots  which  may 
be  either  real  or  complex  conjugates.  Complex  conjugate  roots  denote 
instability  where  the  wave  frequency  is  or,  the  real  part  of  the  com¬ 
plex  frequency,  and  where  one  of  the  waves  grows  at  a  rate  proportional 
to  exp(a.jt)  (o^t  is  the  imaginary  part  of  the  complex  frequency). 
Figure  31  shows  the  dispersion  curves  calculated  for  two  cases  for  the 
WSC  using  (6.2.3).  the  parameters  used  for  case  A  correspond  to  treat¬ 
ing  the  current  as  a  single  broad  flow;  L  (the  current  half  width)  =  75 
km;  Hj  (depth  of  the  shelf)  =  250  m,  Hg  (depth  at  about  5°E)  =  2500  m, 
V  *  15  cm  s’1,  and  f  *  0.000142  s"1.  The  small  Rossby  number,  0.014, 
would  imply  stable  flow,  and  the  results  in  Fig.  31a  support  this.  The 
upper  curve,  with  frequency  which  increases  slowly  with  wavenumber  from 
an  initial  value  of  o=  0.82,  corresponds  to  a  double-Kelvin  wave 
modified  by  the  flow.  The  two  lower  curves  coincide,  both  of  them 
having  a  non-dimensional  frequency  a  =  0.007  for  k  >  1.3.  All  three 
waves  are  stable,  and  their  periods  of  about  1  day  and  72  days  do  not 
correspond  with  major  spectral  features  of  the  WSC.  Similar  calcula¬ 
tions  using  V  *  50  cm  s’1  also  failed  to  identify  any  unstable  waves  or 
any  stable  waves  with  periods  prominent  in  the  fluctuations  of  the 
current.  This  is  notable  because  50  cm  s"1  represents  an  extreme 
velocity  estimate  for  the  WSC,  and  therefore  a  maximum  likelihood  for 
instability. 

The  second  case,  B,  is  representative  of  the  stability  at  local 
scales  characteristic  of  the  lateral  structure  of  the  flow.  The  para¬ 
meters  were  L  *  20  km,  *  250  m,  Hg  s  1250  m,  and  V  *  15  cm  s’1;  the 
computed  Rossby  number  was  0.053.  Figure  31b  shows  the  dispersion 
curves  for  this  case.  Again,  no  instability  is  indicated.  Two  roots 
coincide  at  a  *  0.026  for  k  >  2,  while  the  root  corresponding  to  the 
third  wave  grows  only  slightly  from  an  initial  value  of  o*  0.67.  In 
addition,  neither  of  the  periods  of  about  1  and  19  days  associated  with 
these  calculated  stable  waves  corresponds  with  prominent  fluctuations 
in  the  flow,  and  the  periods  differ  markedly  from  the  6  days  observed 
in  the  broad  peak  in  the  U-component  spectra.  Using  a  value  of  50  cm 
s’1  for  V  (Rossby  number  »  0.176)  also  suggested  no  instability.  How¬ 
ever,  one  of  the  stable  waves  did  have  a  period  of  about  6  days  for  a 
wavelength  less  than  about  100  km.  As  evidenced  by  the  current  records 
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Figure  31.  Dispersion  diagram  for  the  Hi  Her  and  Mysak  (1977)  baro- 
tropic  instability  model  for  two  di77erent  current  widths 
(L)  and  two  different  velocities  (V).  a)  Case  A:  L  *  75 
km,  V  *  15  cm  s’1  (Rossby  number  0.014),  V  *  SO  cm  s’1 
(Rqssby  number  0.047).  b)  Case  B:  L  *  20  km,  V  *  15  cm 
s’1  (Rossby  number  0.053),  55  cm  s"1  (Rossby  number 

0.176).  Non-dimensional  frequency  o  is  the  dimensional 
frequency  <*»  divided  by  f.  Non-dimensional  wavenumber  k  is 
the  dimensional  wavenumber  1  multiplied  by  L. 
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at  A  and  C,  a  lateral  velocity  difference  of  50  cm  s-*  was  very  infre¬ 
quent,  and  so  it  is  unlikely  that  this  result  has  any  significance  in 
accounting  for  any  major  features  of  the  variability  of  the  WSC. 

In  summary,  barotropic  instability  is  not  a  major  contributor  to 
the  variability  of  the  WSC,  as  determined  from  results  of  scaling 
considerations  and  of  a  model.  The  mean  velocities  and  horizontal 
shears  are  too  small  to  produce  this  kind  of  instability  given  the 
length  scales  characteristic  of  the  current. 

6.2.2  Baroclinic  instability 

Perturbations  due  to  this  type  of  instability  draw  energy  from  the 
available  potential  energy  stored  in  the  water  column,  expressed  by 
tilted  isopycnals.  Buoyancy  fluxes  effect  the  transfer  of  this  poten¬ 
tial  energy  into  kinetic  energy  of  the  perturbations.  A  simple  way  of 
expressing  the  situation  is  that  a  flow  may  be  baroclinically  unstable 
when  the  vertical  velocity  shear  is  greater  than  the  geostrophic  velo¬ 
city  based  on  the  sloping  isopycnals.  The  buoyancy  flux  acts  to  level 
the  Isopycnals.  One  Indication  of  a  buoyancy  flux  in  the  proper  direc¬ 
tion  to  feed  energy  into  perturbations  is  that  the  mean  lateral  turbu¬ 
lent  heat  flux,  u'T',  is  down-gradient.  For  the  WSC,  iPT1,  as  deter¬ 
mined  from  the  current  and  temperature  records,  is  negative  (down- 
gradient),  and  so  Is  at  least  in  the  proper  direction  to  transfer 
energy  from  the  mean  flow  into  perturbations. 

Details  of  the  derivations  of  necessary  conditions  for  the  occur¬ 
rence  of  baroclinic  instability  may  be  found  in  LeBlond  and  Mysak 
(1978)  and  Gill  (1982)  among  others,  and  only  the  result  is  summarized 
and  discussed  here.  The  conditional  equation  is: 


2  2 
|  Qy  [|Z|/!U-c|]  dz  +  (f/N)2  [  |  Z  |  / 1 U  -  c  |  j  U2 


(6.2.4) 


where  Z  Is  a  complex  amplitude  of  the  perturbation  streamfunction,  c  is 
the  complex  phase  speed,  f  is  the  local  value  of  the  Coriolis  parame¬ 
ter,  N  Is  the  buoyancy  frequency,  Uz  is  the  vertical  gradient  of  the 
velocity  of  the  mean  field,  Qy  »eT  -  ((f/N)^uz)z  is  the  cross-isobath 
gradient  of  mean  potential  vorticity,  and  sT  is  the  topographic  beta- 
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effect.  Thus  it  is  necessary  that  either  the  first  or  the  second  term 
of  (6.2.4)  be  negative  for  instability  to  be  possible.  As  enumerated 
by  Wright  (1981)  there  are  four  possibilities:  1)  Qy  changes  sign 
within  the  fluid;  2)  somewhere  in  the  fluid,  the  sign  of  Qy  is  opposite 
to  that  of  Uz  at  the  surface;  3)  somewhere  in  the  fluid,  Qy  has  the 
same  sign  as  that  of  Uz  at  the  bottom;  and  4)  the  sign  of  Uz  is  the 
same  at  the  surface  anfl  the  bottom.  Errors  in  the  observed  vertical 
profiles  of  Qy  and  U2  may  lead  to  incorrect  conclusions  as  to  whether 
or  not  any  of  the  above  possibilities  is  actually  fulfilled.  In  order 
to  argue  with  any  confidence  the  question  of  whether  or  not  baroclinic 
Instability  is  possible  in  terms  of  necessary  conditions,  it  is  thus 
very  Important  that  the  vertical  profiles  of  Qy  and  Uz  be  accurate. 
The  data  set  for  the  WSC  is  insufficient  to  establish  the  required 
profiles  of  Qy  and  Uz*  and  so  no  attempt  can  be  made  here  to  a  priori 
predict  whether  or  not  the  flow  is  unstable  in  the  baroclinic  sense. 

Schott  and  Bock  (1980)  examined  low-frequency  fluctuations  of  the 
Norwegian  Current  from  the  standpoint  of  energy  Interaction  terms. 
They  calculated  vertical  profiles  of  barotropic  and  baroclinic  insta¬ 
bility  terms,  and  they  used  them  in  conjunction  with  a  kinetic  energy 
budget  equation  for  fluctuating  motions  to  determine  the  magnitude  and 
direction  of  transfer  of  energy  between  the  mean  and  fluctuating 
fields.  In  order  to  do  so,  they  required  extensive  hydrographic  meas¬ 
urements  as  well  as  current  measurements  from  a  diamond  array  of  moor¬ 
ings.  A  similarly  comprehensive  data  set  is  not  available  for  the  WSC, 
so  that  it  is  not  possible  to  address  the  baroclinic  (and  barotropic) 
Instability  problem  from  energy  transfer  considerations,  either. 

Another  approach  is  to  consider  a  simple  baroclinic  model  of  the 
WSC  and  to  Investigate  the  stability  of  traveling  wave  perturbations  in 
the  model.  The  model  used  Is  of  a  two-layer  geostrophlc  flow  in  a 
channel  with  a  sloping  bottom.  Smith  (1976)  used  this  configuration  to 
show  that  baroclinic  Instability  could  explain  fluctuations  in  the 
Denmark  Strait  overflow  with  a  period  of  about  2  days.  Mysak  and 
Schott  (1977)  applied  this  configuration  to  explain  the  presence  of  a 
spectral  peak  at  2-3  days  In  the  flow  of  the  Norwegian  Current.  In 
view  of  the  positive  results  obtained  In  these  two  Investigations  and 
of  the  similarity  between  the  physical  settings  for  the  WSC  and  the 
Norwegian  Current,  this  model  likewise  is  applied  to  the  WSC. 
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Mysak,  Johnson,  and  Hsieh  (1981)  subsequently  examined  the  valid¬ 
ity  of  using  a  channel  model  to  investigate  instabilities  of  coastal 
currents  in  cases  where  there  is  actually  an  open  boundary  on  the 
seaward  side  and  a  shelf  leading  to  a  coastal  wall.  They  found  that 
either  an  outer  boundary  consisting  of  a  quiescent  ocean,  or  an  inner 
boundary  of  a  shallow  sloping  region  rather  than  a  rigid  wall,  gave 
rise  to  a  new  class  of  baroclinic  instabilites  due  to  interaction 
between  waves  of  differing  cross-stream  structure.  Even  so,  the  domi¬ 
nant  baroclinic  instability  was  that  associated  with  the  gravest  mode. 
In  the  WSC  case,  only  the  gravest  mode  instability  is  investigated. 

I  describe  the  model  briefly  below  and  refer  the  reader  to  Mysak 
and  Schott  (1977)  for  details  of  the  derivation  of  the  dispersion 
relation  for  waves  in  the  flow.  Basically,  the  derivation  seeks  tra¬ 
veling  wave  solutions  of  the  linearized,  quasi -geostrophic  potential 
vorticlty  equations  for  each  of  the  two  layers,  subject  to  the  condi¬ 
tion  of  no  normal  flow  at  the  channel  walls.  The  resulting  dispersion 
relation  is: 


acn  +  ba  +  c 


(6.2.5) 


where  a  ■  K2(K2  +  ♦  F2), 

b  -  -k((Vx  +  V2)a  ♦  K2(Vx  -  V2)(F2  -  Fx)  ♦  T(K2  ♦  F^), 

C  »  k2(aV1V2  +  K 2(VX  -  V2)(V1F2  -  V2FX)  +  (K2  ♦  Fx)). 

For  the  model,  the  width  of  the  channel  is  L,  the  upper  and  lower  layer 
thicknesses  are  h^  and  h2  respectively,  and  the  bottom  slope  is  s, 
positive  upward  toward  the  shoreward  wall.  The  horizontal  velocity 
scale  is  U,  and  reduced  gravity  is  denoted  by  gr.  The  other  parameters 
are  Fj,  F2  «  f2L2/grhj^,  f2L2/grh2,  the  internal  Froude  numbers  for  each 
layer;  T  «  sfL2/Uh2,  a  topographic  parameter;  Vj,  V2  are  non-dimen¬ 
sional  velocities  equal  to  V]/U,  v2/U  where  vj  and  v2  are  actual  velo¬ 
cities  in  the  respective  layers;  k  is  the  dimensional  wavenumber  1 
multiplied  by  L;  and  K2  ■  k2  ♦  rm2,  a  wavenumber  parameter  where  n  is 
the  mode  number.  The  non-dimensional  frequency  is  o*  «/(U/L),  the 
frequency  divided  by  an  inverse  time  scale.  Solving  the  dispersion 
relation  (6.2.5)  yields  two  values  for  a  which  may  be  real  or  complex 
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conjugates.  The  latter  case  corresponds  to  a  baroclinic  instability, 
and  one  of  the  waves  will  grow  at  a  rate  proportional  to  exp(o^t).  On 
a  dispersion  diagram,  the  line  along  which  the  two  roots  coalesce  indi¬ 
cates  the  frequency-wavenumber  range  over  which  instability  is 
possible. 

Calculations  were  made  for  two  channel  widths,  75  km  (Case  A)  and 
20  km  (Case  B),  corresponding  to  current-width  and  core-width  length 
scales,  respectively.  In  each  case,  the  upper  layer  depth  h^  was  taken 
as  250  m,  based  on  vertical  density  profiles  which  showed  the  strong 
pycnocline  above  100  m  to  have  decayed  away  for  the  most  part  by  250  m. 
For  case  A,  the  lower  layer  depth  h2  was  2000  m,  typical  of  the  depth 
of  the  outer  slope  region,  and  for  case  B,  h2  was  1000  m,  typical  of 
the  depth  In  the  region  of  moorings  A  and  C.  In  each  case,  calcula¬ 
tions  were  made  using  two  velocity  sets:  one  set  was  (U,vj,v2)  * 

(25.15.4)  cm  s*1,  and  the  second  was  (25,25,4)  cm  s"*.  Results  depend 
only  on  the  values  of  v^  and  v2,  and  using  a  common  scaling  velocity  U 
allows  results  to  be  plotted  on  the  same  o-k  diagram.  The  value  of  4 
cm  s'1  for  v2  was  determined  from  the  direct  current  data  from  C5-984 
(1980-81).  The  value  of  15  cm  s"1  for  v^  Is  typical  of  the  mean 
current  in  the  upper  250  m.  The  value  of  25  cm  s~*  for  vj  was  chosen 
in  order  to  investigate  the  effect  of  larger  vertical  shear.  A  velo¬ 
city  histogram  from  C5-84  Indicates  that  about  80#  of  the  speed  meas¬ 
urements  were  25  cm  s’*  or  less  (14#  In  the  range  20-25  cm  s"*),  so 
that  25  cm  s’*  seems  a  good  choice  for  an  upper  limit  which  might  still 
be  expected  to  occur  a  significant  part  of  the  time.  Reduced  gravity 
was  0.003  m  $"*,  determined  from  density  profiles,  and  the  bottom 
slopes  were  0.031  and  0.010  respectively,  for  cases  A  and  B. 

Figure  32a  shows  the  dispersion  curves  for  case  A  with  (vj,v2)  * 

(15.4)  cm  s'1  and  with  larger  shear  (25,4)  cm  s“*.  Baroclinic  insta¬ 
bility  Is  possible  for  both  sets  of  values.  For  the  smaller  shear,  the 
wave  with  the  maximum  growth  rate  (indicated  by  the  dot)  has  a  wave¬ 
length  of  34  km  and  a  period  of  3.8  days.  The  range  of  periods  for 
Instability  Is  about  3.4-4  days.  For  the  larger  shear  the  maximum 
growth  rate  is  at  wavelength  41  km  and  period  3.4  days.  The  instabil¬ 
ity  range  Is  about  2.9-3.9  days.  Thus,  while  doubling  the  value  of  the 
shear  Increases  the  range  over  which  Instability  may  occur.  It  does  not 
greatly  affect  the  wavelength  and  period  of  the  maximally  unstable 
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Figure  32.  Dispersion  diagram  for  the  Hvsak  and  Schott  (1977)  baro- 
cllnlc  Instability  model  for  two  different  current  widths 
(L)  and  two  different  velocity  pairs  (vi,v?).  The  pair 
(vi.v?)  are  the  upper  and  lower  layer  velocities,  respec¬ 
tively.  Dot  denotes  maximum  growth-rate  wave.  The  scale 
velocity,  U,  Is  25  cm  s'1,  a)  Case  A;  L  ■  75  km;  bottom 
slope,  s  *  0.031;  (vi,v?)  as  shown,  b)  Case  B:  L  *  20 
km,  s  ■  0.01.  Non-dimensional  frequency  o  is  the  dimen¬ 
sional  frequency  u divided  by  U/L.  Non-dimensional  wave- 
number  k  is  the  dimensional  wavenumber  1  multiplied  by  L. 
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wave.  However,  the  maximum  growth  rate  associated  with  the  larger 
shear  Is  about  twice  that  associated  with  the  smaller  shear. 

The  periods  mentioned  above  are  shorter  than  the  period  of  about  6 
days  corresponding  to  the  broad  peak  In  the  U-component  spectra.  The 
wavelengths  are  one  and  one-half  to  two  times  the  horizontal  scale  of 
approximately  20  km  estimated  from  direct  current  measurements  and  from 
hydrographic  data.  Thus,  barocllnlc  Instability  of  the  current-width 
mean  flow  does  not  appear  to  account  for  the  spectral  peak  at  6  days. 
However,  it  Is  quite  possible  that  this  type  of  Instability  does  make  a 
contribution  at  slightly  shorter  periods,  but  Insufficient  to  cause  a 
distinct  spectral  peak.  While  this  may  simply  be  a  result  of  smearing 
of  the  energy  about  a  period  of  6  days  due  to  nonstatlonarlty.  It  Is 
more  likely  that  the  broad  spectral  peak  represents  contributions  from 
more  than  one  source.  This  Is  discussed  In  Section  6.3. 

Calculations  for  case  8,  corresponding  to  a  core-width  length 
scale  of  20  km  and  shallower  bottom  slope  of  0.01  yielded  Instability 
over  a  much  wider  range  of  wavelengths  and  periods  than  did  case  A. 
This  Is  a  result  of  the  lesser  slope.  Recall  that  the  term  Qy  In  the 
stability  condition  equation  (6.2.4)  contains  a  positive  term  6y  which 
corresponds  to  the  restoring  effect  due  to  topography.  Reducing  the 
bottom  slope  reduces  the  value  of  0-p,  and  therefore  lessens  the  con¬ 
straints  on  this  type  of  Instability.  This  effect  Is  reflected  In  the 
dispersion  curves  for  case  B  (Fig.  32b).  For  the  same  smaller  and 
larger  shear  values  as  used  In  case  A,  the  maximum  growth  rate  waves 
respectively  have  wavelength-period  pairs  of  (63  km,  7.5  days)  and  (84 
km,  7.3  days),  and  again  the  difference  Is  due  to  the  smaller  restoring 
force  exerted  by  the  shallower  slope.  The  wavelengths  greatly  exceed 
the  observed  horizontal  length  scales,  but  the  periods  are  only 
slightly  greater  than  6  days.  The  comments  at  the  end  of  the  discus¬ 
sion  of  case  A  as  to  whether  or  not  there  is  a  discernible  contribution 
to  the  spectra  generally  apply  here  also. 

Case  B  Is  essentially  a  hypothetical  case  based  on  an  underlying 
assumption  that  a  velocity  core  scales  the  instability  problem  cor¬ 
rectly.  For  the  WSC,  this  Is  not  necessarily  so,  especially  with  re¬ 
spect  to  the  vertical  shear.  In  general,  the  vertical  shear  has  the 
same  sign  across  the  entire  current  width.  Thus  it  is  probably  more 
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appropriate  to  consider  barocllnlc  Instability  In  terms  of  the  entire 
width  of  the  general  northward  flow  rather  than  of  a  single  core  width. 

To  conclude,  barocllnlc  Instability  of  the  mean  flow  of  the  WSC  Is 
a  possible  mechanism  for  producing  variability  In  the  flow  at  time 
scales  of  3-4  days  and  wavelengths  of  35-40  km.  However,  there  Is  no 
distinct  spectral  Indication  of  the  contribution,  although  any  such 
Indication  may  be  masked  by  unsteadiness  or  other  unexplored  sources. 

In  order  to  proceed  further  In  assessing  the  contribution  by 
barocllnlc  Instability  to  the  variability  of  the  WSC,  It  Is  necessary 
to  complement  model  studies  with  an  energy  budget  approach.  This 
approach  requires  extensive  current  data  coupled  with  hydrographic 
data,  so  that  the  magnitude  and  the  direction  of  the  various  energy 
transfer  terms  can  be  determined.  Such  data  are  not  available  for  the 
WSC,  and  therefore  the  model  results  can  only  be  considered  suggestive. 

6.3  Coastal  trapped  waves 

Coastal  trapped  waves  have  been  advanced  as  an  explanation  for 
low-frequency  variability  of  currents  and  sea  level  In  various  coastal 
regions  of  the  world.  (Nysak,  1980a,b).  These  waves,  trapped  over  the 
continental  shelf  and  slope,  are  another  possible  source  of  velocity 
fluctuations  with  tine  scales  on  the  order  of  a  week.  Two  broad  cate¬ 
gories  of  trapped  waves  are  possible.  In  a  stratified,  constant-depth 
ocean,  one  type  of  free  wave.  Internal  Kelvin  waves,  can  propagate 
along  a  coastal  boundary.  In  a  homogeneous  ocean  with  shelf  topogra¬ 
phy,  a  second  type  of  free  wave,  continental  shelf  waves,  can  also 
exist  and  propagate  along  the  boundary.  Generally,  both  stratification 
and  topography  are  present  and  give  rise  to  hybrid  waves  with  charac¬ 
teristics  of  both  Kelvin  and  shelf  waves.  These  hybrid  waves  are 
called  coastal  trapped  waves  (Gill,  1982).  Details  of  the  properties 
of  these  waves  have  been  studied  by  (among  others)  Allen  (1975),  Wang 
and  Mooers  (1976),  Clarke  (1977),  Huthnance  (1978),  and  more  recently. 
Brink  (1982a,b)  and  Kawabe  (1982). 

Kelvin  waves  owe  their  existence  to  the  presence  of  a  coastal 
wall.  They  stem  from  a  balance  between  the  Coriolis  force  and  the 
pressure  field  with  the  constraint  of  no  normal  flow  at  the  lateral 
boundary.  Thus  In  the  northern  hemisphere,  this  balance  requires  that 
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Kelvin  waves  propagate  with  the  coast  on  the  right,  and  they  decay 
exponentially  away  from  the  coastal  wall. 

Continental  shelf  waves  are  a  form  of  topographic  Rossby  wave,  and 
owe  their  existence  to  both  rotation  and  offshore  variations  in  depth. 
Shelf  waves  conserve  potential  vorticity  and  thus  are  confined  to 
traveling  along  contours  of  f/H,  so  that  trapping  occurs  over  the 
steeply  sloping  topography  of  shelf  and  slope  regions.  Their  motion 
consists  essentially  of  a  series  of  horizontal  eddies  of  alternating 
sign  propagating  along  the  shelf/slope  margin.  They  are  also  right- 
bounded  in  the  northern  hemisphere  (Huthnance,  1978). 

Hysak  (1980a, b)  provides  both  an  extensive  reference  list  and  a 
comprehensive  review  of  observations  of  coastal  trapped  waves  and  of 
the  theories  accociated  with  various  topography/stratification  models. 

6.3.1  Generation  of  shelf  waves 

The  generation  of  continental  shelf  waves  requires  an  input  of 
vorticity  over  the  shelf  region,  and  this  is  usually  accomplished  by 
the  wind.  However,  other  mechanisms  not  related  to  atmospheric  forcing 
such  as  scattering  by  topographic  irregularities,  nonlinear  interac¬ 
tions,  longshore  variations  in  topography,  incident  Rossby  waves,  and 
variable  discharge  at  the  coast  may  also  excite  shelf  waves  (Mysak, 
1980a). 

Adams  and  Buchwald  (1969)  and  Gill  and  Schumann  (1974)  investi¬ 
gated  the  generation  of  shelf  waves  by  the  wind.  Gill  and  Schumann 
(1974)  solved  the  forced  vorticity  equation  for  long,  non-dlvergent 
shelf  waves.  They  found  that  an  oscillating  alongshore  wind  stress  in 
conjunction  with  the  coastal  boundary  provided  the  most  efficient  means 
of  adding  vorticity  over  the  shelf.  Alongshore  wind  stress  In  the 
proper  direction  produces  an  onshore  Ekman  transport  in  the  surface 
layers.  Hass  is  conserved  by  an  offshore  transport  in  the  deeper 
water,  and  It  is  this  deep  flow  over  the  changing  shelf  topography 
which  adds  vorticity  through  vortex  stretching.  If  the  direction  of 
the  wind  stress  reverses,  the  directions  of  the  shallow  and  deep  flows 
reverse,  and  vorticity  of  the  opposite  sign  Is  produced.  Thus  an 
alternating  wind  stress  field  could  produce  vortices  of  alternating 
sign  which  would  propagate  with  the  coast  on  the  right.  Gill  and 
Schumann  (1974)  also  showed  that  vortex  stretching  or  shrinking  due  to 
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either  the  curl  of  the  Mind  stress  (Ekman  pumping)  or  to  changes  in 
atmospheric  pressure  (inverse  barometer  effect)  was  at  least  an  order 
of  magnitude  less  important  than  that  due  to  variations  in  the  long¬ 
shore  wind  stress. 

The  discussion  so  far  has  concerned  the  generation  of  free  waves, 
that  is,  waves  which  are  generated  and  then  propagate  away  from  the 
generating  region.  Gill  and  Schumann  (1974)  showed  that  if  the  gener¬ 
ating  disturbance  moves  along  the  coast,  then  forced  shelf  waves  with 
time  and  length  scales  comparable  to  the  disturbance  may  also  result, 
in  addition  to  the  free  waves.  If  the  disturbance  moves  in  the  direc¬ 
tion  opposite  to  free  wave  propagation,  the  forced  and  free  waves  may 
travel  in  opposite  directions  along  the  shelf. 

For  all  but  very  special  cases,  the  governing  equations  defy 
analytical  solution,  and  so  numerical  methods  have  been  applied  to  the 
problem  of  calculating  the  properties  of  trapped  waves.  Brink  (1982b) 
used  a  numerical  model  to  compute  free  coastal  trapped  wave  frequencies 
and  cross-isobath  modal  structure  for  several  locations  along  the  Peru 
coast,  and  he  compared  his  results  with  direct  observations.  Off  Spits¬ 
bergen,  the  bottom  topography  (Fig.  3),  with  steep  slope  and  shallow 
shelf,  appears  quite  suitable  for  supporting  coastal  trapped  waves.  1 
therefore  examine  the  possibility  of  these  waves  appearing  in  the  WSC 
by  applying  Brink's  (1982b)  model. 

6.3.2  Description  of  the  model 

The  following  presentation  of  the  model  Is  summarized  from  Brink 
(1982b).  The  model  system  Is  linear  and  inviscid  with  constant  rota¬ 
tion,  stratification  Is  horizontally  uniform,  and  the  bottom  topography 
does  not  vary  In  the  longshore  direction.  The  system  of  equations  for 
momentum,  continuity,  and  conservation  of  density  reduce  to  an  equation 
for  the  conservation  of  potential  vorticity: 

Pxxt  +  Pyyt  +  *  8tt) (Pz^)zt  *  (6.3.1) 
subject  to  boundary  conditions: 

pzt  +  (N*/9)Pt  *  0  at  z*0,  (6.3.1a) 
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(f2  +  3tt)pzt  +  N2MPxt  +  fPy)  a  0 

at  z  *  -h{x). 

(6.3.1b) 

pxt  +  fPy  a  0  at  x«0. 

(6.3.1c) 

p  0  as  x  -4.  -  , 

(6.3. Id) 
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where  x,  y,  and  z  are  the  onshore,  alongshore,  and  vertical  coordi¬ 
nates,  p  is  the  perturbation  pressure,  f  is  the  Coriolis  parameter,  g 
is  gravity,  N  is  the  Brunt-Valsala  frequency,  h(x)  is  the  local  depth, 
and  subscripts  denote  differentiation.  The  boundary  conditions 
express:  a)  a  free  surface,  b)  no  normal  flow  at  the  bottom,  c)  no 
normal  flow  at  the  coastal  wall,  and  d)  coastal  trapping.  Searching 
for  free  coastal  waves,  a  propagating  pressure  perturbation  with  wave- 
number  1  and  frequency  u: 

p  «  p(x,z)  exp(iut  +  ily) 

Is  assumed,  so  that  the  system  (6.3)  becomes  a  two-dimensional  eigen¬ 
value  problem,  and  numerical  techniques  can  be  employed  to  solve  it. 
The  problem  Is  solved  on  a  grid  of  25  horizontal  points  by  17  vertical 
points  proportionally  spaced  across  the  breadth  and  local  depth  lim¬ 
its.  Parameterlzatlons  of  stratification  and  bottom  topography  are 
profiles  of  Brunt-Valsala  frequency  squared  and  depth  with  distance 
offshore.  For  a  given  wavenumber,  1,  the  model  seeks  a  corresponding 
frequency,  u,  by  the  technique  of  resonance  Iteration.  That  is,  the 
model  assumes  forcing  at  wavenumber  1  and  Iterates  to  find  u>  for  which 
the  response  Is  a  maximum.  This  (u,l)  pair  specifies  one  point  on  the 
dispersion  curve  for  a  free  wave  of  particular  mode  which  would  be 
supported  by  the  stratification/topography  combination.  The  model  also 
produces  a  matrix  of  the  pressure  distribution  p(x,z)  from  which  the 
alongshore  velocity  distribution  v(x,z)  Is  derived,  and  thus  the 
pressure  and  velocity  structures  over  the  shelf  and  slope  may  be 
mapped. 

The  bottom  profile  used  here  Is  essentially  that  in  Fig.  3,  but 
with  a  flat  shelf.  Two  Brunt-Valsala  profiles,  derived  from  CTD  data 
obtained  In  September  and  November,  1977,  and  differing  mainly  In  the 
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upper  200  ra,  were  used.  Subsequent  computations  showed  that  the 
pressure  structures  obtained  for  both  stratification  profiles  were 
quite  similar,  but  that  phase  speeds  were  about  25%  smaller  using  the 
more  weakly  stratified  November  profile.  This  is  consistent  with 
earlier  findings  that  phase  speeds  of  shelf  waves  tend  to  increase  with 
Increasing  stratification  (Mysak,  1967;  Wang,  1975).  Weak  stratifica¬ 
tion  characterizes  the  region  much  of  the  year,  so  only  the  results 
computed  using  the  November  profile  are  presented  here. 

6.3.3  Results  of  the  model 

Figure  33  shows  dispersion  curves  for  the  first  two  trapped  wave 
diodes.  The  third  mode  was  also  computed,  but  typical  phase  speeds  were 
only  a  few  centimeters  per  second,  and  periods  ranged  from  about  50-200 
days,  so  that  this  mode  is  unlikely  to  contribute  to  variability  with  a 
time  scale  of  a  week  or  so.  The  curve  for  the  first  mode  wave  shows 
that  for  a  period  of  about  3  days,  the  corresponding  wavelength  is 
nearly  2100  km,  and  the  phase  speed  is  about  8.2  m  s*1.  Longer  periods 
correspond  with  longer  wavelengths.  The  first  mode  velocity  structure 
is  essentially  barotropic.  From  a  practical  standpoint,  two  factors 
rule  out  the  applicability  of  the  first  mode:  1)  shelf  waves  are 
usually  generated  by  synoptic  scale  weather  systems  and  therefore 
reflect  the  typical  length  and  time  scales  of  about  1000  km  and  1  week, 
respectively,  and  2)  the  longest  wavelengths  which  might  be  expected 
are  basin-scale  wavelengths  (~1000  km),  but  these  have  periods  slightly 
less  than  1.5  days. 

For  the  second  mode,  a  period  of  14  days  corresponds  to  a 
wavelength  of  about  1000  km,  and  a  period  of  4  days  to  about  260  km. 
In  this  range,  the  second  mode  Is  nearly  non-dispersive  with  northward 
phase  and  group  velocities  both  of  about  75  cm  s'1-.  This  period  range 
corresponds  well  with  the  broad  spectral  peak  in  the  U-component 
spectra,  and  the  wavelength  range  is  compatible  with  the  physical 
setting.  In  addition,  the  decay  time  and  decay  length  scales  inferred 
from  the  model  are  about  10  days  and  1000  km,  so  that  waves  generated 
on  the  east  side  of  the  basin  are  able  to  persist  long  enough  to  affect 
the  flow  off  Spitsbergen. 

In  order  to  make  comparisons  with  observations,  longshore  wavenumber 
information  Is  necessary  to  define  a  trapped  wave  of  given  period  from 
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Figure  33.  Mode  1  and  2  dispersion  curves  for  the  Brink  (1982b) 
coastal  trapped  wave  model  using  stratification  and  topo¬ 
graphy  applicable  to  the  West  Spitsbergen  Current  west  of 
Spitsbergen  at  79°N.  The  non-dimensional  frequency  o  is 
the  dimensional  frequency  <*>  divided  by  f.  The  non-dimen¬ 
sional  wavenumber  k  is  the  dimensional  wavenumber  1  multi¬ 
plied  by  the  shelf/slope  width  L  »  150  km. 
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the  dispersion  curve.  Usually,  records  of  sea  level  height  at  the 
coast  or  current  records  from  moorings  separated  in  the  longshore 
direction  are  used  to  determine  the  longshore  wavenumber  (e.g..  Smith, 
1978).  In  the  absence  of  such  data,  it  is  only  possible  to  guess  at 
the  wavelength  in  an  effort  to  arrive  at  an  associated  wave  period. 
One  approach  is  to  consider  the  most  efficient  forcing  mechanism,  which 
is  the  forcing  of  a  cross  isobath  flow  in  the  deep  layers  by  Ekman 
convergence  or  divergence  at  the  coast.  Approximately  850  km  of  open 
water  separate  Norway  from  Spitsbergen,  while  Spitsbergen  has  a  western 
coastline  about  350  km  long.  Thus,  any  trapped  waves  generated  locally 
near  Spitsbergen  would  have  length  scales  comparable  to  the  length  of 
the  coastline.  From  the  dispersion  curve,  a  second  mode  trapped  wave 
with  a  wavelength  of  363  km  has  a  period  of  5.8  days,  which  is  consist¬ 
ent  with  the  center  of  the  broad  spectra  peak  in  the  cross-isobath 
velocity  spectra.  This  is  not  sufficient  grounds  to  conclude  that 
coastal  trapped  waves  of  the  above  characteristics  account  for  the 
variability  around  a  period  of  about  6  days,  but  the  results  certainly 
suggest  the  possibility  that  they  are  a  feature  of  the  WSC. 

Figure  34  shows  the  velocity  structure  over  the  upper  slope  and 
shelf  computed  by  the  model  for  a  second  mode  trapped  wave  with  wave¬ 
length  500  km  and  period  7.7  days.  This  structure  is  representative  of 
the  frequency-wavenumber  range  of  Fig.  33.  The  velocity  structure  is 
nearly  barotropic,  consistent  with  model  studies  conducted  by  Brink 
(1982b)  and  with  findings  by  Allen  and  Romea  (1980)  which  indicated 
that  with  increasing  latitude,  the  characteristics  of  trapped  waves 
approach  those  of  barotropic  shelf  waves.  In  comparison,  while  there 
is  some  degree  of  shear  in  the  mean  flow  of  the  WSC  (see  Table  5),  the 
current  measurements  show  that  the  fluctuations  tend  to  be  mostly 
barotropic  in  nature.  The  velocity  structure  indicates  a  high  speed 
region  close  to  the  coast,  a  node  over  the  shelf,  a  region  of  high 
velocity  over  the  shelf  break  which  decays  offshore,  and  another  node 
over  the  slope.  This  modal  structure  implies  that  the  variability  due 
to  the  second  mode  should  be  largest  near  the  shelf  break,  decrease 
offshore,  and  then  increase  again  seaward  of  the  node  over  the  slope. 

Spectra  of  alongshore  velocity  show  comparable  energy  content 
between  A1-98V  and  C1-93V  in  the  3-10  day  range.  In  addition,  the 
velocity  fluctuations  of  A1-98V  and  C1-93V  are  negatively  correlated, 
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Figure  34.  Mode  2  perturbation  velocity  structure  for  the  Brink 
♦  (1982b}  coastal  trapped  wave  model  west  of  Spitsbergen 

along  79°N.  Wavelength  and  period  are  500  km  and  7.7  days. 

!  Contours  are  in  arbitrary  units. 
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and  coherence  calculations  between  them  yield  phase  differences  between 
90°  and  180°  for  periods  longer  than  about  3  days.  A  node  between  A 
and  C  would  in  fact  yield  such  a  negative  correlation  and  comparable 
energy  levels.  However,  while  the  mean  alongshore  velocity  at  D  is 
highest  of  all  the  mooring  locations,  the  energy  content  of  its  spec¬ 
trum  shows  levels  comparable  to  those  at  A  and  C  over  similar  time 
scales.  In  the  3-10  day  range,  the  energy  is  comparable  between  A1-98U 
and  C1-93U,  but  both  are  a  decade  greater  than  D1-98U.  Observations 
thus  do  not  bear  out  strong  trapping  of  fluctuating  energy  at  the  shelf 
break,  even  though  the  alongshore  mean  velocity  is  highest  there. 

The  model  results  show  that  the  physical  setting  west  of  Spits¬ 
bergen  can  support  free  coastal  trapped  waves  with  a  period  which, 
assuming  a  reasonable  wavelength,  corresponds  with  a  broad  peak  in  the 
cross-isobath  velocity  spectra.  The  model  predicts  a  nearly  barotropic 
velocity  structure,  which  is  slightly  different  than  shown  by  direct 
velocity  measurements.  The  model  also  shows  a  node  in  the  velocity 
structure  over  the  slope  region,  consonant  with  the  observed  and  calcu¬ 
lated  correlation  between  A  and  C.  In  spite  of  the  highest  mean  veloc¬ 
ity  occurring  at  D,  direct  observations  do  not  bear  out  the  strong 
concentration  of  energetic  fluctuations  over  the  shelf  break  (compared 
with  farther  offshore  as  indicated  by  the  model). 

6.4  Lateral  movement  of  banded  flow 

Historical  evidence  for  banded  flow  in  the  USC  was  presented 
earlier  in  Section  1,  and  observational  evidence  for  banded  flow  from 
the  1976-77  deployments  was  presented  in  Section  3.  Lateral  movement 
of  such  a  flow  back  and  forth  across  the  line  of  moorings  is  one  way  of 
producing  the  negative  correlation  between  currents  at  A  and  C.  Over- 
layed  plots  of  longshore  velocities  A1-98V  and  C1-93V  show  the  negative 
correlation.  With  few  exceptions,  major  fluctuations  in  the  records 
are  coincident  but  opposite  in  sign.  Moveover,  they  are  out-of-phase 
even  when  the  velocity  fluctuations  are  irregular  in  time.  This  last 
observation  suggests  the  lateral  structure  may  be  relatively  invariant, 
in  agreement  with  the  evidence  from  hydrographic  data.  In  addition, 
the  average  phase  difference  between  the  cross-isobath  velocity  fluctu¬ 
ations  of  A1-98U  and  C1-93U  Is  approximately  zero  for  periods  longer 
than  about  3  days.  These  results  support  the  concept  of  a  definite. 


Invariant  flow  structure  meandering  back  and  forth  over  the  slope.  The 
propagation  of  trapped  waves  is  one  way  to  account  for  the  lateral 
movement,  but  not  necessarily  the  sole  explanation,  since  the  fluctua¬ 
tions  do  not  appear  to  be  strictly  periodic  or  of  consistent  time 
scales.  Observations  have  shown  that  the  Kuroshio  may  occupy  different 
paths  (cf.  Stommel  and  Yoshida,  1972)  which  appear  related  to  topog¬ 
raphy,  and  the  flow  of  the  WSC  over  the  slope  may  represent  another 
example. 


6.5  Atmospheric  forcing  and  short  term  variability  of  the 

West  Spitsbergen  Current 

In  Section  5,  it  was  shown  that  annual  variation  in  the  wind  field 
could  account  for  annual  fluctuations  of  the  WSC  if  the  response  were 
in  the  form  of  Rossby  waves.  Higher  frequency  wind  forcing  can  be 
expected  to  generate  higher  frequency  fluctuations  (bounded  by  the 
inertial  frequency).  However,  at  these  shorter  time  scales,  the  nature 
of  the  oceanic  response  will  also  be  different,  as  was  pointed  out  in 
Section  5.  I  now  discuss  the  possible  role  of  atmospheric  forcing  in 
relation  to  fluctuations  of  the  WSC  with  time  scales  of  a  few  days  to  a 
few  weeks.  But  note  that  the  data  needed  to  realistically  assess  the 
degree  of  wind  forcing  on  time  scales  this  short  (at  least  daily 
values  of  local  winds  and  atmospheric  pressure)  are  not  available  for 
the  area  under  study. 

Willebrand  (1978)  found,  in  a  study  of  the  North  Atlantic  wind 
field,  that  uost  of  the  wind  variance  is  associated  with  eastward 
traveling  cyclones  with  time  scales  of  2-10  days.  This  is  consistent 
with  Vowlnckel  and  Orvig's  (1970)  observation  that  4-6  cyclones  per 
month  traverse  the  Greenland  and  Norwegian  seas  in  winter.  Philander 
(1978)  showed  that  eastward  traveling  disturbances  could  not  directly 
excite  free  Rossby  waves.  Instead  the  oceanic  response  is  trapped  in 
the  vicinity  of  the  disturbances,  and  therefore  should  show  similar 
space  and  time  scales.  Philander  (1978)  also  pointed  out  that  the 
forcing  by  these  cyclones,  which  are  characteristic  of  mid-  and  high 
latitudes,  penetrates  to  the  ocean  bottom.  Thus,  bottom  topography  can 
assume  an  important  role  in  the  oceanic  response  to  the  passage  of  the 
disturbances  by  introducing  topographic  waves  with  spatial  or  temporal 


characteristics  which  differ  from  those  of  the  forcing;  the  oceanic 
response  will  reflect  the  characteristics  of  both  these  contributions. 

The  generation  of  shelf  waves  by  Ekman  convergence  or  divergence 
at  a  coastline  has  already  been  discussed.  A  second  possible  source  of 
waves  is  atmospherically  forced  cross-isobath  water  movement  over  the 
section  of  steep  continental  slope  between  Spitsbergen  and  Norway, 
leading  to  the  generation  of  bottom  trapped  topographic  Rossby  waves 
confined  to  the  slope.  The  phase  relation  between  temperature  and 
cross-isobath  velocity  fluctuations  as  calculated  from  C5-984,  where 
temperature  leads  by  about  90°,  is  in  the  proper  sense  to  indicate 
bottom  trapping  (Thompson  and  Luyten,  1976).  However,  the  fluctuating 
energy  is  greater  at  C5-384  than  at  C5-984,  in  disagreement  with  bottom 
trapping,  which  speaks  against  this  explanation  unless  the  C5-384 
record  primarily  represents  variance  originating  in  the  upper  ocean. 
Deeper  measurements,  e.g.,  at  600-700  m,  are  required  to  shed  further 
light  on  the  possibility  of  bottom  trapping. 

The  separation  of  forced  from  free  response  poses  a  problem.  The 
periods  and  wavelengths  of  free  second  mode  coastal  trapped  waves  are 
similar  to  those  of  the  wind  field,  and  the  cross-isobath  velocity 
spectra  peak  broadly  at  about  the  same  time  scale  as  the  winds.  Figure 
35  shows  a  spectrum  for  the  meridional  wind  stress  at  75°N  7.5°E  which 
was  calculated  from  the  daily  NWS  pressure  maps.  The  shape  is  not 
identical  to  the  U-component  spectra,  but  it  is  apparent  that  there  is 
also  a  broad  concentration  of  energy  In  the  wind  field  at  time  scales 
which  correspond  with  fluctuations  of  the  current. 

The  lack  of  appropriate  data  hamper  reaching  any  definite  conclu¬ 
sion  about  the  true  relationship  between  wind  and  flow  of  the  WSC  at 
time  scales  of  a  few  days  to  a  few  weeks.  However,  it  seems  likely  that 
the  mid-frequency  fluctuations  of  the  WSC  reflect  the  complicated 
response  of  the  ocean  to  the  eastward  passage  of  cyclones  across  the 
slope  and  shelf  on  the  eastern  side  of  the  Greenland  Sea. 


7.  Sumroar 


Current  and  temperature  records  have  been  used  to  calculate  the 
volume  and  heat  transports  of  the  West  Spitsbergen  Current,  and  to 
Investigate  the  temporal  and  spatial  variability  of  the  flow.  In  1976- 
77,  four  moorings  were  spaced  laterally  in  the  flow,  and  from  1977  to 
79  and  during  1980-81,  moorings  were  maintained  at  the  same  location  as 
one  of  the  1976-77  moorings  in  the  main  core  of  the  flow. 

Over  the  upper  continental  slope  region,  lateral  spatial  scales  of 
10-20  km  characterized  the  flow,  and  while  the  current  was  in  general 
strongly  correlated  vertically,  there  were  episodes  of  days  to  weeks 
during  which  the  vertical  shear  of  alongshore  velocity  was  large  or 
even  reversed.  Current  meters  separated  by  only  about  15  km  across  the 
flow  were  strongly,  but  negatively  correlated.  Mean  cross-isobath 
velocities  were  offshore  and  generally  about  10-202  of  the  alongshore 
mean.  A  broad  range  of  time  scales  of  variability,  from  tidal  to 
annual,  appears  in  the  current  records. 

During  1976-77,  volume  transports  ranged  from  1.4  Sv  southward  to 
11.9  Sv  northward,  and  the  year-long  mean  was  5.6  Sv.  Of  this,  the 
barocllnic  and  barotroplc  components  were  1.0  and  4.6  Sv,  respectively. 
The  seasonal  pattern  Indicated  is  for  maximum  transport  in  late  fall 
and  winter,  minlumum  in  February-March,  and  reduced  transport  (less 
than  maximum)  during  spring. 

The  mean  heat  transport  during  1976-77  was  46  TW  northward  with 
extrema  of  8  TW  southward  and  122  TW  northward.  The  seasonal  pattern 
was  similar  to  that  of  the  volume  transport.  The  mean  alongshore 
turbulent  heat  transport  was  3.2  TW.  In  terms  of  heat  flux,  the  along¬ 
shore  mean  advective  flux  was  almost  one  order  of  magnitude  greater 
than  the  alongshore  turbulent  flux,  and  at  least  two  orders  of  magni¬ 
tude  greater  than  the  cross-isobath  flux.  The  cross-isobath  turbulent 
heat  flux  was  down-gradient  (offshore),  in  the  proper  direction  to  feed 
energy  from  the  mean  flow  Into  fluctuations. 

Bridged  current  records  three  years  long  clearly  showed  the  strong 
seasonal  fluctuations  of  speed  and  temperature  which  occurred  during 
the  period  1976-79.  There  appeared  to  be  decreasing  trend  in  the 
transport  temperature,  but  no  significant  Interannual  differences  in 
the  volume  transport  over  this  three-year- interval. 
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The  mean  wind  stress  curl  over  the  Green land -Norwegian  sea  region 
appears  to  be  Important  In  driving  the  northward  transport  of  the  West 
Spitsbergen  Current.  During  1976-79,  the  wind  stress  curl  was  about 
one-half  that  In  1971-72,  and  the  mean  velocities  of  the  West  Spitsber¬ 
gen  Current  over  those  respective  periods  were  similarly  related. 
Topography  likely  plays  an  Important  role  in  steering  and  causing 
convergence  in  the  flow.  In  addition,  dynamic  considerations  imply 
that  fluctuations  of  the  wind  stress  curl  can  only  cause  fluctuations 
in  the  mean  flow  at  time  scales  longer  than  about  one-half  year.  Thus, 
the  general  condition  in  the  Green land-Norwegi an  sea  region  seems  to  be 
a  slowly  varying  Sverdrup  balance,  responding  to  low-frequency  fluctua¬ 
tions  of  the  wind  stress  curl. 

The  West  Spitsbergen  Current  exhibits  much  high-frequency  (but 
subtldal)  variability,  and  many  possible  causes  were  investigated. 
Spectra  of  alongshore  current  records  showed  fairly  uniform  energy 
levels  at  periods  longer  than  about  about  four  days.  At  shorter  peri¬ 
ods,  energy  fell  off  sharply  except  for  a  strong  peak  at  the  semi¬ 
diurnal  tidal  frequency.  Spectra  of  cross-isobath  current  records 
peaked  broadly  about  a  period  of  six  days  and  decreased  by  three  orders 
of  magnitude  to  the  annual  period.  Shorter  than  six  days,  cross- 
isobath  and  alongshore  spectra  were  quite  similar. 

Model  results  Indicated  that  barotropic  instability  in  the  West 
Spitsbergen  Current  was  highly  unlikely  because  of  generally  low  Rossby 
numbers  for  the  flow,  and  thus,  barotropic  instability  is  probably  not 
a  source  of  variability  which  appears  in  the  cross-isobath  current 
spectra  at  a  period  of  about  a  week. 

Baroclinic  Instability  Is  one  possible  contributor  to  variability 
at  periods  somewhat  less  than  a  week.  For  parameters  corresponding  to 
the  width  and  speed  of  the  West  Spitsbergen  Current,  model  results 
indicated  that  the  flow  may  exhibit  baroclinic  instability  over  a  range 
of  wavelengths  and  periods  of  30-50  km  and  3-4.5  days,  respectively. 
For  vertical  shear  characteristic  of  the  West  Spitsbergen  Current,  the 
fastest  growing  wave  had  a  wavelength  of  34  km  and  a  period  of  3.8 
days.  Ooubllng  the  shear  nearly  doubled  the  growth  rate,  but  changed 
the  corresponding  wavelength  and  period  by  only  about  20%. 

The  topography  and  stratification  west  of  Spitsbergen  along  79°N 
may  support  a  variety  of  coastal  trapped  waves.  However,  only  the 


second  mode  wave  Is  a  likely  candidate  for  accounting  for  variability 
at  periods  of  about  one  week.  For  periods  longer  than  1.5  days,  first 
mode  wavelengths  are  longer  than  the  basin  scale  length  of  about  1000 
km,  while  third  mode  waves  are  clearly  unimportant.  Corresponding 
wavelength  and  period  ranges  for  the  second  mode  trapped  waves  are 
260*1000  km  and  4*14  days.  This  wavelength  range  Is  physically  reason¬ 
able  for  the  region,  and  the  period  range  corresponds  with  the  range  In 
the  cross-isobath  spectra  In  which  energy  Is  concentrated.  The  second 
mode  trapped  wave  velocity  structure  Is  essentially  barotroplc  over  the 
shelf/slope  region,  and  contains  a  node  over  the  slope.  Both  of  these 
features  are  consistent  with  the  current  measurements.  Alongshore 
wavelength  Information  Is  necessary  to  confidently  assess  the  role  of 
the  second  mode  trapped  wave  concerning  variability  of  the  West  Spits¬ 
bergen  Current. 

In  winter,  4-6  cyclones  per  month  traverse  northeastward  across 
the  region,  and  at  high  latitudes,  forcing  extends  to  the  bottom.  It 
Is  highly  likely  that  these  traveling  disturbances  contribute  to  the 
variability  of  the  West  Spitsbergen  Current  at  periods  corresponding 
with  the  frequency  of  passage.  It  Is  also  likely  that  the  oceanic 
response  Is  modified  by  the  Influence  of  strong  topographic  features. 
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8.  Future  Investigations 

There  are  two  different  directions  in  which  to  progress,  and  the 
data  requirements  are  somewhat  different  for  each.  One  Is  the  investi¬ 
gation  of  the  short-term  variability  and  associated  dynamics  of  the 
WSC,  and  the  other  is  the  maintenance  of  a  long-term  monitoring  program 
to  measure  the  transports  into  the  Polar  Basin  in  order  to  assess  the 
Implications  for  the  heat  budget. 

Needed  for  further  Investigation  of  the  short-term  variability  are 
alongshore  wavelength  information,  local  wind  stress  data,  and  sea 
surface  height  records  at  some  point  on  the  west  coast  of  West  Spits¬ 
bergen  In  conjunction  with  current  measurements.  Ideally,  periodic 
hydrographic  sections  of  closely  spaced  stations  would  complement  the 
measurements.  These  types  of  data  are  necessary  for  differentiating 
between  possible  dynamic  sources  of  variability  In  the  flow  and  for 
determining  the  Importance  of  local  winds  in  generating  that  variabil¬ 
ity.  All  are  obtainable  with  the  existing  technology  In  oceanographic 
Instrumentation.  There  ere  many  existing  studies  of  boundary  currents, 
and  further  detailed  study  of  the  dynamics  of  the  short-term  variabil¬ 
ity  of  the  WSC  presents  still  another  unique  addition  to  that  branch  of 
physical  oceanography. 

One  of  the  key  factors  In  the  success  of  a  long-term  monitoring 
program  In  the  WSC  Is  refining  the  transport  estimates  of  the  flow  and 
linking  the  transport  to  some  Index  such  as  the  flow  at  a  single 
location.  Refining  the  estimates  requires  an  accurate  characterization 
of  the  lateral  and  vertical  structure  of  the  flow,  and  based  on  the 
present  work,  this  would  require  an  array  of  moorings  about  10  km  apart 
spanning  from  the  coast  of  Spitsbergen  westward  to  about  3°E,  with  at 
least  some  deep  current  meters  to  guage  the  bottom  flow.  Linking  the 
transport  to  a  single  Index  Is  simply  a  matter  of  correlation  calcula¬ 
tions.  The  present  work  represents  a  step  in  this  direction,  and  is 
probably  sufficient  for  making  coarse  estimates  of  long-term  transport 
fluctuations.  However,  It  Is  Important  that  future  experimental  work 
take  Into  account  the  lateral  structure. 

Such  an  extensive  program  of  measurements  Is  an  expensive  and 
formidable  undertaking,  and  the  use  of  remote  sensing  techniques  is  one 
way  of  reducing  the  amount  of  ground-based  data  required.  The  strong 


temperature  signature  of  the  WSC  and  the  variation  of  dynamic  height 
across  the  flow  suggest  that  remote  sensing  of  surface  temperature  and 
sea  surface  height  is  a  possible  way  of  defining  the  lateral  and  along¬ 
shore  structure  of  the  flow,  particularly  if  it  is  possible  to  obtain 
ground  truth  Information  by  simultaneous  earth-based  and  space-based 
monitoring. 

The  work  reported  in  this  investigation  made  limited  use  of 
models.  As  better  data  become  available  for  Inserting  into  model 
equations,  models  will  become  better  descriptors  of  the  flow.  As  such, 
models  have  valuable  potential  which  should  not  be  neglected  in  opti¬ 
mizing  the  design  of  field  experiments. 

The  WSC  is  rich  in  both  spatial  and  temporal  variability,  and  it 
has  some  consequence  for  the  budgets  of  the  Polar  Basin.  Future  inves¬ 
tigations  thus  have  the  great  potential  to  be  fruitful,  interesting, 
and  Important. 
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